INTER MITWIRKUNG VON L. COLLATZ - H. GORTLER- K. KLOTTER - H. NEUBER 
V. TOLLMIEN UND C. WEBER » HERAUSGEGEBEN VON FR. A. WILLERS, DRESDEN 


nn — — — —— — —  —— 


AND 37 SEITE 1—80 HEFT 1/2 JANUAR/FEBRUAR 


AUS DEM INHALT: 


"R. Zurmühl, Behandlung der Plattenaufgabe 
nach dem verbesserten Differenzenverfahren / 
EN P. Beckmann, W.Franz, Asymptotisches Verhalten 
er der Zylinderfunktionen in Abhängigkeit vom kom- 
| : plexen Index / A. M. Sen Gupta,Stresses in certain 
thin elastic plates rotating about normal axes / 
G.R. Verma, ApplicationofDirac’s delta-function 
in isolated force problems ofsemi-infinite elasticsolid 
of isotropicand non-isotropic materials /[E.Koppe, 
Die dicke Platte mit nichilinearer Spannungsvertei- 
lung / V. Vodicka, Biegungsschwingungen in zu- 
sammengesetzten Stäben / J. Münch, Beiträge zum 
Quellsenkenverfahren für die Berechnung von 
Überschallströmungen/C. Heinz,Reflexion ebener 
Druckwellen an einer festen Wand 


KLEE EM TIERE ER UN SEES 
B-U=C -H BES: P.R-E GH EINES FEEN 
EINGEGANGENE BÜCHER 
NEFAT-C: HIER: 1 Ser 


| 


BR ADEMITE.-VERLAG GMBH - BERLIN 


ee EEE 


Z. angew. Math. Mech. Bd. 37 Nr. 1/2 S. 1-80 Berlin, Januar/Februar 1957 


: - Nachrichten BE ne Dee eines ey es 80 


2 \ 
% .. 0 Haren 2 


3 Eerannt w. ne Asymptotisches Verhalten de 
der Zylinderfunktionen in Abhängigkeit vom komplexen “x 5 : 


Index... ad le ner > aa a au 


A.M.Sen Gupta: Stresses in certain thin elastic plates 
a JE3 aboMt narmal.axes. a RE AN SEEERT ARE 


. R. Verma: Application of Dirac’s delta-funktion 
in isolated force problems of semi-infinite elastic solid of 
'isotropic and non isotropie materials... .. 34 


schaften zu Berlin 


'Koppe: Die dicke Platte mit nichtlinearer a Ki A Bu er Ba 
Spannungsverteilung enenreheenneenunene 38 Erhard Schmidt ee 


v Voditka: Biegungsschwingungen in zusammen- 
.  geseizten Stäben 220er erieeee nenn Al 


Er Münch: Beiträge zum Quellsenkenverfahren für die 
Berechnung von Überschallströmungen a  SSLC 


ec Heinz: Reflexion ebener Druckwellen an einer festen 63 
Wand un... re! Re ER. 


Georg Hamel, Helmut Hasse, H.L. Schmid 
und Kurt Schröder 


> 


AUS DEM INHALT DES HEFTES 5/6 Ba. 5 


‘ Othmar Zoubäle 


Kleine Mitteilungen. W. Uhlmann: Zur = £ EEE € SE 
 Fehlerabschätzung bei Interpolationspolynomen..... . . 73 Über ein SEELE EEE für ge 
3x Br Everling: Eine Verallgemeinerung des Horner- ae erfand kreis me a a et 


‚schen SCHEMASCHST TE eo Behand EEE 74 


. A. Aiserman u. F.R. Eontmechefi Stabilität - 
der Gleichgewichtslage in einem nicht-holonomen 


rungen 


Josef Weier 
Über offene Euklidische Mengen 


# 


% Be rreeningen en eTrÄR an T lee 75 ne 


_ Eingegangene Bücher. MN IN 79 ge ee a a Tee 


mationen im projektiven R, 
 ANSCHRIFTEN DER AUTOREN DIESES HEFTES 


Dr.-Ing. R.Zurmühl, Darmstadt, Martinstr. 75 


SDE-P. a E 
RESEAN Franz, Münster (Westfalen), Schloßplatz 2 


Prof. A. M. Sen Gupta, P. ©. Mesra, Ranchi (Bihar-India) 
Prof. G. R. Verma, Pilani, Rajasthan (India), Birla College of 
Science 
E. Koppe, Göttingen, Max-Planck-Institut für Strömungsfor- 
schung, Böfttingerstr. 6/8 
Dr. V. Vodicka, Pilsen, Stalinstr. 32 
Dr. J. Münch, München 2, Walter-van-Dyck-Platz 1 
Dr. C. Heinz, Lörrach, Haagner Str. 32a 
Dr. W. Uhlmann, Hamburg 13, Harvestehuderweg 10 
Dr. W, Everling, Aachen, Templergraben 55, Math. Institut 
Lehrstuhl C 


Prof. M. A. Aisermann, Ya Institut für Automatik und 


Lothar Berg 


Über die Resolvente einer singulären Integral- 
- gleichung 


Wilhelm Specht 
Die Lage der Nullstellen eines Polynoms 


Rolf Reissig 


Über die Beschränktheit der Lösungen einer 
“ nichtlinearen Differentialgleichung 


Bestellungen durch eine Buchhandlung erbeten 


AKADEMIE-VERLAG-BERLIN 


Telemechanik der Akade- 
Prof. F.R. Gantmacher, Moskau mie d. Wiss. der UdSSR 


Prof. Dr. H. Schubert, Halle (Saale), Am Krähenberg 5 
Prof. Dr. W. Lichtenheldt, Dresden A 27, Einsteinstr. Ta 
Prof. Dr. L. Schiller, Weilburg (Lahn), Frankfurter Str. 19 | 
Dipl.-Ing. A. Schubert, Dresden A 20, Thomas-Mann-Str. 24 Wir bitten alle Manuskriptsendungen direkt an den Her- 
Prof. Dr. H. Heinrich, Dresden A 27, Friedrich-Hegel-Str. 31 ausgeber Prof. Dr. Dr. h. c. Fr. A. Willers, Dresden A 36, 
Prof. Dr. W. Albring, Dresden A 20, Südhöhe 9 Kauschaer Str. 6, zu schicken. Zu jeder Arbeit erbitten wir 
Dipl.-Math. K.-H. Bachmann, Dresden N 23, eine kurze Zusammenfassung des Inhalts, die von uns in 

Wilder-Mann-Str. 42 b die englische, französische und russische Sprache übersetzt 
Prof. P. Görlich, Jena, Zeiss-Werke wird. DerVerlag behält sich für alle Beiträge das Recht der 
Prof. Dr. Fr. A. Willers, Dresden A 36, Kauschaer Str. 6 Vervielfältigung, Verbreitung und Übersetzung vor. 


Prof. Dr. G. Bürgermeister, Dresden A 20, Tiergartenstr. 86 R R 
Obertäg. N. Klamka, Berlin-Karlshorst, Hentigstr. 25a Die Autoren erhalten von den Hauptaufsätzen 75, von den 


Dr. K. Magnus, Freiburg (Breisgau), Unter Linden 7 Kleinen Mitteilungen 10 Sonderdrucke ohne Berechnung. | 


ee WW 


ar 
u 


5 ® hr SR 2 EN ST I u 


4 P3 a T 


{713 N ANILNG- E Bu Ne DT DE 


THEMATIK UND MECHANIK 


INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 


Band 37 Januar/Februar Anl: Heft 1/2 


Behandlung der Plattenaufgabe nach dem verbesserten 
| Differenzenverfahren er 


Von Rudolf Zurmühl in Darmstadt 


für die Plattenaufgabe Adu = p/K mit rechteckigen Rändern werden Differenzenformeln höherer 
Genauigkeit nach dem Mehrstellenverfahren entwickelt, die sowohl das Einarbeiten der Randbedingungen als 
auch die Berechnung der zweiten Ableitungen, also der Biegemomente mit hoher Genauigkeit ermöglichen. Die 
Verwendung der Formeln wird für die verschiedenen Lagerungsfälle angegeben und an einigen Zahlenbeispielen 
erläutert. Kin Vergleich mit bekannten strengen Lösungen zeigt die Güte des Verfahrens. 


Assuming the plate to be rectangularly bounded, the plate problem AAu = p/K is treated by means of an 
improved caleulus of finite differences which permits both the treatment of the boundary conditions and the 
determination of the 2nd derivatives, ü. e., the endbing moments, with high accuracy. The application of the 
formulae to the various possible ways of supporting the plate is discussed and illustrated by numerical examples. 
A comparison with known exact solutions demonstrates the accuracy of the method. 


Pour le probleme de plaque AAu = p/K avec des bords rectangulaires des formules de differences d’une 
ewactitude superieure sont developpees d’apres la methode des plusieurs lieux (Mehrstellenverfahren) permettant 
la consideration des conditions aux bords ainsi que la calculation de la deuxieme derivee, donc la calculation 
tres exacte des moments de flexion. L’application des formules est indiquee pour les situations differentes et 
elle est illusiree au moyen de quelques exemples numeriques. Une comparaison avec des solutions connues et 
rigoureuses demontre la qualit& du procede. 


B z3apnaue umaru6a maactuu AAdu=p/K C NPAMOYTOAbBHbIM KOHTYyPOM BEIBONATCH 
PA3HOCTbHbIeE COOTHOMIEHHSI IOBBIIIEHHOH TOYHOCTH, KOTOPbIe TO3BOANAMT YHOBJIETBOPHTL 
KpaeBbIM YCIOBHAM, A TAK5KE BbIJHCIHIHTBb BTOPbIE IPONH3BONHBIE MH, CJIEIOBATEIIBHO, UH3TH6AIOIIME 
MOMEHTbI € 60AbIH0Ä TOYHOCTbIO. IlpuMmeHeHHe hPopMyA YRAa3bIBAaTca NAH Pa3ımYHBbIX CIIOCO6OB 
KPelllleHUA IIIACTUHKH H UIIIIOCTPUPYETCA YUCJIOBBIMH IPHMepamu. CpaBHeHHe C H3BECTHBIMH 
TOYHbIMH PeIleHMAMM TIOKa3bIBaeT IIPeuMmyecTBa 3T0oTO MeTona. 


1. Einleitung 


Für die Plattenaufgabe, also die Berechnung von Verformungen und Beanspruchungen einer 
dünnen elastischen Platte unter senkrechter Belastung, sind strenge Lösungen in Form von 


_ Reihenentwicklungen nur in begrenzter Anzahl verhältnismäßig einfacher Lagerungsarten be- 


kannt. In allen übrigen Fällen ist man auf Näherungslösungen angewiesen. Nun stellt die Auf- 
gabe an die bekannten Näherungsverfahren der praktischen Mathematik ungewöhnlich hohe 
Anforderungen, indem hier in der Hauptsache nicht die Verformung, also die Lösung der Rand- 
wertaufgabe selbst interessiert, sondern ihre zweiten partiellen Ableitungen, denen die Biege- 
beanspruchungen proportional sind. Das gleiche trifft auch auf die in vielem verwandte Scheiben- 
aufgabe zu, auf die sich die hier entwickelte Berechnungsmethode gleichfalls anwenden läßt. Der 
besagte Umstand hat zur Folge, daß Verfahren, die sonst gute Näherungswerte liefern, für die 
Plattenaufgabe versagen können, indem sie zwar die Funktionswerte gut, ihre zweiten Ablei- 
tungen dagegen nur unzureichend wiedergeben. Das gilt vor allem für die Verfahren, bei denen die 
Abweichungen der angenäherten von den strengen Funktionswerten möglichst klein gemacht 
werden, wie die Fehlerquadratmethode und die Kollokationsmethode!). Aber auch das bekannte 
Ritzsche Verfahren arbeitet nicht mehr befriedigend. Damit kommt dem auf anderem Prinzip 
beruhenden Differenzenverfahren besondere Bedeutung zu, bei dem die in der Differentialgleichung 
auftretenden Differentialquotienten durch Differenzenquotienten ersetzt werden. Es ist daher 
auch schon frühzeitig für die Plattenaufgabe herangezogen worden, und zwar vor allem in der 
bekannten Schrift von Marcus [2], wo die Differentialgleichung 4. Ordnung in zwei Gleichungen 
2. Ordnung aufgespalten und diese nach dem Differenzenverfahren bearbeitet sind. Indessen 
sind die Ergebnisse bei nicht zu enger Maschenweite, also erträglichem Arbeitsaufwand nur mäßig 
genau, und sie befriedigen auch nur bei verhältnismäßig einfacher Lagerungsart. Bei verwickelter 
Lagerung hat Marcus das Verfahren mit Reihenentwicklungen kombiniert. 

Nun sind auf dem Gebiete des Differenzenverfahrens in den letzten Jahren wesentliche Ver- 
besserungen erzielt worden, insbesondere in der Richtung des sogenannten Mehrstellenver- 


1) Vgl. etwa [4], S. 130/31, 148/51, 184/51, 184/87; [6], 8. 428. 
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Steifigkeit K unter senkrechter Last der Dichte p = p(, y) lautet die Differential 


zue ı oliene Kragen zu Kla 

"bedingungen und die Berechnung der zweiten Ableitun, 

_ friedigenden Lösung gebracht werden. Ich beschränke mich 

eckigen Stücken zusammengesetzte Berandung, schließe also 
 berandete Platten us. RT ?. 


n zu 


ne ee, Die Plattenaufgabe ee 
Für die Durchsenkung u = u(x, y) der Mittelfläche einer dünnen elast 


AAu = 4, play) u. ae 


© | ee. 
mit dem Operator A = FIR + ae’ also AAu = Un + 2 Ugayy + Uyyyy Aus der Ireloru SEE 
leiten sich die eigentlich interessierenden Schnittgrößen, nämlich die Biegemomente m,, My 
Drillmoment m,, und Querkräfte q,, q,, bezogen auf die Längeneinheit, folgendermaßen her: 


m, ee en 
m, = —K (Uyy + HU) 
Er RU 17) 12 NEE ni 

Ge = HK Una + I OR REEL Er 
I = —K (ug + Uyaz) 7 


Außer den Querkräften interessieren noch die an Rändern x = konst. bzw. y = konst. angenome- 
nen sogenannten Randkräfte 


; Ir =—K ur +@ — u) Uzyy) ’ (5), 

%y=—K U + 2 — u) Un) 
die aus den Querkräften durch Zuschlag sogenannter Ersatzscherkräfte hervorgehen. e 
Zur Differentialgleichung treten Randbedingungen. Wir nehmen zueinander senkrechte 


Randlinien an, die in Richtung x = konst. und y = konst. verlaufen. Für den Rand x = konst. 
kommen folgende Randbedingungen in Betracht: | 


1. eingespannter Plattenrand u=0, „=0 as LS 
2. frei-aufliegender Plattenrand u=0, m, =0, also u, = U IE EEE: 
3. freier Plattenrand m,=0, ao’. Fuu,=0 (8) 
,.=0, also u. +2 W)uy=0l ' . 

Dazu kommt noch für die freie rechtwinklige Plattenecke: 
My 0; "also use ON N 


Vielfach wird die in den Randbedingungen (8) auftretende Querdehnungszahl u vereinfachend zu 
Null angenommen. Wir werden hier alle Betrachtungen für beliebigen u-Werte durchführen. 
Erst in den Zahlenbeispielen in Nr. 9 ist durchweg u = 0 gesetzt worden. 

Für die Lagerungsarten bedienen wir uns folgender Symbole: 


eingespannter Plattenrand I.) 
frei aufliegender Plattenrand u.=0% 
een freier Plattenrand m,=0, 9,=0 


3. Das Mehrstellenverfahren der Plattengleichung 
Zur Behandlung unserer Aufgabe nach dem — gewöhnlichen oder verbesserten — Differen- 
zenverfahren wird die — z. B. rechteckig berandete — Platte mit einem quadratischen oder recht- 
eckigen Netz überzogen. Wir verwenden im folgenden ausschließlich quadratische Netze mit 


Be 
nahe ok dem Ditforen venvorfa e 


2 ei 


sich der Plattenrand darin nicht ein, so kann man sich 


E) s elfen Fa man zwei ee Fälle durchrechnet, auf die das quadratische Netz paßt, 
un die Ergebnisse auf den dazwischen liegenden Fall interpoliert. Betrachtet werden dann die j 
= ‚Funktionswerte u(&;, ;) bzw. deren Näherungswerte u,, in den Gitterpunkten x, y, des Netzes. 3 


-/ 


Beim gewöhnlichen Differenzenverfahren ersetzt man die in der Differentialgleichung auf- 
retendlen Ableitungen durch Differenzenquotienten der Näherungswerte ı,,. Indem man dies e 

 anallen inneren Punkten und gegebenenfalls auch — unter Einarbeitung von een — a 
an Randpunkten durchführt, erhält man bei linearer Differentialgleichung ein lineares Gleichungs- x 
system für die Unbekannten u,,, das — etwa durch Elimination nach dem Gaußschen Algo- La, 
rithmus?) — gelöst wird. nr 
: Beim sogenannten Mehrstellenverfahren wird der Umstand ausgenutzt, daß von der Ba 
Lösungsfunktion u(z, y) nicht allein diskrete Werte u,, vorliegen, sondern über die Differential- A 
_ gleichung — wenigstens bei hinreichend einfachem Bau dieser Gleichung — an jeder Gitterstelle ea 

auch noch eine bestimmte Wölbungseigenschaft in Form von Ableitungen, bei der Plattenauf- N 
gabe die Werte AAu;, = P(&;, Y)[K = p;rK. Indem man diese Werte nicht nur an dereinen 

gerade betrachteten Stelle verwendet, sondern auch noch an mehreren dem Gitterpunkt i,k % 

benachbarten Punkten r, s, erhält man genauere Formeln der Bauart 


Li 


Drau Jedi Ay AAu,, + höhere Taylorglieder N, 2.202700 


Dabei erstrecken sich die Summen über gewisse Punkte des Gitters, über je einen „Stern“ für die 
Werte u,, und AAu,,. Die Koeffizienten a,, und A,, werden so bestimmt, daß bei Taylorent- 
wicklung der Werte ı,, und AAu,, an der Stelle i, k die Gleichung in möglichst hohen Taylor- 
gliedern erfüllt wird. Wieder erhält man, indem man die Mehrstellenformel sowie ähnlich gebaute 
3 Randformeln zur Berücksichtigung der Randbedingungen auf alle in Betracht kommenden 
= Gitterstellen der Platte anwendet und die Differentialgleichung in der Form AAu,;,; = p;,|K 
_ heranzieht, ein lineares Gleichungssystem für die Unbekannten u;;. 


4 Die Mehrstellenformel der Plattengleichung lautet bei quadratischem Gitter der Maschen- 


} weite A in leicht verständlicher Schreibweise: 

2 en Ba 
E VB U Brot 1 1 

2 [4] 1-0 3 —0 I|u=5|1 2 1|4Munı a. 
192° 0 —2 1 = 1 

r 1 1 1 


| Diese im folgenden als Hauptformel bezeichnete Näherungsgleichung ist von beträcht- 

; licher Genauigkeit; ihr Fehler ist von der Ordnung A®. Sie besagt, daß der Wert u,, des betrachteten 

y Gitterpunktes i, k den Koeffizienten 36, die 4 rechts und links, oben und unten benachbarten 
Werte ı;;1,%, U,xıı je den Koeffizienten —10 erhalten usf. Auf der rechten Seite erhält 'der 
Wert AAu,, h? den Koeffizienten 1, die benachbarten Werte je den Koeffizienten 1/2. Der 
Bau der Formel wird durch einen ‚Stern‘ nach Bild 1 veranschaulicht, wo jeder Punkt einen 
Funktionswert u, jeder Kreis einen Wert AAu andeutet. 


Bild 1. Stern der Hauptformel A Bild 2. Innere Punkte O, randnahe Punkte x und Randpunkte » 
einer rechteckig berandeten Platte 


Durch Formel A sind in einer rechteckig berandeten Platte alle im Abstand 2 h und mehr 
vom Plattenrand gelegenen Gitterpunkte als sogenannteinnere Punkte erfaßbar, Bild 2. Nicht 
erfaßt werden außer den eigentlichen Randpunkten noch die vom Rand um A entfernten rand- 


2) Vgl. etwa [6], 8. 91 ff. 
1 * 


4 _  Zurmühl, Behandlun 


nahen Punkte. Für diese letzten bei 


‚wird. Die sehr genaue Hauptformel A des Mehrstellenverfahrens ist nur sinnvoll in Ver 


5 Be 


entwickeln, welche zugleich die Randbe \ bı 
Dies führt somit zu der ersten der beide: folgende ı Aufgaben, mi 
Arbeit befassen will: nn AT 
1. Entwicklung genügend genauer Randformeln, d. h. Näherungsformeln, v 
bedingungen mit gleicher oder doch fast gleicher Genauigkeit wiedergeben, wie sie ( 
formel aufweist; : TE >27 
2. Entwicklung entsprechend genauer Ableitungsformeln, also Näherungsformeln zur 
rechnung der als Biegemomente interessierenden zweiten Ableitungen u,,, U,, aus den gew« 
nen Näherungswerten u;,. I A N ib: 
Mit der befriedigenden Lösung beider Aufgaben steht und fällt eine sinnvolle 
des genaueren Differenzenverfahrens auf die Plattenaufgabe. Denn ein Blick auf Bi 
stark die Anzahl der Randpunkte und randnahen Punkte die der inneren überwiegen H 
so mehr, je gröber das Netz gewählt wird, was mit dem genaueren Verfahren ja gerade 


a 


mit Randformeln entsprechender Genauigkeit. Die Genauigkeit beider Formeln muß aufein: 
abgestimmt sein. Das Gleiche gilt für Formeln zur Berechnung der zweiten Ableitungen. Ihre 2 


Berechnung als gewöhnliche Differenzenquotienten würde bei der im Mehrstellenverfahren mög- 


lichen groben Gittereinteilung zu ganz schlechten Werten führen, womit der Genauigkeitsgewinn ö 
des Mehrstellenverfahrens wieder völlig verloren gehen würde. r. 

Beide Aufgaben, die Entwicklung von Randformeln und von Ableitungsformeln, sind im 
Grunde ein und dieselbe, nämlich Aufstellung von Näherungsformeln für bestimmte Ableitungs- 
größen, sei es am Rande oder im Innern des Gebietes. Beide lassen sich nach dem Mehrstellen- 
prinzip, also durch Heranziehen der durch die Differentialgleichung gegebenen Wölbungsgrößen 
AAu an mehreren Stellen, in höchst befriedigender Weise lösen, wie die Ergebnisse der Zahlen- 
rechnungen am Schluß der Arbeit bestätigen werden. 

Rand- und Ableitungsformeln haben somit auch einen ähnlichen Aufbau wie ihn die Haupt-- 
formel in (10) aufweist. Bezeichnet allgemein D[u] die fragliche Ableitungsgröße, z. B. u, oder u,, 
oder dergleichen, so lautet hier der Ansatz 


Dful;,; = I a,u,, + 3% A,, AAu,, + höhere Taylorglieder ... . . . . (11), 


wobei sich die Summen wieder über passend zu wählende Sterne erstrecken. 

Wieder werden die Koeffizienten a,,, A,, so bestimmt, daß möglichst weitgehender Taylor- 
Abgleich erzielt wird. Es zeigt sich, daß die Genauigkeit — bei entsprechendem Aufwand wie bei 
Formel A — infolge verlorener Symmetrie gegenüber der Hauptformel nur um eine A-Potenz 
geringer wird, also von der Ordnung A® ist. 


4. Randformeln 


Bei Aufstellen der Randformeln in der Form (11) zeigt sich zunächst, daß der Abgleich der 
Taylor-Glieder durch Funktionswerte u,, ganz bestimmter um die Randstelle gelegener Gitter- 
punkte erreicht wird, soweit die erste Summe in (11) in Betracht kommt. Die abgeglichenen 
h-Potenzen verteilen sich auf die am linken Rand gelegenen Gitterpunkte nach folgendem Schema, 
wobei Symmetrie bezüglich der waagerechten Mittellinie besteht: 


WANN 6 
Ye Wer eye 
0 Mae ae 
a N 
His bieß 
6 


Die mit jeder h-Potenz hinzutretende Punktzahl entspricht dabei gerade der Anzahl der 
abzugleichenden Glieder, so daß auf diese Weise, allein mit der ersten Summe von (11), Rand- 


®) In der oben zitierten Arbeit von Collatz [5] wird noch eine außerhalb der Pla 
Punktreihe mitgenommen. Die inneren randnahen Punkte sind dann der Hauptformel er 
äußeren werden durch die Randbedingungen erfaßt. Dies ist dort einfach möglich, weil nur die beiden Randbedin- 
gungen 4, — 0 und 4, = 0 betrachtet und auch nur durch einfache Differenzenquotienten angenähert werden 
womit die Außenpunkte praktisch nicht als zusätzliche Unbekannte auftreten, sondern zu den inneren in ganz 
einfacher Beziehung stehen. Sollen aber die Randbedingungen genauer erfaßt werden, so wird durch en 
were a en Sn a ge sondern nur rechnerische Bedeutung zukommt, die 
er Unbekannten und damit der eitsbedarf ga; ächtli Ö i hi 
von der Einführung äußerer Punkte konsequent aha Perärhtich nlöbn Au 


u er eieeee 


:npri es Anden en Far Porec en höher SERIEN 
g des „Sternes“ zu gewinnen. Wie weit man auf diese Weise Gitter- 

nkte einsparen, also u-Werte durch AAu-Werte ersetzen kann, ist von Fall zu Fall zu unter- 

ı1en. 

Er , Hier und im folgenden soll durchweg gefordert werden, daß die Ausdehnung des Sternes 
2 = nicht über 4 h nach beiden Richtungen hinaus geht, wie oben im Schema a Dies ist 
_ mit Rücksicht darauf erwünscht, daß man für eine Platte be- N. 
stimmter Ausdehnung andernfalls, bei breiter gebauten For- An 
meln, zu feinerer Gittereinteilung übergehen müßte, womit sich 
der Arbeitsaufwand immer gleich um ein Vielfaches erhöht. 
Denken wir etwa an eine quadratische Platte der Seitenlänge a 
mit Randwerten u=0 und beliebiger Belastung, so sind bei 
einer Maschenweite h = a/4 insgesamt 9 Unbekannte u,, ent- 
sprechend den 9inneren Punkten zu berechnen, bei ha a/d 
_ hingegen 16 Unbekannte, was einer Erhöhung der Eliminations- 
arbeit auf etwa das 6-fache gleich kommt. Es ist also von größ- 


ter praktischer Wichtigkeit, mit möglichst wenig Gitterpunkten re 

auszukommen. +92 
3 Im vorliegenden Falle gelingt die Einsparung der über die = Banane ar % 
4 h-Begrenzung hinausreichenden Gitterpunkte — drei unter i 

Berücksichtigung der Symmetrie — durch Aufnahme entsprechend vieler AAu-Werte, wobei 


ein Abgleich bis h® erreicht wird. Durch Hinzunahme eines vierten AAu-Wertes läßt sich so- 
gar noch ein weiteres Gitterpunktpaar einsparen, und man erhält für die im folgenden auf- 
geführten Randformeln einen Stern nach Bild 3, wo wieder die in der Formel auftretenden 
 AAu-Werte durch Kreise gekennzeichnet sind. 


Gebraucht werden folgende Randformeln für einen Rand x = konst.: 


N ER 


2 8 2 
16 —80 16 P) 
> Rı] 60u,h=| —_149 [336] _ 14 32 —31u—-10.8 0|AAuht-+ (AM), 
: —— 2 
E- 16 —80 16 
E 2 8 2 


t — 180 ugy h? = 


13 
w—|0 124 —72 |AAuht + (h?), 
ai 
—l 
uw—|1 0 7|AAuM + (A), 
—l 
3 
u—-|6 24 —18 |44u ht — 
3 
3 
u-+u- 3 —12 21|AAuf + (A?) . 
3 


[R3] 247, W = 
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Durch diese Formeln sind die Randpunkte und randnahen Punkte erfaßbar, soweit sie um 
mindestens 2 h von einer Ecke entfernt liegen. Der Punkt, auf den die Formel angewandt wird, 
ist in ihr durch Umrahmen hervorgehoben. Im einzenen werden die Formeln bei den verschie- 


denen Lagerungen folgendermaßen angewandt: 
Rand | Randbedingung | Formel | Angewandt auf 


Eingespannt 


Ur = 0 Bi Randnahen Innenpunkt 


Frei aufliegend 


ee) R2 Randnahen Innenpunkt 
K 
Frei 
Uge + MUyy — 0 R2, R2a | Randnahen Innenpunkt 
o0—e Ugez + (2 — U) Uryy — 0 R3 Randpunkt 


5. Eekformeln 


Besondere Behandlung verlangen die Plattenecken. Die hier gelegenen vier Punkte Nr. 1 
bis 4 der Abb. 4 sind durch die bisher aufgeführten Randformeln nicht erfaßbar, es werden also 
besondere Eckformeln benötigt, die teils symmetrisch, teils nicht symmetrisch gebaut sind. Beim 
Abgleich der Taylor-Glieder nach einem Ansatz (11) zeigt sich wieder eine besondere Verteilung 
der abgleichbaren h-Potenzen auf die Gitterpunkte nach untenstehendem Schema, soweit die 
erste Summe in (11) in Betracht kommt. 


SOHmmwe | 5 
wm oIlo 


IS ger) 


2 4 
Bild 4, Lage der Eckpunkte 


Soll die Eckformel nun wieder auf eine Ausdehnung von 4 hin beiden K i l 
beschränkt bleiben, so zeigt die Rechnung, daß dann en Abgleich bis }® 
schen Falle ‚möglich ist, nicht dagegen im allgemeinen nicht symmetrischen, bei dem noch 
wenigstens ein Punkt in der Entfernung 5 h von der Ecke aus in x- oder y-Richtung hinzugenom- 
a Be muß. Um doch die Formelausdehnung von 4 h bei einer Genauigkeitsordnung von 
. ee Au können, haben wir uns bei den nicht symmetrischen Formeln El und E2 für 
ne ne u, =D bzw. u,, = 0 eines Kunstgriffs bedient, der darin besteht, daß wir in 
ee ie „triviale Randbedingung U,, = 0 hinzugenommen haben, die aus u=0 = konst. 
Re RR bzw. frei gelagerten Randes x = konst. hervorgeht. Die Formeln 
en so eine Linearkombination der Randwerte u, bzw. u,, mit u,,, die Null gesetzt 
Een nt ee exakt den gewünschten Randbedingungen gleichwertig, da die Bedingung 
he = n erungsfunktion nur angenähert erfüllt sein wird. Doch befriedigen die hier- 
aan Er A = a durchaus. — Den gleichen Zweck erfüllen die symmetrischen 
Heuichen en Fu die sich auf den Eckpunkt selbst beziehen. Sie verwenden außer der 
en er Iinsung ,=0 bzw. u, = 0 die triviale zusätzliche Bedingung u, = 0 
EEE Ken Re Kunstgriff zur Behebung der angeführten Schwierigkeit besteht 
en : a: Re u, bzw. u,, an mehreren Randpunkten gleichzeitig verwendet, 
Pe a arten tellen. Die so erhaltenen Mehrstellen-Randformeln RM 1 und RM 2 

gens keine AAu-Werte mehr enthalten, dürfen nur in Verbindung mit anderen Rand- 
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formeln verwandt werden, damit die Randwerte eindeutig festgelegt sind. — Lediglich für die 
Bedingung 7, = 0 in der Nähe der freien Plattenecke ist man auf eine Formel der Ausdehnung 
5 h der Genauigkeitsordnung A® oder aber auf eine der Ausdehnung 4 h der Ordnung h? angewiesen 
(Formeln E3 und E3°). — Die Sterne der Eckformeln sowie der Mehrstellen-Randformeln, denen 
die ecknahen Punkte 3 und 4 zugänglich sind, zeigen die Bilder 5 und 6. 


Bild 5. Stern der Eckformeln E, ES, EES Bild 6. Stern der Mehrstellen-Randformeln RM 


Auf diese Weise entstehen folgende Formeln, deren Verwendung in den verschiedenen Eck- 
Lagerungsfällen anschließend erläutert wird. 


| ES1 | —. a 
SERIE 8 0 
= oO 10 —ı24 0 8 —2|u—-|0 12 0|AAuht = 60 (u? + u) h+ (R), 
x 
— 149 [448] —ı24 ie. 1 0 
u, 1 10 
397 11 
Bea AA 12 
ss 2 w—|12 102 12] 4Au rt = — 180 (uß} + u) RR+(R?), 


| 3296 — 1586 


—1041 377 


=) 
—148 248 —148 
od 85388. "360... —148 


160 —388 248 


} 
1} 
L=z= — 
R,, [280 [156] 85 —148 


—20 
—64 152° —64 

+4: 325 —664 360 —64 ül 

—564 1120 —664 152 —20 

290 —564 325 —64 z 


44 
322 —l4 4 | Aduf-+ 
32 


7 
1 —10 7 |14uM= 72 (PA+T@)R+(M), 
1 


BESI | 
= 
4 
x [77 
| EES2 
ar ale nu sy, 
' (6) —311 —192 —156 AAu hk = — 60 (u + u) h? 
| +), 


148 [592] —_192 


—176 148 —3ll 
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—228 141 


159 
—712 984 —600 40 


84 
l 

1653 —2544 1782 —600 141|, __|12 —120 84|AAuh = 360 u@) h2 + (M?). 
—2412 3912 —2544 984 —228 s 


1384|—2412 1653 — 712 159 


152 —304 152 1 
238 212 —72 —4 —lEl, —|-5 68 —85| 440 M = 300 u@ h + 36009 %® + (R}). 
o 
Ur —1375 [1264| —880 288 —17 1 


428 —168 118 —4 —14 E> 


—83 


—664 5 


20 36 —106 —I1|,_2|-1 292 —151| AAu = — 180 u 124-900 upn2 


1996 752 —68 5 +(h}), 


—540 451 —106 —11 


[118] —308 324 —152 18 
—21 77 


E3 5. —ı110 55 
—220 440 —220 
' 320 —836 750 —320 86 ut: u— 
a [306] 1238 —1808 1268 —438 46 > 
—_ 45 —286 475 - 390 86 
0 0 
— [49 —102 53 |4Aum—u.|—2 6 — |AAu nt = 247, RB -4 (RP). 
0 0 


—104 114 —56 11 


—70 —238 112.910» 


35 —104 114 —56 1 


Ugz h? + (h?). 


HDD - 
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Mit Hilfe dieser Formeln lassen sich die verschiedenen Ecklagerungsfälle in folgender Weise 
behandeln: 


Eekrand | Randbedingung | z . Formel | Angewandt auf 
Ur a Mr = 0, in =: 
| y=0 ind 1 Punkt 1 
= 
07 
Uxx N) ee) 
ES2 Punkt 1 
uy=0 ind 
%yy 
U 28 
2 %=0 in? El oder RM1 Punkt 1 
] 07 UY=0 in? El oder RMI Punkt 1 
zei--09;-- %, = 0 in 2 EES1 Punkt 4 
x 
DE 02 Ur=0 in3 E2 oder RM2 Punkt 1 
Be er ggg EES2 Punkt 4 
Ux 4 
Ur =Uy=0 in2 EES2 Punkt 1 
I 2 
Fa 82 Ury = 0 in 2 EES3 Punkt 2 
a | en a a A 
DE ER 003 E3 oder E3° Punkt 3 
u 149 [9 y=-0"m4 ES3 Punkt 4 
Bei voller Symmetrie: ES3 Punkt 3 = 4 


Yxy=0. Uxx=Ugy=0 


—-4=0,9,=0 


a a 


Bild 7. Behandlung der Platte in Ecknähe: Innere Punkte O nach Haupt- 
formel, Randpunkte | |nach Randformeln, Eckpunkte [| 
nach Eckformeln 


Die Behandlung einer rechteckig begrenzten Platte in Ecknähe insgesamt veranschaulicht 
Bild 7. Auf alle inneren Punkte wird Hauptformel A angewandt, auf Rand- und randnahe Punkte 
die Randformeln und auf Eck- und ecknahe Punkte die Eckformeln. Den randnahen Punkten 
sind dabei je nach Lagerungsart die Bedingungen u, = 0, u, —=0 oder m, = 0 zuzuordnen, 
den Randpunkten für freien Rand die Bedingung q, = 0 mit Ausnahme der freien Ecke, die 
durch u,, = 0 erfaßt wird. 
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6. Einspringende Ecke 
Verhältnismäßig einfach ist die einspringende Ecke behandelbar. Hier lassen sich die Rand- 
formeln bis an die der Ecke benachbarten Randpunkte heranführen, diesenoch mit erfassend. Ledig- 
lich für die drei Punkte 1 bis 3in Bild 8 sind besondere Formeln aufzustellen, die aber keine Rand- 
bedingungen, sondern die Differentialgleichung anzunähern haben. Für den Eckpunkt1 gibt es eine 
symmetrische, für die Nach- w 
barpunkte 2 und 3 eine sym- 


metrische und antimetrische 2 RUTESTERFER 
Formel, deren Sterne Bild 9 ee © 
zeigt. Die drei Formeln lauten: re 
[o) oe 0 
03 e ee ©® 
: eo ae 
(@Xs 
En oe ® 
ö ® 


Bild 8. Einspringende Ecke 


— 
8 —|1mumM=0+(M), 
a 


74 


= 
—ı 4.-7la2Reneıı 
an 


1 
—1 0.0 |AAuM=0+(R). 

0 
oo ao. .oa.®®©®© Anwendung der Formeln: 
® © ® © © ® ® Formel | Angewandt auf 
I ae 1S1 Punkt 1 

IS2 oder IA Punkt 2 
® ® le] ® IA oder IS 2 | Punkt 3 
® © [e_® Die gesamte Behandlung in der Nähe einer einsprin- 

enden Ecke zeigt Bild 10. 
© © Be ö 
Bild 10. Behandl tt k i 

® © ee Tuc e 


nach Randformeln RB 
[1 nach Inneneckformeln IS und IA 


7. Berechnung der zweiten Ableitungen u,,, U,, 


ee u;, der Biegefläche durch Auflösen des linearen Gleichungssystems 
ermittelt worde j i i 1 
en sind, hat man aus ihnen die zweiten Ableitungen u,,, u,, zu berechnen, aus 


m 1, nach GI. (2 die Ableitungen lassen 
Mehrstellenp so nach einem Ansatz der Form (11) Näherungsformeln 

Er enauigkeit gewinnen. Bei gleichem Aufwand wie bisher wird eine Genauigkeitsordnung Re 

Be zu RP, Uyy h?, also die Ordnung h? für die Größen u,,, u,, selbst erreicht. — Für die Innen- 

_ punkte ergibt sich zunächst die Hauptformel en pn = 


TREE ME ee i 
4 aM 0 044 2” af e 


7 —196 [378] —196 alu—]2 — 2 pam. + 
4 —M4 40 7 —24 4 i 
b A 2 —4 = : 

die auch u,,, liefert, indem man die x- mit der y-Richtung vertauscht. Für die Randpunkte liegen die 
Formeln bereits vor als R2 und R3, die lediglich nach den Größen u,, h? und u,, h? aufzulösen sind. 
Für die randnahen Innenpunkte, also die vom Rand um h entfernten lauten die Formeln: 


ern Bern ke ge S Pe Ace 
‚ergeben. Auch für 


Fe E eyme 


7 


.ä 


Br ER 1 
93.4, 1.110: Z46 SuR'6 9 
[BR1] — 180 v,,.n2 = —129 [150] 9 —M 3lu+|o 11 —3 [44m + (R), 
a ee a? 


—ll 


2 a 
a Da8 4ü —12 
6 [474] —08 48 —s|u+|1ı 0 7]4A4uMt + (MM). 
= = 336 4202-12 
2 et 


[BR2] — 180 4,7 = 


= Eine größere Anzahl von Ableitungsformeln werden für die Eckpunkte und ecknahen 
; Punkte 1 bis 4 in Bild 4 gebraucht. Hier wird der Genauigkeitsgrad h® dadurch erreicht, daß man 
F. in die Formel einen der bekannten Randwerte einbezieht. Auf diese Weise entstehen folgende 
E neuen Eckableitungsformeln BE, in denen der Punkt, für den die Ableitung berechnet wird, 
1 durch Umrahmen hervorgehoben ist wie schon in den Formeln BR. Zu jeder Formel ist überdies 
2 durch eine kleine Skizze angedeutet, auf welche Randbedingung die Formel bezogen ist. In der 
Formel ist die verwendete Randableitung, die laut Randbedingung Null wird, mit aufgeführt, um 
die Formeln auch noch für eine allgemeinere Aufgabenstellung verwendbar zu machen, worauf 
a wir in Nr. 8 zurückkommen. 


s BEI 
4 ) 
£ Yx 7120 Uns h? = ee eh 
| -0- 
% 0 — 330 u@h + (h?), 
/ [882] 
| 2 
dw — 720 uyR? = 56 0|44uM + 
ax | 2 + 300 u® A + (h?). 


—7 14 —7| AAuht+- 
+ 180 u® h® + (A?), 


56 —140 126 —56 14 u+ 


—2231 14728] —2806 


21 —70 91 —56 


Y%x  — 1800 u? = 


E2I F — 
12 4 =02 5 
Z Eee 16 —|„ı|ı2 116 12|414u #— 
= 971 [1432] Ida 297223 5 — 180 u®. h? + (A”), 
en 100) 321 16.4 
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BE5 5 
U — 360 gr h?= w+| —24 48 —24 | 4Auht + 
ex 0 + 1500 u® Ah -+ (h”), 


24 —48 24 
—96 192 —96 
192 -480 432 —192 48 


1883] —2304 408 64 —öl 


72 °-—240 312 —192 


24 —48 24 


BE6 
96 : 192 —96 0 
— 360 Ugı h? = 192. —480 432 —19 48 |u+| —24 48 —24 | AAuMt + 
—19%2  +576 —672 384 —96 0 + 1500 u® h + (RP), 
Uxx [2147] —3120 1392 
Ux 
BE7 
0 
— 720 ig = 2472 AAu Mm — 
Ba x 300 «Ph + (7), 
-- >. 
Ux 
BES 
0 
— 1800 ug, h? = —714 — | Adult 
r 0 + 1800 u h? + (A?) , 
Uxx a 
[BEI —104 
> —208 1016 —208 175 
Le _ 3 31620 A 584 146 |, +2 —46 —2|AAu A + 
—3516 [7048 |—2056 —1168 292 25 + 180 u h? + (A). 
De An —544 1440 —1709 824 —161 


Die Anwendung dieser Formeln ist für die verschiedenen Eckverhältnisse im folgenden zu- 
sammengestellt, wobei die Punkte wieder nach Bild 4 numeriert sind. Ein * hinter einer Formel- 
nummer besagt, daß die Formel unter Vertauschen der x- und y-Richtung anzuwenden ist. 


Formel zur Berechnung von 


Rand Punkt Nr. 
Mn Uyy 
1 BEI BEI* 
2 01 0!) 
3 BE5 0 
nn 4 0 BE5* 
1 BE3 BE3* 
2 0 0 
3 0 0 
er 4 0 0 
il BE9 BE9* 
| 2 0 0 
N 3 Ugz = —MUyy | BE8* 
Lo 4 BES | Uyy = — MUgr 
1 BEI BE2 
2 BE6 0 
3 BE5 0 
4 0 0 
1 BEI | BE2 
2 BE6 1) 
3 BE5 0 
enmene 4 BE7 Uyy = — Ulgg 
1 BE3 BE4 
2 0 0 
3 0 0 Ari 
(et 4 BES Uy— — Un Hy = 0. 
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Es sei hier noch auf eine andere Möglichkeit zur Behandlung der Plattenaufgabe hingewiesen, 
nämlich auf die unmittelbare Berechnung der Ableitungen z,,, u,, ohne Verwendung der Funk- 
tionswerte u. Dazu macht man für die Aleitungen selbst je einen Mehrstellenansatz der Form 


3 Gün,t+ 2% A,, Adu,, = 0 + höhere Glieder 


I br, Uy,,+ % B,, AAu,, = 0 + höhere Glieder a 


£ Dazu treten wieder Randformeln ähnlicher Bauart, in denen die Randwerte u,, U, My My 
9,,q, eingearbeitet sind. Die Randordinaten u selbst treten nicht auf. Man erhält so zwei 
Systeme linearer Gleichungen für die interessierenden Werte u,,, , Uy,,,. Da zur Erhöhung der 
Anschaulichkeit indessen meist auch die Form der Plattenfläche u(«, y) selbst interessieren wird, 
haben wir diesen Lösungsweg hier nicht weiter verfolgt. 


8. Platte unter Einzellasten 


Auch die unter einer oder mehreren Einzelkräften belastete Platte ist mit Hilfe des Differen- 
zenverfahrens behandelbar, indem man die von den Einzellasten herrührenden Singularitäten ab- 
spaltet. Es genügt hier die Betrachtung einer einzigen Last P. Man zerlegt dann die Durch- 
senkung u in einen von der Einzellast herrührenden Anteil u,, der die Singularität erfaßt, aber die 
Randbedingungen nicht befriedigt, und den restlichen Anteil u, einer unbelasteten Platte unter 
Randbedingungen, die zusammen mit den Randwerten von u, die Randwerte der wirklichen 
Platte ergeben. Es wird also gesetzt 


DENT UL 1 2 en ae ee N) 
mit 
12 22 1 RE nr EN 
Als u, empfiehlt sich (vgl. etwa [1], S. 254) 
12 T z 
2 Eh ae er 
ee ee A EN EARREIN 


mit beliebig gewähltem r,. Nach Berechnung der interessierenden Randwerte von u, erfolgt die 
Berechnung von u, nach dem Differenzenverfahren unter p = 0 mit den negativen Randwerten 
von ı,. Mit Rücksicht auf diese Sonderform der Plattenaufgabe sind in allen Randformeln die 
Randwerte, die sonst Null zu setzen sind, mit aufgeführt. Auch in der verwandten Scheiben- 
aufgabe AAu = 0 treten von Null verschiedene Randwerte auf. 


9. Zwei Zahlenbeispiele 


Zur praktischen Erprobung des Verfahrens wurde eine Reihe von Beispielen durchgerechnet, 
für die exakte Lösungen in Form von Reihenentwicklungen bekannt sind. Davon seien hier zwei 
Beispiele wiedergegeben, bei denen verwickeltere Randbedingungen vorkommen, nämlich die 
dreiseitig gelagerte [7] und die über Eck eingespannte [8]*) quadratische Platte unter konstanter 
Flächenlast p = konst. Einfachheitshalber ist die Rechnung mit verschwindender Querdehnung 
u = 0 durchgeführt. In beiden Fällen wurde die Maschenweite h = a/4 bei einer Plattenseiten- 
länge a gewählt. Für Durchsenkung u und Biegemomente m,, m,, die bei der quadratischen 
Platte von der Form 


= 2 
„=m,pa 


4 

u=u,, me—m.D 0, m 

sind, werden die dimensionslosen Faktoren u, m,, m, berechnet. Nur für das erste Beispiel sind 

den Näherungswerten die exakten Werte der Reihenentwicklungen vollständig gegenübergestellt 

worden. Im zweiten Beispiel wurde mit Rücksicht auf die sehr verwickelte Reihenauswertung 

lediglich der Wert der maximalen Durchsenkung u der freien Plattenecke mit dem aus der Arbeit 

[8] entnommenen Reihenwert verglichen, zumal diese Arbeit die Werte der Biegemomente nicht 
explizit angibt. 

Um die Anwendung der angegebenen Formeln zu erläutern, sind im folgenden die linearen 

Gleichungssysteme aufgeführt. Die Formeln sind dabei teilweise mit passenden Faktoren multi- 

pliziert, damit die Diagonalelemente des Gleichungssystems von angenähert gleicher Größen- 


4) Auf das zweite Beispiel der über Eck eingespannten Platte hat mich freundlicherweise Herr Professor 
Karas, Darmstadt, hingewiesen. Gerade dieses Beispiel erwies sich für die Durchbildung des Verfahrens von 
außerordentlichem Nutzen, indem es die richtige und sachgemäße Art der Rand- und Eckbehandlung deutlich 
machte. 
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ür di i iminati i ei Seiten ist überall 
ordnung werden, was für die numerische Elimination ratsam ist. Bei den rechten Sei jeräl 
der Faktor (h/a)* = 1/256 weggelassen worden und ein Faktor 10 hinzugefügt. Die Werte u 
gehen also aus der unmittelbaren Lösung der Gleichung durch Multiplikation mit 1/2560 hervor. 


1. Beispiel: Dreiseitig gelagerte Platte 
Die Lage und Numerierung der benutzten Gitterpunkte zeigt Bild 11. Unter Ausnutzung 


der Symmetrie erhält man 8 unbekannte Gitterpunkte. Die Koeffizienten und rechten Seiten 


des Gleichungssystems zur Berechnung von u,, zeigt folgende Tafel 1, in der auch die jeweils be- 
nutzten Formeln mit angegeben sind. 


8 7 
ee a = 


Bild 11. Gitterpunkte der 
dreiseitig gelagerten Platte 


Tabelle 1. Koeffizienten des Gleichungssystems für die dreiseitig gelagerte Platte 


Formel |Nr. | il N 2 3 4 5 6 7 8 r 
R2 1 34,8 — 24,8 — 27,72 16,48 8,16 —2,88 —0,87 0 7,8 
ES2 2 — 15,86 36,24 4,32 — 16,30 —0,44 3,28 0,11 —0,37 15,0 
A 3 — 10 —4 36 —20 — 10 —4t 1 2 30 
R2 4 8,24 —13,84 — 27,12 42,96 8,24 — 13,84 —0,59 1,24 1.8 
R2 5 8,16 —2,88 | 27.72 16,48 34,80 | —24,80 —9,69 4,96 7,8 
RM2 | 6 en 11-4. P16,0. 04, 208 32,0 11.4, 14 16:00 a 
R3 M 16,0 —b6,4 —47,4 32,0 38,4 —27,2 5,6 — 12,8 18 
EES2 | 8 15,2 —10,4 — 15,6 —4N —19,2 48,8 —31,1 38,0 0) 


Auflösung des Gleichungssystems ergibt die u-Werte der Zahlentafel 2, denen die exakten 
Werte der Reihenentwicklung gegenübergestellt sind. Die prozentualen Fehler sind auf den 
maximalen Wert u, bezogen. 


Tabelle 2. Durchsenkungen % der dreiseitig gelagerten Platte. Vergleich mit den 
exakten Werten 


Nri ud 08 3 anne Bng 8 
z..10° 4,5262 3.2513 7,7439 5,5427 9.6421 6,8860 | 10.7223 7,6465 
Beakı 4.5412 | 32659 7.1549 5.5515 9.6457 6.8905 | 10.7260 7.6546 
Fehler | —0150 | —o1s | —o1o | —ooss | —ooss | —ooas | 0037 | —oosı 
%, —0,14% | —0,14% | —0,10% | —0,08% | 0,08% | 0,04% | 0,04%, | 0.08%, 


Aus den u-Werten errechnen sich die Biegemomente m,, m, der Zahlentafel 3, denen wieder 
die exakten Werte gegenübergestellt sind. Die prozentualen Fehler beziehen sich auf den maxi- 
malen Wert m, des Punktes 7. Die Biegemomente an den drei aufliegenden Seiten sind Null. 


In der Tafel sind die Formeln zur Berechnung der Werte m, = — u,, m, = — U,, mit auf- 
geführt. 
Tabelle 3. Biegemomente m,, m, der dreiseitig gelagerten Platte 

Nr. ji 2 3 + 5 6 7 8 
Formel BR2 BE3 B BRI: BR2 BE3 R2a BES 
My: 10° 42,1814 33,6808 73,0928 56,7547 91,7237 70,0633 102,4837 77,4506 
Exakt 42,1329 33,8442 73,1535 56,8880 91,7056 70,1808 102,3536 77,732] 
Fehler ‚0485 —,1634 —,0607 —,1333 ‚0181 —,1175 ‚1301 —,2815 
9%, 0,05% | —0,16% | —0,06%, | —0,13% 0,02% | —0,12%, 0,13%, | —0,28% 
Formel BR1 BE4 B BR2 BR1 BE4 R2 — 
m, + 10° 22,8436 16,8456 21,6232 15,4924 13,4564 9,5822 0,0515 0 
Exakt 23,1979 17,1016 21,6689 15,4631 13,3419 9,4546 0 0 
Fehler —,3543 —,2560 —,0457 ‚0293 ‚1145 ‚1276 ‚0515 — 
9%, —0,35% | —0,25% | —0,05% 0,03%, 0,11% 0,13%, 0,05% en 
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Die Ergebnisse sind offenbar über Erwarten gut, zumal wenn man den durchaus tragbaren 
Arbeitsaufwand in Betracht zieht, zu dem die ungleich viel mühsamere Reihenauswertung in 
gar keinem Verhältnis steht. 


2. Beispiel: Über Eck eingespannt — freie Platte 


Lage und Numerierung der Gitterpunkte zeigt Bild 12. Die u-Werte liegen symmetrisch 
zur Hauptdiagonalen. Bei den Biegemomenten sind die Werte auf und unterhalb der Diagonalen 
angegeben. Aus ihnen gehen die oberhalb der Diagonalen durch Vertauschen von m, und m, her- 
vor. Die Randpunkte sind hier mit Rücksicht auf die am eingespannten Rand auftretenden Biege- 
momente besonders bezeichnet. 


Tabelle 4 
Koeffizienten des Gleichungssystems der über Eck eingespannt-freien Platten 
Formel | x | 1 2 3 4 5 6 7 8 9 10 N r 
ES1 1| 448 —24,8 0 3,2 1,6 0 02 | —04 0 0 12 
R1 2| — 80 35,2 —14,4 —8,0 4,8 0 08 | —0,1 0 0 12 
A 31 —2 —20 36 —i —20 —2 2 2 2 0 30 
E1 4| —15,2 18,2 —3,6 63,2 —44,2 14,4 —8,4 520 00 15 
R2 5| — 1,44 16,40 | —27,72 | —13,84 43,04 | —12,40 2,48 | —10,28 8,72 0,13 7,8 
EES2 6 0 30,4 —15,6 508 —38,4 59,2 46,4 | —62,2 29,6 | —17,6 0 
EES1 a 24 —36 48 24 —24 84 —78 Dar 10:5 0 
R3 8| —32 32,0 —47,4 —16,8 54,4 —13,6 —6,4 43 | —4,8 1,5 18 
ES3 9 0 —29,6 36.0 49,6 —77,6 16,0 | —29,6 17,0 31,2 |—19,0 | 138 
EES3 |10 0,4 —12,0 17,82 19,68 | —50,88 89,11 |—14,24 | 33,06 |—48,24 | 13,84 m2 


Bild 12. Gitterpunkte der über Eck Bild 13. Biegemomente mx der über Eck eingespannt- 
eingespannt-freien Platte freien Platte 


Tabelle 5. Durchsenkungen 2% der über Eck eingespannt-freien Platte 


"Nr. |. | 2 


u. 4 
2-10 | 09350 | 2,241 17,5230|4,2350| 12,3266 | 20.2334 | 5,6660 | 16,6616 | 27,2652 | 35,9315 
I | ! | \ 

Der exakte Wert der maximalen Durchsenkung an der freien Plattenecke ist in [8] mit 
10°. u = 36,19 angegeben, der Fehler beträgt also —0,26 = —0,72%,. Auch dies kann in Anbetracht 
der Schwierigkeit der Aufgabe als gutes Ergebnis bezeichnet werden. Für die Biegemomente sind 
keine Vergleichswerte bekannt. Ihre Näherungswerte sind in Zahlentafel 6 zusammengestellt, 


und die m,-Werte in Bild 13 in perspektiver Darstellung über der Platte aufgetragen. 
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Tabelle 6. Biegemomente m,, m, der über Eck eingespannt-freien Platte 
Ka EEE MT 


Nr. 1 2 | 3 Fk 6 7 
IR. SE ER SEN RR BER EI a 2 tn, 
Formel. | Bl BR B BEI BR2 .BE9 BE7 
710° | 10,354 | —32,407 | 4,228 | —56,403 | 7,358 | 17,432 | —79,350 
ormel BE2 BR2 B BE2 BRI BEI ıh 
10 | 10354 | 2288 | 4228 | 3077 | 8692 | 17432 | 0 
BEN ER RENT N U SEE 

8 9 10 | 00 01 02 04 HE: 
ee ee BSP TE FERN VEEEEEREEREEREEEEEEEN VEREEREEEEEREEEREEEEEE 

R2a BES “= = BES R2 BES BE6 
_ 9,615 | 35,913 ) 0  \-AaL216 |-117,038 |-ı84,169 | —241,932 

) 0 0 | 0 0 0 ) 0 


10. Behandlung allgemeinerer Aufgaben nach dem Mehrstellenverfahren 


Das unserem Vorgehen zugrunde liegende Prinzip ist nicht auf die vorliegende Plattenauf- 
gabe beschränkt, sondern auf allgemeinere Rand- und Eigenwertaufgaben anwendbar. Für die 
allgemeine Randwertaufgabe der Differentialgleichung 

La] fe, u) Pre >= ln 
mit Randbedingungen der Form 
M][u] = 0 N: 


worin L[u] und M[u] lineare partielle Differentialausdrücke in u(x, y) bedeuten, lauten die den 
Ansätzen (10) und (11) entsprechenden Gleichungen. des Mehrstellenverfahrens 


D a,U,, + & A,, L[u]l,, = 0 + höhere Taylor-Glieder 
Mlu]; « = & b,,4,, + & B,, L[u],, + höhere Taylor-Glieder . 


Reh): 
1 


In dieser Weise läßt sich z.B. auch die schiefwinklige Platte nach Einführen schiefwinkliger 
Koordinaten angreifen, wobei die Gleichungen den Sinus des Achsenwinkels als Parameter ent- 
halten. —- Der vorliegende Beitrag möge in diesem Sinne auch als ein Beispiel für die Anwendung 
des Mehrstellenverfahrens auf partielle Rand- und Eigenwertprobleme überhaupt angesehen 
werden. 
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 Asymptotisches Verhalten der Zylinderfunktionen Be, 
In Abhängigkeit vom komplexen Index | 3 


“ 


Von Peter Beckmann und Walter Franz, Inst. in Münster 


Es wird für festes komplexes Argument untersucht, welche Sattelpunkte in den verschiedenen Gebieten der | Br 
komplexen Indexebene das asymptotische Verhalten der Zylinderfunktionen bestimmen. Dadurch ergeben sich 
in der Indexebene Verzweigungsschnitte für die asymptotischen Darstellungen. Auf diesen liegen die Nullstellen. u 


For constant value of the complex argument it is investigated which of the saddle-points in the different parts Y- 
of the complex index-plame determine the asymptotic behaviour of Bessel functions. It is found that the asymptotic I 
representations are related to the branch cuts in the index plane. It is on these curves that the zeroes are situated.  F 


Pour l’argument complexe et fixe on examine quels points culminants dans les regions differentes dans \le BE» 
plan del’index complexe determinent Vallure asymptotique des fonctions cylindriques. Par cela r&sultent des aLr. 
branches pour les representations asymptotiques dans le plan l’index. Sur cellesci se trouvent les points zero. - 


IIpıu hukcupoBaHHOM KOMNJIEKCHOM APTYMEHTE HCCHENYIOTCH CEIOBbIeE TOUKH B IIIOCKOCTH En 
KOMINIEKCHOTO HHNECKCA, ONPEeHENAIMIIME ACHMIITOTHYecKoE HOBeNeHHe PyHKIMÜ Beccenn. B » 
INIOCKOCTH HHNEKCA IIOJTyYeHbI IIMHMH PA3BeTBJIeHUA ACHMIITOTHYECKHUX BbIPa;KeHNÄ, HA KO- . 
TOPbIX PACHOJNO>KeHbI UX HYJIH. 
‘ 
b 


1. Einleitung : 
Durch Auswertung der Sommerfeldschen Integraldarstellung der Zylinderfunktionen Nr 


- Z,(e) = ii Io, 0, x) dx 


nach der Sattelpunktsmethode ergeben sich auf bequeme und übersichtliche Weise brauchbare 

asymptotische Formeln für große Werte des Arguments oder des Index. Die verschiedenen 
 Zylinderfunktionen unterscheiden sich nur durch den Verlauf des Integrationsweges im Un- 
‘endlichen der x-Ebene. 


Da die Lage der Sattelpunkte durch eine mehrwertige Funktion bestimmt wird, sind 
durch eine genaue Untersuchung des Integrationsweges die wirklich überschrittenen Sattelpunkte 

auszusondern. Für feste komplexe Werte des Index hat Emde [1] die Lage der Integrations- 
3 wege in der y-Ebene in Abhängigkeit vom Argument untersucht. In der Beugungstheorie ist 
& es jedoch notwendig, das asymptotische Verhalten der Zylinderfunktionen von festem kom- 
F plexen Argument als Funktion des Index zu kennen und insbesondere daraus die Lage der 
” Nullstellen zu ermitteln. Zunächst muß man den Verlauf der Integrationswege untersuchen 
Wir werden uns jedoch nicht auf eine Untersuchung der Verhältnisse in der y-Ebene beschränken, 
sondern unsere Ergebnisse direkt in die »-Ebene übertragen. 


2. Die Integraldarstellung für komplexe v» und o 
Die Sommerfeldsche Integraldarstellung der Zylinderfunktionen lautet: 


Zu) =. [etemrrirzay. a TR RL 


w 
Der Integrationsweg verbindet jeweils zwei Täler des Integranden, in denen ein negativer Real- 
teil des Exponenten für die Konvergenz des Integrals sorgt. Ihre Lage ergibt sich aus der Kon- 
vergenzbedingung. Diese lautet mity=u-+ iv: 


Re(-iesing+tivg)<o0 für NEE 2). 
Wir spalten sin y auf in die Summe von zwei Exponentialfunktionen und erhalten mit o = re‘? 
Re ee tor + üense 421”) <0 RE Re N 

r 


Für v— + ©o wird das Verhalten dieses Ausdruckes durch die jeweils größere der beiden Ex- 
ponentialfunktionen bestimmt. Für die Lage der Täler ergibt sich 


v—>+% tg>u>s+p EL ETAGE 
.D ——0o S—-p>u>—5— mad ee A 


2 
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Bild!) 1 zeigt die Lage der Konvergenzbereiche sowie den allgemeinen Verlauf der Integrations- 
wege für verschiedene Zylinderfunktionen. 
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Bild 1. Integrationswege für die Sommerfeldsche Darstellung 
der Zylinderfunktionen 


3. Die Anwendung der Sattelpunktsmethode 
Die Lage der Sattelpunkte ist bestimmt durch: 


0) ER : 
N TE ST TER 
eos a ET 
0 
Wir wollen mit x, den Sattelpunkt bezeichnen, für den gilt 
ru ae 
Entsprechend der Mehrwertigkeit des arc cos erhalten wir zwei Klassen von Sattelpunkten: 
VA, EEE ee ln 
und 
X. = —- + 2Iin. Bir. suoie er > IA er (6b). 


Bei der Sattelpunktsmethode wählt man als Integrationsweg die Fallinien des Sattelpunktes, 
d. h. die Linien auf denen vom Sattelpunkt aus gerechnet der Betrag des Integranden möglichst 
rasch abnimmt. Da die Höhenlinien des Integranden gegeben werden durch: 


Re(—iosiny-+ivy) = const 
lauten die Gleichungen der Fallinien über die Sattelpunkte: 
Im —iosiny+ivg)=Im(—iosiny,. +ivy,) = const. 


Dies ist aber auch die Gleichung der Steiglinien. Daher müssen wir durch eine Untersuchung 
des Integranden in der Umgebung der Sattelpunkte die Fall- und Steiglinien voneinander trennen. 
Im allgemeinen werden die Fallinien eines einzelnen Sattelpunktes nicht mit dem Integrations- 
weg in der Sommerfeldschen Integraldarstellung übereinstimmen. Daher ist gegebenenfalls 


der Weg so über mehrere Sattelpunkte zu führen, daß die Fallinien nacheinander durchlaufen 
die richtigen Täler verbinden. 


!) Diese und alle anderen Figuren wurden gezeichnet für 9726: 


— lo ie u 
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4. Die Lage der Sattelpunkte bei variablem » 
Bevor wir den Verlauf der Fallinien näher untersuchen, wollen wir die durch die Funktion 


Y.= Aalc 0087 N IE ee 5 6 ee EEE (7a), 


BELOICOS Tr A N ae AARAU) 


vermittelte Abbildung der v-Ebene auf die y-Ebene betrachten. Sie gibt bei gegebenem »v die 
Lage der Sattelpunkte an. 


Der Punkt v»= 0 wird in die Punkte y = > + mn abgebildet. Die Gerade durch die 
Punkte v = o und » = — o wird gegeben durch 
vol, a N N On A ES (8a). 
Für die Abbildung erhalten wir: 
Er a ft > (8b). 


Dem Abschnitt 0<t< + 1 entsprechen in der x-Ebene die Strecken 5 < <0,D =, 
mod 2x und die Strecken O<u< en) ‚ mod2r; dem Abschnitt —1<t=<0 ent- 


5 3a 


sprechen die Strecken 5= u <a, v=%0, mod 2r und zr<u<z 5° v= 0, mod2r. Für 


die restlichen Teile der Geraden wird y komplex. Mit 
cosy = cos (u +iv) = cos u Cof v— i sin u ©inv 
kommt 
Em 17 u=2mn = rkoliit. a BANSILcE (9a). 
ee wann FIR IHN „ea (9b). 


Dies sind Geraden parallel zur imaginären x-Achse. 
Für die auf der Geraden (8a) senkrechte Gerade durch den Punkt v = 0 gilt 


u ee ET St (10), 
Für die Abbildung erhalten wir: 
ee) u=,5 mod2r v= — MEinlil 
RN: 
0 u=—; mod2r v= Uröinlt 
a) u=5 mod?2r v= XEinlt 
er BEÜLLDE 
v<0O u=—; mod2r v= — ArEinit 


In der Umgebung des Punktes v = 0 setzen wirv=teir, t>0, und x =; Imn-+eer, 
Durch Entwicklung von cosy um den Punkt y = 5 + mr folgt 

P eiv.— peiv g ellako) em" N. le at. un 128), 

“=y—latg)—ıma RE ZERE RR  FL OD), 


In der Umgebung des Punktes » = D sind die Bilder der Geraden v = Lei? ebenfalls Geraden. 
Sie schließen mit der positiven u-Achse den Winkel & (12b) ein. 

Für 00 brauchen wir von den beiden Exponentialfunktionen des cosy nur die jeweils 
größere zu berücksichtigen. Aus 


fer = ne (ee Heut). WE NEN 
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erhalten wir so te 


vaio um 9 —y mod 27 | 
5 = 140% | 


Dv—>— oo um y—_p mod 27 
In Bild 2 ist die Abbildung der v-Ebene auf die y-Ebene dargestellt. 
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Bild 2. Abbildung der v-Ebene auf die -Ebene durch x = are cos an Bild 3. Verlauf der Kurven Im (—io) sing — (x — ma) cosy) = 0 


Die Zahlen geben den jeweils abgebildeten Quadranten der v-Ebene an 


5. Der Verlauf der Fallinien 


Wir betrachten zwei Sattelpunkte, y, und %,, die verschiedene Höhe haben. Von x, aus 
gesehen kann die Fallinie über diesen Sattelpunkt rechts oder links an x, vorbeiführen. Es 
kann aber auch der Fall eintreten, daß die Fallinie über den zweiten Sattelpunkt führt. Dann 
ist sie für x, eine Steiglinie. Die Bedingung dafür, daß eine Fallinie über zwei Sattelpunkte 
führt, lautet: 


Im —-iosny +ivg)=Im-iesng Firg). .: 2...) 

Zwischen x, und x, können die Beziehungen bestehen: 
2 Y=2mn—yı m= 0,1, 2,2... 2 Be 
%“%=2nn + n=l, 233,20: Do Re 


Im Fall (16b) liegen beide Sattelpunkte auf derselben Seite der u-Achse. Setzen wir (16b) in 
(15) ein, so erhalten wir: 

Im —iosng +ivyg)=Im(-iosng iv +2nivn) ..2..dı2. 
Diese Bedingung wird nur von rein imaginären » erfüllt. 


Im Fall (16a) liegen die Sattelpunkte auf verschiedenen ‚Seiten der u-Achse. . Aus (15) und 
(16a) folgt: 


Reife (sing —.[x mi] eösy))=0 Mi 

Auf der reellen x-Achse wird diese Bedingung erfüllt, wenn 
sin u — (U mm)eosu, 0. ne ee 
wu mE EIND. 


Die Lösungen dieser Gleichung bezeichnen wir mit u„. In ihrer Umgebung liefert (18) 
Re (er [& — U„) (U, — ma) sin u„, + (x — u„)? sin u, 


= oo cos u, (2 — [u — ma]) + --- ) = (52 20 
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Bezeichnen wir die Phase von (x — u„) mit o, so ist für u, #mz 


ET a N pe (21a), 
o= 5 mo ee (21b) 
und für .=mz 
Delete ae Er re (22a), 
o= an +0 — 2 ARD RATE MA (22b). 


Dies sind drei drei gegeneinander um z/3 geneigte Kurven durch den Punkty= mr 
Für große Werte des Imaginärteils von y ist y cosy das bestimmende Glied in (18). Für 
v — + 00 erhält man 


Bild 8 Bild 9 
Bild 4—9. Schematischer Verlauf der Fallinien für verschiedene Lagen der Sattelpunkte 


3 = RR  Eianger Mai Me 
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Für v—— oo kommt analog 
nee Der) =o re ee 
ee er 2 & „a2An): 


Nach (23b) und (24b) nähern sich die Kurven asymptotisch den Bildern der reellen v-Achse, 
In Bild 3 ist der Verlauf der Grenzkurven dargestellt. Die Schnittpunkte der Kurven mit ver- 
schiedenem m liegen auf den Bildern der imaginären v-Achse. a 

Um die Fallinien von den Steiglinien zu scheiden, entwickelt man den Exponenten von 
(1) in der Umgebung des Sattelpunktes. Bezeichnet man den Abstand vom Sattelpunkt mit 
ds e‘°, so erhält man 


— Vosiny,ı ir + sing aspe\&tr +2 ei A ee Ho 
Setzt man für die Phase von sin y zur Abkürzung ®, so erhält man für die Richtung der Fall- 
linien 
a 1) 9 Pe Re 
2 2 
Die Bilder 4—8 zeigen den Verlauf der Fallinien für verschiedene Lagen der Sattelpunkte auf 
den Grenzkurven. Ausgehend von Bild 8 erhält man den allgemeinen Verlauf der Fallinien. 


Wandert der Sattelpunkt S in Bild 8 auf der Grenzkurve nach links unten, so berührt 

die Fallinie II nicht mehr S, sondern läuft oberhalb vorbei und verbindet die Täler A und B 
(Bild 6). Verläßt jetzt S die Grenzkurve nach links, so läuft die Fallinie I links an S’”’ vorbei 
und verbindet TalC und B. Überschreitet S die Grenzkurve, welche das Tal C über s mit 0 
verbindet (mn—=0 im Bild 3), so läuft die Fallinie von C nach A. Wandert der Sattel dagegen 
nach rechts, so mündet I zunächst im Tale B’. Überschreitet der Sattelpunkt die nächste 
Grenzkurve (Bild 3), so verbindet die Fallinie das TalC mit dem rechts von B’ liegenden 
Tal, usw. Erreicht der Sattelpunkt das Bild der imaginären v-Achse, so läuft die Fallinie nach 
X= + ©o und führt über alle Sattelpunkte unterhalb der reellen x-Achse (Bild 9). Wandert 
der Sattelpunkt weiter nach rechts, so löst sich die Verbindung sämtlicher Täler auf in die 
Verbindung je zweier benachbarter 


B4 
GL; | Täler. Bei jedem Überschreiten einer 
| Grenzkurve führt jetzt die Steiglinie, 
L, die bei reinimaginärenvnachy =— oo 
2 ao > 5 | 3 K57 läuft, auf einen um x weiter rechts 
Q, liegenden Berg. Liegt der Sattelpunkt 
| innerhalb des von den Grenzkurven 


durch die Punkte 0, S’,rz und S in 
Bild 8 gebildeten Gebietes, so verbin- 
det die Fallinie die Täler C und B. 
(Hier ist (15) nicht hinreichende Be- 
dingung, wie z. B. aus Stetigkeitsgrün- 
den folgt.) 

Aus Symmetriebetrachtungen er- 
gibt sich das Verhalten in den übrigen 
Gebieten der y-Ebene. 

Zusammenfassend sind die Ver- 
hältnisse in Bild 10 dargestellt. Im 
Bereich 1 verbindet die Fallinie jeweils 
zwei benachbarte, auf verschiedenen 
Seiten der reellen x-Achse gelegene 
Täler. Im Bereich 2 verbindet die 
Fallinie zwei benachbarte, auf der- 
N selben Seite der reellen g-Achse ge- 
legene Täler. Liegt der Sattelpunkt im Bereich 3, so verbindet die Fallinie zwei Täler auf ver- 
schiedenen Seiten der u-Achse, deren Talsohlen in u-Richtung einen Abstand größer als x haben 
Die gestrichelten Kurven teilen die Fallinien, welche die vom Sattelpunkt aus gesehen nach rechts 
laufen, von solchen, die nach links laufen. 


Bild 10. Gebiete der x-Ebene, die sich durch ihren Fallinienverlauf unterscheiden 
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6. Übertrag in die v-Ebene 


Wir können jetzt angeben, über wieviele und welche Sattelpunkte die Integrationswege der 
verschiedenen Zylinderfunktionen führen. Für die praktische Auswertung unserer Ergebnisse 
übertragen wir die Grenzen der Gebiete 1, 2 und 3 von Bild 10 in die v-Ebene. 

In der Umgebung des Punktes » = o ist, w 0; inv—=0cosy, entwickeln wir nach Po- 
tenzen von y, und erhalten mit (22a) für die Grenzen 


N 
Tr Ir „a“ +4n+2nn) 


v.o+ a Fe ee ne 7 (27). 
Für | > r ist |Im (g,)| > 1 und damit 
0 
sing» + 1cosx;; X, = atc cos : s + ilog = DE LEE (28). 


Aus der Gleichung der Grenzkur- 
ven folgt: 


Im b (og - 1) 0 (29). 


Man setzt» = te'Y und erhält: 
Dee 
eh Lu 

Diese transzendente Gleichung er- 

gibt für jedes {> r/2 genau drei 

Werte von y: 


g>uy>Ga>m>5 


und 
TC u 23 
>23 3 > pp >—T. Bild 11. Gebiete der v-Ebene mit verschiedenem Verlauf der Fallinien 
Für die vom Punkte v» = — o ausgehenden Grenzen ist lediglich y durch y» + m zu ersetzen. 


Bild 11 zeigt die Abbildung der Gebiete 1,2 und 3 von Bild 10 auf die v-Ebene. Die Punkte v;; 
und »_; liegen zwischen + i|o| bei rein imaginären og und + 1,5 ilo| bei reellen 0. 


7. Die Nullstellen der Zylinderfunktionen 
Jeder der Sattelpunkte liefert einen Beitrag 


Ben 7 
J) —oSInX; tive iz 
V eh re re. in; (31). 


TEOSIN X; 


Sein Betrag wird im Wesentlichen bestimmt durch das Glied iv X, Im Exponenten. Liegt einer 
der Sattelpunkte in einem der in Bild 10 mit 3 bezeichneten Gebiete, so führt der Integrations- 
weg über die Sattelpunkte x, und x, sowie einen oder mehrere der gegenüber x, um ein ganzzahliges 
Vielfaches von 2 verschobenen Sattelpunkte x, + 2/r, mit + für Im) <0 und — für 
Im (v) < 0. Die Beiträge der Sattelpunktex, + 2 Irr können daher gegenüber dem Beitrag des 
Punktes z, vernachlässigt werden. Auch in den jetzt verbleibenden Gebieten, in denen der Weg 
über mehr als einen Sattelpunkt führt, liegt im allgemeinen der eine Sattelpunkt höher als der 
andere. Die Bedingung dafür, daß bei Sattelpunkte die gleiche Höhe haben, lautet 


Rel-iosing tivg) =Re(-ieosing + iv) A. Der‘ 
Zwischen x, und %, bestehen die Beziehungen 
Kae re ii oder u a 


Wir diskutieren lediglich den Fall yg = — xı; der Fall yg = — %ı + 2 läßt sich durch die Trans- 
formation y* = Yı FW xX% = Xa + v auf den ersten zurückführen. 


Zunächst betrachten wir die y-Ebene. Für kleine x, ergibt sich anälog zu (22a) 


neo Be tan Se. 
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Bezeichnen wir die Phase von x, mit ß, so erhalten wir: 


Re ee N in 
BP=z”% ap (34) 


Dies sind drei um 7/3 gegeneinander geneigte Kurven durch den Punkt y = 0. Für große Werte 
des Imaginärteils von „=u-+ to ist 


Py40o mean DStP 0. re nen er Ci. 


een u= (an —1)5— ne Le ee 
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.- == oo oes[ 
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Bild 12. Verlauf der durch die Punktey = mr gehenden Kurven 
Re (—ie (sing — (x — ma) cosx)) = 0 


Die Kurven nähern sich asymptotisch den Bildern der imaginären Achse. In Bild 12 ist der Ver- 
lauf der Kurven dargestellt. Die von den Punkten y=r und y = —r ausgehenden Kurven 
entsprechen den Fällen = — yı +2r mit =, —2r und y = % 


Wir übertragen die Kurven in die v-Ebene und erhalten analog zu (27) für die Umgebung 
von Mo: 


r |xs[® s (2 ++ =) 
2 Ä 


vs0+ ie) 
Für den Verlauf im Unendlichen der v-Ebene ergibt sich analog zu (30) mit » = t ei? 
ey  — urn ee re 
log — —1 


Diese transzendente Gleichung ergibt für jedes 2> r/2 genau drei Werte von vy:p <y, < 


37 7 
lo und > 


Für die Kurven durch v» = — po ist lediglich y durch y + x zu ersetzen. (Bild 13.) 


Der Integrationsweg der Hankelfunktion 1. Art führt in der Umgebung der in Bild 13 mit 
h, und h_ı bezeichneten Kurven über zwei Sattelpunkte etwa gleicher Höhe. Auf den Kurven 
selbst liegen die Nullstellen. Für die Hankelfunktion 2. Art gilt dasselbe für die mit Ah, und h, 
bezeichneten Kurven. 

In den Gebieten der v-Ebene, in denen der Integrationsweg jeder der beiden Hankelfunk- 
tionen über nur einen Sattelpunkt führt, läuft der Integrationsweg der Besselfunktion über zwei 
Sattelpunkte. Dies ist der Fall, wenn die Sattelpunkte y, und — y, beide in dem in Bild 10 mit 


> 
2 ’ 


| 
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1 bezeichneten Gebiet der y-Ebene liegen-(siehe auch Bild 11). Die Sattelpunkte haben die 
gleiche Höhe auf der Kurve (32), die den Punkt » = o unter dem Winkel m + 5 verläßt (jin 


Bild 13). Tritt diese Kurve in das Gebiet2 (Bild 11) ein, so führt der Integrationsweg der Bessel- 
funktion über drei Sattelpunkte: y,, — x, und 20 —x,. Der Sattelpunkt 2” — x, liegt rechts 


Bild 13. Nullstellenkurven der Zylinderfunktionen und Gebiete mit verschiedenen 
asymptotischen Darstellungen 


von der Nullstellenkurve von H\(o) tiefer als y, und — x, so daß x, und —x, das asymptotische 
Verhalten der Besselfunktion bestimmen. Links von der Nullstellenkurve von H}”(o) ist das Ver- 
halten bestimmt durch die Sattelpunkte x, und 27 — x, die auf der negativ reellen »-Achse 
Beiträge gleicher Höhe liefern. In der Umgebung der in Bild 13 und j bezeichneten Kurve müssen 
für die asymptotische Darstellung der Besselfunktion 2 Sattelpunkte berücksichtigt werden. Auf 
7 selbst liegen die Nullstellen. 


8. Asymptotische Formeln für die Zylinderfunktionen 


Die folgenden Formeln gelten für — > <arg(o)< > Die Formeln für größeres oder 
kleineres arg (o) ergeben sich mit Hilfe der Umlaufsrelationen für die Zylinderfunktionen [2]. 


a) Die Hankelfunktionen 
Die Punkte» = o und» = — e sind Verzweigungspunkte für die asymptotische Darstellung 
der Hankelfunktionen. Als Verzweigungsschnitte wählt man die Kurven, auf denen die Null- 
stellen von H®(o) bzw. H‘®(o) liegen (h, und h_ı bzw. h, und h_ 2). 
In der Umgebung von » = + o werden die Nullstellenkurven mittels des positiven 
reellen Parameters s dargestellt durch 
A + für H%’(e) 
tvwotel®se ° — für H®%g) 


Für den Verlauf im Unendlichen der v»-Ebene gilt asymptotisch mit +» =! evunde=re? 


wart 


bey. een: 9) 
log = —— 1 
mit 


0<y<, für H9@, 


5 <y<0 für H®(o) (s. Bild 13). 


Die Nullstellenkurven von H®(o), für 9 > 0, bzw. von He), fürg <0, schneiden die reelle 
»-Achse in den Punkten + v,. Bild 14 zeigt |vo|/\ol als Funktion von Q. 
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Mit v/o = cos y lauten die asymptotischen Formeln für H“®(o) und H}”(e) im Bereich I 
und II (Bild 13) 


2: DEREN: 
H®(p) m EST A "een <a el aller 
ee er (20) 
a1 
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Bild 14. Schnittpunkte der Nullstellenkurven der Hankelfunktionen mit der 
reellen v-Achse in Abhängigkeit von der Phase des Arguments 


mit 

VERNEREIETR! PETE EDIT ARE 
im Bereich III (Bild 13) 
2 > |/ Zi eiesingtivz, —n=< al = We ee 
zesiny 

H2 (0). HIN 0) = 20 Sa en a ee ee a 

mit 
Im(y) <0, —a.< Reyy< m. 7.2, Pelle 


(Die Differenz zwischen H$”(e) und H}?(o) ist die stärker als exponentiell verschwindende doppelte 
Besselfunktion.) 


Im Bereich IV (Bild 13) 


H&%Xo) ER IR —tesingtivx, —ız< et) ER re (42a), 
zesinx 
Hilo) > er He (0) Se 
mit 
Img) >0, —2r <Reg)<0 urn ee 


Auf den Nullstellenkurven ist jeweils die Summe der in den angrenzenden Gebieten geltenden 
Darstellungen zu nehmen. 


Mit den für |Im(y)|> 1 zulässigen Näherungen 
sing» ticosy, 1m tilg 


vereinfachen sich die Formeln. 
In den Bereichen I und IV ist 


ar Ei 
He) a 5 Re ie ne 


und in den Bereichen II und III 


Dr EION\ ES 
ji 2 Var [ek 
Ole: (43b), 


u u ee a see 


Z., angew. Math. Mech. s 5 r 3 5 
ER N ee eine 


Für Hf(e) gilt in Bereich I und III 


„2 ee: 
oo) ni) (ee) Pe 
BE a, re rede) 
und in Bereich II und IV 
2 eo” 
H® ae en ee 
(o) er (44b) 


b) Die Bessel/funktion 


Die asymptotische Darstellung der Besselfunktion hat ebenfalls einen Verzweigungspunkt 
bei v = o. Den Verzweigungsschnitt führt man längs der Nullstellenkurve der Besselfunktion 
(j in Bild 13). In der Umgebung von » = o lautet deren Parameterdarstellung 


eg- gs. „T.. lie FEEUNETZ REINE 


Links von der Nullstellenkurve von H{”(e) im 2.Quadranten der v-Ebene fällt die Nullstellenkurve 
der Besselfunktion asymptotisch zusammen mit der negativ reellen v-Achse. 
Das asymptotische Verhalten der Besselfunktion ergibt sich mit 


Io) = 5 (HE) + Ho) 


links von Ah, und h, (Bild 13) aus den asymptotischen Formeln für dieHankelfunktionen. Außer- 
halb der Nullstellenkurve der Besselfunktion braucht hier nur die jeweils größere der beiden 
Hankelfunktionen berücksichtigt zu werden. 


Rechts von A, und h, ist die Besselfunktion stärker als exponentiell klein und wird asympto- 
tisch gegeben durch 


Zi Sr : 
ES —tesiny-+ iv 
J,(0) + Vans esinx EIER a) 
mit 
Img) >0, — ra <Rey)<+n. 
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Stresses in certain thin elastic plates rotating about normal axes 
By A.M. Sen Gupta!) in Ranchi Bihar (India) 


Es wird eine direkte Methode entwickelt, um das Problem gewisser gleichmäßig rotierender dünner elastischer 
Platten mit krummliniger Begrenzung zu lösen. Zur Illustration der Methode wird das Problem einer dünnen 
rotierenden Platte von der Form eines regulären krummlinigen Polygons allgemein gelöst. Die entspre- 
chenden Probleme für einige dünne Platten mit Begrenzungen von einfacherem Typ werden ausführlich 
behandelt. 


In this paper a direct method for solving problems of certain thım elastic plates having curvilinear boundaries 
and rotating steadily about normal awes is developed. To illustrate the method the problem of a rotating thin 
plate having a regular curvilinear polygonal boundary, in general, is solved. Discussions of the corresponding 
problems of thin plates with some simpler types of boundaries are given in detaul. 


Dans cet article une methode directe pour resoudre des problömes de certaines plaques €lastiques et fines avec 
des: bords curvilignes tournamt continuellement sur des axes normaux est developpee. Afin d’expliquer la 
methode le probleme d’ume plaque fine et rotative avec um bord curviligne et polygonal est r&solu en gen£ral. 
On traite en detail des discussions des problemes correspondants concernant des plaques fines avec quelques 
types de bords similaires. 


B Hacrosmeii pa6oTe u3taraerca UpAMOoH MeTon pemieHun 3anau TeopuM TOHKUX INTACTUH 
C KPUBOJIHHeÄHBIM KOHTYPOM, BPalmaloIlmxcaH C TIOCTOAHHOÜ CKOPOCTBbIO BOKPYT Ocefi, HOPMAJI- 
BHBIX K CPenHHHOU MOBepxHocTu. DB KayectBe IIpuMepa pellaetcH B o6IMeM BHNe 3amaya 
BpaııeHnA TOHKOÜ IIIACTUHKH, KOHTYP KoTopoi mpencTaBiser co6ofi IpaBaJIbHbIa KRPUBOJI- 
NIHejÄHLIÜ MHOTOYTOJIHHK. ‚usa mpocrelimx BMUMOB KOHTypa COOTBETCTBylolan 3anaıa pacc- 
MaTpuBaetcn 6oJlee IIONPOÖHO. 


1) Professor in Mathematics, Birla Institute of Technology, P. O. Mesra, Ranchi (Bihar-India). 
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ee 
Nomenelature 


The following nomenclature is used in this paper. 
x, y = Rectangular coordinates. 


£&,n = orthogonal curvilinear coordinates. 
22, 99, &J = average normal and shear stresses. 
Er Ep &, — normal and shear strains. 


ar? up? 72Y 


= DUB &n average normal and shear stresses in curvilinear coordinates. 
z-iysz—r—-ig. 
E+in Befr. 
angular velocity. 

density of the material. 

— Poisson’s ratio. 

— Young’s Modulus. 

Modulus of rigidity. 

+. 

a plane harmonic function. 
Conjugate harmonic functions. 
stretch ratio. 

angle between the curve & = constant and the x-axis. 
87 ,.02G 8. 

real positive constants. 

a positive integer. 


(| 


o 


Bougbbr mann 
S 
= 
. 


>> 


1 er 


Ss. 


3 0° 0° 
— the Laplace’s operator Fr + De 


Lu} 


we 
® 


“real part of”. 

“Imaginary part of”. 
1+%(n+1%, Q=2/A(n+|]). 
1+%(a+lD), G=-Am+2. 


= 


II 


Introduetion 


A new direct method of solving problems of steady rotation of certain thin elastic plates 
about normal axes has been discussed in this paper. If the boundary of the plate be given by 
the curves & = constant or „= constant where + in = f(x + iy) , it is shown that the stresses 
can be expressed explicitly in terms of these coordinates and by suitable choice of a potential 
function, the stresses can be easily determined. For simpler problems of thin plates with holes, 
under no body forces, a similar method was employed by Sen?). But in problems of a plate 
rotating in its plane about a normal axis, the reversed mass accelerations- are to be treated as 
body forces, and the method requires further development. One advantage of the method is 
that once the general forms of the stresses have been deduced, the unknown potential functions 
involved can be written down without much difficulty and in many cases solutions are obtained 
in finite terms. To illustrate the method the problem of a thin rotating plate having a boundary 
in the form of a regular curvilinear polygon of n-sides (n > 2) is solved. 

For suitable values of n and the other constants involved in the equation, it is found that 
the curvilinear polygon takes the form of a two-cusped epicycloid, a dumbbell or a cog-wheel 
(cf. Rothman?). 

It is found that the method can be successfully employed in solving problems of rotating 
thin disks bounded by such simpler curves as ellipses etc. 

It may be mentioned in this connection that the tentative method of this paper differs 
from the function-theoretic method of Muskhelishvili*) who solved plane problems of elasticity by 
suitable conformal mapping of the region on a unit circle. It is believed that the problems con- 
sidered in this paper have not been previously solved by any investigator. 


Method of solution 


Let the plate rotate about a normal axis passing through the origin of coordinates with a 
constant angular velocity ®. Assuming the reversed effective forces as the body forces the pro- 


*) Muskhelishvili, N. I. Some basic problems of the Mathematical theory of Elasticity. Translated by 
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blem can be treated as a statical one. The plate is assumed to be in a state of generalised plane 
stress. Thus, the equations of equilibrium, satisfied by Xx, yy and xy are 


27 Ze 57 EI x 
La er e ER 
om BE ER ae 
de Toy ER 
The stress-strain relations are 
7, (2 — 0yy) 
WE (4.002). Mel a ee Bee (1.2) 
a 
ey = u" 
The compatibility relation is 
VER ER ER (1.3). 


a a Away 
If we put SH Er Ta A (1.4), 
the compatibility equation (1.3) can be expressed as 


a U a 2 7 ae a a 3 (1.5) 


Asa solution of (1.5), we can write 


2-9 — zu RO Er ee I (1.6). 
Using relations (1.4) and (1.6), we get, from the equations (1.1) 
ger ar 
am e “10 zei 0 
N (1.2); 
Et } RR _ 
or vila+ 3 +0eorr + = —N) 
Similarly from the equations (1.1), we get 
ep) ae 
Rn im +- (1+0)oo r + 3 n U Re (1.8) 
4 0y 
2, 
and , 4m ul ee (1.9). 


The equations of equilibrium (1.1) and the consistency equations (1.7), (1.8) and (1.9) will 
be satisfied if we write 


40% x 0%, 2 


1 
+ dl +oeorr- + ler) +y. . (110) 


Asa n 4 9% 
y 92 x 92 2 1 ng 
y+- Er +0)ew!r ee N ar 4 2 a. ze — yo (® — )—y (1.1) 
EIERN 1.12 
and RETTET Re I ek er 


in which wand ® are conjugate harmonic functions given by 


B+iY=faH+ig. 
Let us introduce the curvilinear coordinates defined by Z = F(£Ö) 
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so that we have®) 


HB) 
10% 1 a om re 
C& 0 = Toe Im’ ein 0-1 Im RE 
1 er RE Re E 
&j 9) — Ein 02, - In 
ee 
Coj + in u TE 
n Dog de 
and Yız en =). 


Also since the sum of the principal mean Ki is REEL of the choice of rectangular 
axes, we have K 
ze uU = EEE m Se ee (1.14) 
In terms of the stresses X%, 79, X the stress components ££, 7n and £n are given by the relations 
&E = Co? 0 + 2% Sin O Co) © + yy Sin? © 
m- R&n 9 —-297&5n0C5 O+MERO |-::::- (1.15). 
&n= (39 — &R) Sin O Cof 0 + 7 (Cof? 9 — Sin? O) 
Let, F, and G, be two conjugate harmonic functions defined by 


FR, +iG =—i@+ip(% +). 
_ w|[o% 0% 0Y 
nen: = 2 &] u: se 
: ay ray se rEE : 
ee. 
Using relations (1.13), we get, from (1.16) 
73 Fı —= / (Co? 0 — Ein 0) +2 9860) 9 Sin 
BR > (1.17) 


6, - 2 Ein 960) 0 — 5 (Col? O— Ein 9) 


Substituting the values of £z, yy, %y, from (1.10), (1.11) and (1.12) in q. 15) and using (1.17), 
we get 


— 02090 022902 = 
S Jr z T— "+41n N —-2( J? 0 w® 
ET a a re 
— 2(1— 0) o w: Re(2) Re 2) A ae (1.18) 
ie; 02290, , 92209 % 
Bee RE 2. 0 TE . 
sEm= AkER: DE 08 +427209,—2(1 + 0) J? o wX zz 
Ben 
+2(1— 0) o w Re(2?) Re (2) ENTE RELTEE ER (1.19) 
_ 02290, 02200, ® he 
I ee) iur. 2 
8 J?&n dm 8 DE dh +2 (1 — 0) o w® Re(z?) Im Mr: GER ZU 
The above results can, otherwise, be put in the form 
EEE — 5, x HERe) __%z8 
022 0] 027 0.22 
ee) en . 4 277 
a, 7; “1 2(1l +0) Jto wzz 
— 2 (1 — 0) Je w® Re(2) Re ( =) dr Ja Fa a (1.21) 


*) Coker and Filon, A treatise on Photo-Elastieity p. 153, 164 (1931). 


2. th. M' , B . . ß 


a 0 
8 Jtm Br (PA) — 2 2 PD +4 
2 
o a ieh 


oldE &E Om m 
+2.(1— 0) J?o @ Re(2) Re) ®F Kan. 2% 122) 


29 0 „ a dzzop 
ee. efeke, 


+2(1— 0) J2o w: Re(z?) Im 2) UPC. (1.23): 


|-2a + 0) Jto w® zz 


Im =— 


Stresses in a rotating regular eurvilinear polygonal plate of n — sides (n > 2) 
Let the transformation be given by 


2 -Ieexp (DL exp ln) ee er (2.1) 
where0<(n+1)}<1. 
From (2.1), we have 
2 = c[exp (of n +Aexp(n +1)ECof n+1) 7] 
y= c[exp (d)Einn+Aexpn+l)&ECinn+Dn] 
z=r=e[expR2&)+rRexp2(n+1)E+2Iexp(n+2Y)Eloinn]) . .. . (23). 
The curves & = 0 are in general, some regular curvilinear polygons of n-sides having n rounded 
corners or vertices (cf. ee) provided 


n+2 
e Ir 11 
(r)e=o is stationary when ©&in nn =0i.e. when 
N= ar : m=0,1,2,...n (numerically). 
n 


We have 
J? = @ [exp (28) +42 (n + 1? exp2(n+1)E+2/(n+1)ex pn +2)8Cojnn] 
Re(2) = @ [exp 29C02n+%2exp2(n +1)ECof?2(n +1) n 
+2Ixp(n+YVElinm+ 2] 
%) — @[exp (2L) +4: (n +1)?exp2(n+ )E+2An+Dxpn+2)/] 
[leo = @ IP + Q &ofnn] Se 


%2| _9e[P,+Q,Cojn] 
08 Io 


on 
We assume 


ker) — —2c@An©innn 
E=0 


Bexp@)+Dep(2n+1)ö 


BE Teaser (2.5) 
where A, Band D are real constants to be determined. Thus, 
P2,=e2[A{xp2@9+R%(n+1’exp2n+ 1) +Besp (2) 
+A(n +D{@A+B)exp(n+2)E+Dexp(n+2)&}Coinn 
+ Dexp2(n + 1) & Co] 2nn] 
[J2 2Q,lk-o= ce [AP+B+Am+1)2@A+BH+D)ECoinn + DEo| 2nı] 
(2.6). 


2:20] = aßp{1l+2n +1? A+2B 
5 &=0 +1 +Dfn+9@A+B)+@n+YD}ECH m 
+2(n+1)DEoj2nn] 


20] -— — enß(a+1)@A+B+D)Einnn+2DEin2n] 
&=0 


5) Stevenson, A.C. Phil. Mag. Ser. 7. Vol. 34. p. 112 (1943). 
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Substituting from (2.3), (2.4), (2.5) and (2.6) in the equations (1.21) and (1.23), we have, 
from the boundary conditions [&&]e-o = 0 and [Enle=o = 0 $ 


[I (F + iG)l=o = c! [co + a P Han +dDa+ta +nQ+mP+An+D)a}Coinn 
+{5+@,P+A(n+1)a}&o2nn+c6oj3nn 
+1{8,6innn+ S,6in2nn+S&in3nn}+Aln+l)gexpid—n)n 
+ {a +A(n+1)a,} exp(idn) + {a P+Aln +1) (a +a,)} expin+4n 
+{, P+Ian +), +a)}expi@n+4n 
+{, P+ian +), +a}expi@n+4Mn 
+{a, P+A(n+1)a,}expi(dn +4)n+Aln+1)a,expi(ön+4)n] (2.7), 

in which 

o=—2[@ PR +9)A+{2+%n +1)” +6n+4M)}B+2(n+1)%D] 

a=—8APQ—Q{n+6+R%(n+1)Bn+29}B+241{2 (m — 2) —R(n+1(n+2)}D 

= —(2A+B)@+{4n—1) +2%(n +1) —2n—%4}D 
g=(n—DQOD+R®(n+){l1—- om +2(n +) +o)}o wc 

S= (PR —-2nQB+2An@n+3+P)D+(—oAnPPP, +in+1)Qleo:e 

S,=2n{2 +2 (n+V)(n+23}D+Al—o)AnQ@P, +QPeoew®’e 

S=nQD+(Al—o)A%:nn+1)Qowd 

a={ß + (Pi +1Q) +2 +&n}oore 

Ge2Bß +0) AH Aew ec 

g=R{l— oO) +2(n +) +n)}ow@:« 

g—=(l—o)ow cd 

4=2( —0)Q,0 wc 

= (l—o)(®+6n +6)2o wc 

= A-NIQN Eure 

and „a=(l—o)(n+1” Mood. 

We assume 

3 


A, nn De Fa 
1 m=0 
EB (le) + An p.h)exp (ne DL. 2 Be 


in which A,, B,, are real constants to be determined. 
From (2.8) we get 
[I (F+i@)le-o = En +1) Ar + {A, +A(n + 1) Ayexp(inn) + {As+Aln +1) As} exp (in) 
+A,exp(iSnn)+Aln+1)B,expid—n)n 
+iB tn +) Byespin+Mn+{Bn+An +1) B}exp(i4n) 
H{B,+Aln +1) B,}expi@n+M)n+{B;+A(n+1)B,}expißn+4)n 
H{B,+Im+D)B}espidan+Mn+B,expiön+AM)n] (2.9). 
Equating the coefficients of cosines and sines of like multiples of n in (2.7) and (2.9) we 


obtain sufficient number of independent equations to determine the unknown constants. 
Thus, we get 


Fti6= 


0 w*c? 
en eu ie 
B ER 0 w?c 


SIR HDi Fe) fr —3—o(n+l) +2 (n+D)4n+3-+o0)] 


D-URony 3 ne 

A, =An{5—n+on—1)—R(n—1)An+5—o)}o wc 
4,=%Rn{7n+l0 +on—2) +R(n+1l)(Adn+5— o)}o wc 
A; =M®n(n+1){3n+5+o(n—1)}o wc 
B,=(l—o)oo?:c 

B,=/l—o)®n+5)o wc 

B, =# (1 —o) (3n? + 12n + 10) 0 w2 ci 

B, =#®(l—o)(n +1) (n?+8n-+10)o0 wc 
B,=Ml—o)(n+l®@n+5owcd 
B,=#(1l—o)(n+1)ow:c. 
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Values of Q,, F and G being known the stress components EE, nm and En are completely 


determined. 
On the boundary £—=0, £E=0. So that [mmle=o can be very easily determined from (1.14). 
We have 
rn 1 
[Imnle=o = [21-0 = Be 1+0oor ber 
E 0 w* 1 
"P+0QCHnn mu a ee a en mel Zu) 


—A:(n+1){4n+6n +3—o(2n +3)}— (1 —o) (n + 1)°4°] 
zer] DELEELUESOTTELE) FE (2.12). 


For 4 = 0, the curve & = O reduces to a circle. PuttingA = 0 in (2.12), we get 
u 2 | 
me dom ee re (2.13) 


R,(= c) being the radius of the circle in that case. 
This is a well-known result found in all text books on the subject. 


Partieular eases of interest 


Case (i). Whenn = 2andA< 1/3, the curvilinear polygon & = 0 takes the shape of a dumb- 
bell. (Fig1 a) 

Thus, puttingn =2 and A<1/3 in the results of the last section we obtain the components 
of stresses in a thin plate of the shape of a dumb-bell rotating steadily about a normal axis. 


Por nz =2, ad) = z (< Si we have 
= 1; 0 


11'345 else —. {13 &0| 27 + Coj 4 „| 


The peripheral stress [7M]-o in this case, for different valus of A (< = are shown graphi- 
cally in Fig. 1b, o being taken as 1/3. 3 
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Case (ii). When n is large, the curvilinear polygonal boundary & = 0, becomes approximately 
that of a cog-wheel (with shallow teeth). 


L.- i R 
Thus puttingn = 20 andA< 5" the results of section 2, we get approximate values of 


the stress components in a cog-wheel of twenty teeth rotating steadily about a central axis. 
In Fig. 2 the peripheral stress [nle=o in a cog-wheel of twenty teeth are shown graphically 


for different values of} |< ae co being taken as 1/3. It is found that the general character of 


the peripheral stress does not change due to the change in the value of}. For all values of ), the 
maximum and minimum peripheral stresses are occurring at the same places viz., at places where 
3 
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[r]e=o is stationary. But the ratio of the maximum and minimum peripheral stresses diminishes 
ası decreases and reduces approximately to unity as A approaches zero. 


Ye zo/gwic: 


Fig. 3 


A plate in the form of a two-eusped epieyeloid rotating about a normal axis 
Putting n = 2 in (2.2) we have 
ee 2 N 
y=cle©&nn+Ae’Gin3n] 
The boundary & = 0, will now be a two-cusped epicycloid if we take? = 1/3. (Fig. 3) 
2: we get 


= e|&ojn+ 4 &3n] 


Thus, when = Dand i—= 


; 
y=e&inn+ 4 &insn] 


These are the equations to a two-cusped epicycloid, the radius of the rolling circle being 


half that of the fixed circle, the points where n = + = 


Thus, puttingn =2 and} = 1/3 in the results of section 2, we obtain the components of 
stresses in a two-cusped epicycloid rotating steadily about a normal axis. 
In this case, we have 


er 1 
[Imnle=o = i8 1—-)8+&j2M)oote. 


‚ being the cusps of the epicycloid. 


At the cusps this peripheral stress is minimum and equal to half the maximum peripheral 
tress which oceurss at n=0 or n. 


Eingegangen am 9. März 1956. 


Application of Dirac's delta-function in isolated force problems 
of semi-infinite elastic solid of isotropic and non-isotropic materials 
By G. R. Verma in Pilani, Rajasthan (India) 


Probleme, bei denen Einzelkräfte auf die Grenzfläche eines unendlichen isotropen oder anisotropen Halb- 
raumes wirken, werden mit Hilfe der Diracschen Ö-Funktion gelöst. 


In this paper problems of isolated forces acting on the boundary of a semi-infinite solid, composed. of iso- 
tropic or non-isotropic material have been solved by using Diracs ö-function. 


Dans cet article des problemes de forces isoldes agissant sur les bords d’un corps semi-infini et compose de 
materiel isotropique ou non-isotropique ont &t& resolus en appliquant la fonction ö de Diracs. 


3anaya PaBHOBECHA U30TPOHHOTO HJIH aHHU3OTPOIHOTO IOJIyIPOCTPAaHcTBa ION MefiCTBueM 
COCPENOTOYEHHBIX CHI, IIPMJIOPKEHHEIX K TPAHHIE, pemaerTcH NpH IOoMomMm d-hyHkuam IInpara. 
Introduetion 


In the first part of the paper (Verma 54) problems of isolated forces acting on the straight 
unbodary of a semi-infinite plate were solved by using Dirac’s ö-Function. Solutions were 


zZ) / NE, A  , = 
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obtained in compact forms in the case of plates of both isotropie and non-isotropic material. 
In this part the corresponding problems of semi-infinite elastic solids with isolated loads on the 
plane boundary are considered. It is shown how Dirac’s ö-function can be used with advantage 
in the problems of this type concerned with semi-infinite solids of isotropic material as well 
as of those having transverse isotropy. For using ö-functions the solutions must be obtained 
in suitable forms and the deductions of such forms is the first object of this paper. Secondly 
it is demonstrated that the use of ö-function facilitates the solutions of problems when several 
loads act at different points of the boundary. 


Method of solution for isotropie semi-infinite solid 


Isolated normal loads, on the plane boundary of a semi-infinite solid r,,, 7,,, and o, will 
satisfy the Beltrami equations 
1 20 17&@ De 
2r 2 2 ==) 1.1 
re RE Ma re a en 2 


where 
02 SR, : ; 
Me at 32° 9=0,70,+0, n = poissons ratio. 
Also T,,, T,, and o, satisfy the equation 
Olgs | Mys , 90, _ 12 
ee Kr (1,2% 


Since © is a harmonic function, for given normal loads on the plane boundary, we can write 
the solution of the equation (1.1) as 


1 oo 1 00 1 


ae ar "sd wre 
where ® is a harmonic function such that 
[9-0 = [elk=o- 
Now we can write 
1 0N 
— EEE EN a ee 12! 
i 27 02 er 


where 


FF sus En dlm— n) de dn 
IE + NH 


and ö(x) is Dirac’s ö-function defined to be zero if x #0 and to be infinite at z= 0 in such 


a way that f ö(x)de =1. It also satisfies the condition that if f(x) is a continuous function 


Std @— ddr = fa). 


| >= 
_ — kan. 2 = RE re 1.5). 
Therefore N > ED + mM +zr® (1.5) 
Since © is a harmonic function we shall have from (1.2) 
924m,  9=214n9 


it being assumed that © vanishes as ® vanishes. 

From (1.4) and (1.5) > 

z 

® = z = m _213/9 9 

27 2 (zei re 
Ze 


m* 
6 Fe —— 


ST En 


0 1+7 | (2n—3)z a ns er 
3 Zen Eu lie at + Ft te 
Bien Z | en L 22 j 
 2En mE Er (m, Et za ee — 24m — N’ + 2RR 
3* 
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Method of solution for transverse isotropie semi-infinite solid 

In the case of solids having transverse isotropy stress strain relations can be written as 
nt my tm 
0, = Cp8&, + Cı &, FF 658%; 

0%, = let 8) + C;&5 


Te CaYer cha ll I ee (2.1) 
Tu: — (44 Yar > 
1 
Toy 3 (1 — Ce) Yzy 
The on of Eine Ü are | 
004.07, or OTzy DOHREOTE RE OS Dre n200, 99 | 
de au +3 Ze ler a Far RB | 
Also in the usual notations if u, v, w are displacements 
= = arg 2 ze = =” ek, N 


from (2. ni . 2 (2.3) we obtain 


0° ae 
Cı — le a dye 2 ie 44 5 Fra SC * a Sch (O3 + 64) »|- 0a (2 
1 0®v 0? 2: oI1 ou w 
5 (Cı — €2) ar? -F vor ur 44972 war en (Cu + C1e) Ey + (Ge + Cu) 2 =0 ...(25) 
ww Ow ov 
A| nt an) Bl tt ee es N ige woiszeat ser 0, WE (2.6). 
(cf. Elliot 1948). 
Assume 
od od od 
u= dx oy > —= KK a a Ve a a A (2.7) 
Then the equations (2.4) and (2.5) are satisfied if 
dd 0 
Cı dx2 - oy? 4 [eu - K (Cı3 + a = — 0 aa a 0 ae (2.8) 
and (2.6) is satisfied if 
dd Hd od 
[3 + Ca + Kca] + =) +K EP =D... ee (2.9) 
Both (2.8) and (2.9) are satisfied when 
= F Eu eu Sr (2.10) 
Cu Kcua+ (Gs + Cu) 
Then y, and y, are the solutions of 
Cı Ca Y” + [Cıs (2 Ca + C13) — C1g Cl Y + Caa Ca = 0 
and (2.8) and (2.9) take the form 
n 0° ? 
[+ S A) MER .\ (2.11) 
where 
0° 0° 
2 Sure = 
Fa dat oy% 
So 
od, 08, 08, , 0b, oD od 
St en Te 


TOTEN er 


K, and K, are values of K corresponding to the values y, and 3 of y. Then stresses are obtained as 


8D, 
6, = (05, or? => 6293) 0, + ®,) + fr 1 = 2 Zr RK, = >) ee (2.15), 
D, RD 
0, (u. or: + Cı ay? )@ + ®,) -+ C1g En za P) 2 RK: ns ae. (2.14), 
[) D, ®. 
09, — (ig = % r a) (8, Hr D,) = ou (K al Fr3 2 22 168 Dr OO (2.15), 
Ty> (u — &n) = Z (Dia) ee ee (2.16), 
2 a ae ee (2.17) 
Tut + at a Eee oo e B 
REN fü + a 3% a 2+-(+K) a | ent heit ou (2.18). 
08, 0D, 
ek 
en 
Ta le ee a a I Be er er (2.19), 
oV. op. 
N, a +K,) Fe (0.2.8) — Ve. (2.20), 
EN oY. 
Tr tu fü - K,) öy (1 -H-.K,) a RE REN Een e2p% 
From (2.11) 
ad = D, od, 
a ar FE ei: 
and 
020, 0, 0D, 
a N 
so 
0, = (C3; Ki — 63 ws Fi (C33 Kg — € 92 —: a net all na sag: (2.22). 
Now suppose for a load P, at (&, n, 0) 
oo oo l 
De (17-6 K;) | | o en ER IR ae NEE (2.23), 
—00 —o0 ; 
7 MK,) | | ” Fe EEE IN, (2.24) 
. . 2 
where 
2 z? 
Ne ee N Nike 
Yı Ya 
r, and r, are equalonz=0. 
putting 2, = © and Zr we ‚ we get 
Yı Yy Ya 
P, € 1 or, 99m 
0, = (ig Kı — Cs ı) > DE = + (0, Ks 05 Y) - at Fe ee (2.25), 
1 = 
[0,]e-0 = — Pı = (@s Kı — © Yı) 2 (1+K)(C-2r0) 


I, : 
+ (C; Kg, — Cs Yo) Y a Rn (2.26), 
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Ze 1. Adrenalin nn einen nn 
P, 
o= 1 1 = 1 ’ 
27 |(&m K,— sy) + K)— — (as K,— c(3y) (1 + Kı) | 
Yı 2 
oo oo 
d+R,) >f. [Pd 
iz 1 1 na 
20105 Kr — Ca Yı) 1 + Ko) — — (is RK, — 3 y)(l+K&ı) | en 
Yı Ya 


If there are several loads P,, Py, Pa: P„ etc. at (E71 0) (&a 7% 0) * * (Ems Nm 0) ete. 


Then 
1+&) 
en ve 1 
27 @ K,— 13yı) L+Kı) = — (O3 Ka — Ca y)(A+Kı) ) 
2 \ Ya 


2 nt 


N 
2 | jb> Im 


09 .—09 


(+&) n_ P, i j 
2 I 1 > TE | 
22 (en K,— sy) 1 + K)) on — (gg BR — sy) (l + K) eu m=1 
1 


Yy. 


Similarly 
= (R,) np, 
= 1 1 L, ’ 
2 (00 KR, — sy) + K,) Ya, —(g RR — y)(l+K)) zit m-1 
BR TI R,T, 
Kı,—K, n io 
2 1 1 TI 
27 | Kı— ey) + RK) — — (3 KR — ca Yy) (l + K,) —tmeı 
| I Yrı 
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Die dicke Platte mit nichtlinearer Spannungsverteilung 
Von Eberhard Koppe in Göttingen 


(Aus dem Max-Planck-Institut für Strömungsforschung, Göttingen) 


Am Beispiel der dicken Platte mit nichtlinearer Spannungsverteilung wird dargelegt, wie man ausgehend 
von den Invarianten des Problems einen bequemen und allgemeinen Weg zu strengen Lösungen findet. An 
den Seitenrändern der Platte können dabei unabhängig voneinander Normal- und Tangentialkräfte, Rand- 
scherkräfte, Randbiege- und Randdrillmomente vorgeschrieben werden. Für zwei sehr allgemeine und prak- 
tisch wichtige Klassen von Belastungen werden die Verteilungsfunktionen für die Verschiebungen über: die 
Plattendicke ausgerechnet. 


On example of the thick plate with non-linear stress distribution, it will be shown, how to find an easy and 
general way to exact solutions starting from the invariants of the problem. On the boundaries of the plate 
we can prescribe normal und tangential forces, shearing forces, bending and torsion couples, independently of 


one another. The distribution functions for the displacements over the plate for two general and practically im- 
portant classes of loads are calculated. 


Au moyen de l’exemple de la plaque Epaisse avec une distribution des tensions non-lin£aires il est expose 
comment on trouve un chemin confortable et general menant ü des solutions exactes en sortant des invariantes du 
probleme. Independantes les unes des autres des forces normales et tangentiales, des forces de cisaillement margi- 
nale et de torsion marginale y peuvent &tre prescrits aux bords lateraux de la plaque. Pour deux classes de 


charges bien generales et pratiquement tres importantes les fonctions de distribution pour les deplacements sur 
l’epaisseur de la plaque sont calculees. 


Fla IpnMepe-INTacTuUHkH KOHeyHoü TO 


2 IIMHBI C HeJIHHeÜHBIM PAacıpepnerreHnneM HaNpsBKeHnü 
N3IATAeTcA YNOÖHBIH MH O6IMMÜ MeTon, 


NO KOTOPOMYy, UCXONA u3 HHBaAPHAaHTOB 3anayH, MO3KHO 


za 


ee 


lautet in Komponenten 


VE HSEHa nee 


it 
0 


Be= (Die partiellen Abteilungen werden durch |, oder |? bezeichnet ; bei nichtkartesischen Koordi- 
naten kann man |; als kovariante Ableitung auffassen.) 


; ni Dabei ist ® eine Invariante; die a, sind die Komponenten eines Vektors; e?*! ist ein voll- : 
a kommen schiefsymmetrischer Einheitstensor (allgemein mit den Werten 0, Si, = en ‚ wenn g a 
die Determinante des Maßtensors g;;, der Koordinaten ist). Yg Vs 

A e . = 

F E 


Bild 1 


Wenn man von Grenzfällen wie „Quader“ oder „Streifen‘‘ absieht, so ist bei der Platte eine 
Koordinate — die z-Richtung (Plattendicke) — gegenüber den beiden anderen Koordinaten — 
x- und y-Richtung (in der Plattenebene) —, die untereinander gleichberechtigt sein sollen, aus- 
gezeichnet. (Vergleiche Abbildung.) Führt man diesen Gedanken durch, so folgt 


a,=d,, u Ed 


4 

ri 

; Dabei sind yund Y Funktionen, die gegenüber Transformationen von x und y invariant sind. &0 
: sind die Komponenten eines vollkommen schiefsymmetrischen zweistufigen Tensors. Setzt man 
} zur Vereinfachung der Schweibweise 


i 2 ee 2 rt und („-pP)=R, 

4 so erhält man schließlich 

vi = Dj: +. > ezik R A er Y70 

oder bei Ausführung der Überschiebung der e-Tensoren 

"Dit +23 RL SH IR re EN) 


“ mit den Kronecker-Symbolen ö,. 


2. Gleichgewichtsbedingungen 


Bei der dicken Platte kann man nicht mehr annehmen, daß die Spannungen linear über die 
Plattendicke verteilt sind; man kann jedoch voraussetzen, daß die Durchbiegung so gering ist, 
daß die Gleichgewichtsbedingungen am unverformten Körper erfüllt werden können. Ferner soll 
der Zusammenhang zwischen Spannungen und Verzerrungen und Verzerrungen und Ver- 
schiebungsableitungen linear sein. Unter diesen Voraussetzungen erhält man für die Gleich- 


gewichtsbedingungen ohne Volumenkräfte 


RE En) ET N TC, 


M fü 
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und für die Oberflächenbedingungen 
. Ä | au. B Re a) 
zen lo He + 13, )nd)= Pi EEE (3) 


Bei Einsetzen von (la) und mit den Symbolen + %=4A und A+|= D folgt aus den 
Gleichgewichtsbedingungen 


2U—M Joe + OR Ye=0, 

1-3, PP’ + OD RArO k 

Mor —ORL+ BEAAR Een (4). 
ae, 

u or — DAYF=0 

1—2u 


Es sei angenommen, daß die Funktion R, von ® und Y unabhängig ist und der Dif- 
ferentialgleichung 


RSS Arne Tr (5) 


genügt. Auf diese Weise werden aus der Lösungsgesamtheit Klassen von Lösungen ausge- 
sondert, die sich im Rahmen des hier behandelten Problems als hinreichend allgemein heraus- 
stellen werden. 

Zu der Differentialgleichung (5) treten die aus den Randbedingungen ?” = 7’ = 0 bei 
z= + h folgenden Beziehungen 


RE=R&k=0 für z=+h. 


Macht man für R, einen Produktansatz R, = S(, y) T(z) und berücksichtigt die Randbe- 
dingungen, so folgt aus (5) für den mit der Biegung der Platte zusammenhängenden Teil von R, 


Re Sein Ct Dn. |: 5.@, Re (6), 
(ungerade) v=0 2 h 
wobei jedes S, der Differentialgleichung 
2v-+1% 2 = 
AS, | 3 ) N er er (7) 


genügt; entsprechend ergibt sich für den mit den Dehnungen zusammenhängenden Teil von AR, 


R, = 3) cos kn: 2 SS LDEW)! 0m en re ee (8) 
(gerade) „=1 h x 
mit 
u 702 
AS, => m . Sr = 0 ei Pehlan et, other. are BEER (9. 


In den nach dem Ausscheiden von R, noch verbleibenden Gleichungen (4) lassen sich die beiden 


ersten Gleichungen nach x bzw. y integrieren und man erhält — bis auf eine Funktion von z, 
die gleich Null gesetzt wird, — 


All 1) e 2(1—W : z 
ee ui FR: Ta = 7A 0,2 sr 


3. Die Invarianten des Problems 


So wie in den Ansätzen für die Verschiebungen die Tatsache, daß es sich um eine Platte 
handelt, bei der die z-Richtung gegenüber der x- und y-Richtung ausgezeichnet ist, wesentlich 
benutzt wurde, wird nun bei der weiteren Behandlung des Gleichungssystems (10) beachtet, 
daß ® und Y nicht irgendwelche beliebige Funktionen sein können, sondern aus den Invarianten 
des speziellen Problems zusammengesetzt sein müssen. 

Diese Invarianten werden aufgesucht: 

Eine Invariante der Koordinaten gegenüber Transformationen von x und y ist trivialer- 
weise der senkrechte Abstand der Punkte der Platte von der Mittelebene. (Bei rotationssym- 
metrischen Problemen würde man ?= x? + y?, das Quadrat des Abstandes eines Punktes 
von der Mittelachse, hinzunehmen.) 

Die Durchbiegung w* = w*(x, y) der Mittelebene ist eine Invariante der Verschiebungen; 
dagegen sind die Verschiebungen u* und v* in der Mittelebene keine Invarianten. Schreibt 


j 
s 
N 
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man jedoch 
u*=o,+ol,, rt ee re um) 


so sind 9 und & Invarianten, aus denen sich die übrigen Invarianten der Verschiebungen in 
der Mittelebene, z.B. u*|, + v*|, = Ag, die mittlere Dehnung, einfach ableiten lassen. 
Invarianten der Kräfte sind die Kraft senkrecht zur Mittelebene Pa = Aw* und 


die mittlere Spannung S 7£ — Ap in der Mittelebene. Zu diesen reinen Pete one 
k 


treten die Invarianten der äußeren Belastung, die bei senkrechter Last p(x, y) die Form p, Ap, 
A? p,..., allgemein A”p haben. 

Gestützt auf den linearen Fall, wo u=u*—zw*,, v=v*— zw*|, w = w* gilt, 
wird mit diesen Invarianten angesetzt 


DB +PE=— zw + fe) Aw + Ir,@)A’p+P+ me) A’p 


8, — AP} — mw* u 9(z) - Aw* -- St 1,(z) 4’p = >= n,(z) A’p a (12). 


Dabei sind f(z), r,(z) und n,(z) ungerade Funktionen g(z), 1,(z) und m,(z) gerade Funktionen. 
Für die Invariante » aus (11) gilt 


mit S, aus (8) und (9). Alle anderen in den Differentialgleichungen oder Randbedingungen 
benötigten Kombinationen von ® und Y lassen sich aus (12) leicht berechnen. 
4. Ein rekursives Gleichungssystem für die Koeffizientenfunktionen 


Trennt man nach geraden und ungeraden Funktionen von z, so erhält man aus (10) die 
beiden Gruppen von Differentialgleichungen 


ra 20) F— 2 — 2] Aw* + 2(1— u) f> AAw* + Inn d—2m) + Llp 


a ne en a (14) 
und 
2A: g’—2u] Awr + [ff + A— 2) gl Adw* +2 (1— m) p 
+ 20 — ba tn +L,d4—2m]4tHp=0 ....(05) 
sowie 
2(1— u) Ap + [l—2u) mo + no] p 
I malern) +2d—e) m +) dp (16) 
v=0 
und x : 
Spin tmel-2u)n)dp 02m 
v=0 
Aus den Randbedingungen in der Form (3) folgt aus 
N) fürızseh 
und 
Te—! für! zrHh, 


ng 
wenn in gerade und ungerade Funktionen in h getrennt wird 


FEN In +) Ap=0 für zehn... (18) 
und = 
d- Mg —uz]4w* + nf: Adw* + l— ml p+ 
u Klo) bar ar] Akt p u un nd p three eh. 23 us (19) 
nr Y” (m, +n)2’p=® a (20) 
und eE 


np + l—mnep +2 [lu a 


49 Koppe, Die dicke Platte mit nichtlinearer Spannungsverteilung »4.37 a ee 


Aus den Differentialgleichungen (14) bis (17) und den Randbedingungen (18) bis (21) 
sollen nun die Verteilungsfunktionen f(), r,(@), n,(2), 9(2), 1,(2), m,(z) bestimmt und je eine 
Differentialgleichung für w* und p gewonnen werden. Dazu wird für die Differentialgleichungen 
der Ansatz 


AAw* + 2 C,4p=0, Zo+3,1P =0,,. re = (22) 


mit vorläufig unbestimmten Koeffizienten C, und c, gemacht und in die Gln. (14) bis (17) sowie 
in die Randbedingungen (18) bis (21) eingeführt. Auf diese Weise erhält man, wenn man nach 
den A’p ordnet, ein rekursives System von Gleichungen und Randbedingungen für die Ko- 
effizientenfunktionen f, 9, r,, I, m,, n,, das gleichzeitig die Konstanten C, und c, liefert. 


4. Der Spezialfall AAp = 0 


Löst man dieses rekursive System für den speziellen Fall AAp= 0, so ergibt sich für 
die Verschiebungen 


u=o,+ sin 5 Solyas cos >. as —2e Wr, + E j BE — | Ar, 
le en era + ze pl. 
a | re BB a 
Han EB -Apl | 
entsprechend 
=, sin, Sc cos za S,— zw, +[ Amt 
ee 1 Apl, - Sun See er (24) 
und 
ww re een h® x ei 
en ee dam ce ae Ham +28] Ap (3). 


Dabei müssen die Funktionen w*, S,, p und S, den Differentialgleichungen 


em are 
ae ne 
AAw* = Gm 10Ch DE EI (26), 
a 
AS — a So 0 u ce a Dura de Bahn ee ee a (27), 
_ ren 
Jo = ERST AD Per or Er (28), 
7a 
AS, — m =0. 0. Vi (29) 


genügen. 
An den Seitenrändern der Platte sollen fünf unabhänsi i 
z t gige Randbedingungen gestellt 
werden. Es genügt daher, von den Funktionenscharen S, (6) und S, (8) nur eweils aie er 
Funktion mitzunehmen. Mit Hilfe weiterer Funktionen der Scharen S, und S, könnte man 
die Bedingungen an den Seitenflächen der Platte noch schärfer fassen und etwa vorschreiben 


wie sich die Kräfte oder Momente in z-Richtung verteilen. 
6. Vergleich mit einer Näherungstheorie 


Es ist interessant, die obigen Ergebnisse der exakten Theorie fü 3 
“E nt, obig orie für den Fall AAp = i 
denjenigen zu vergleichen, die sich aus der von M. Schäfer [1] vereinfachten und 


a ee) a ne a I a nn 


ee E ae 
7 


ME ea | 
e Platten ergeben. Hierzu werden die Gleichungen 


a Dt N en 
RE ee N | 
pe Ha 874,3 

ei | 


70% # 
AS, gg Ser ljol= 


geschrieben. Die Näherungstheorie [1] unterscheidet sich durch die eingeklammerten Faktoren 


BrDie Übereinstimmung ist praktisch vollkommen; zumal ıw* hier direkt die Durchbiegung der 


‘von der exakten Theorie: 
E a 1! 


ei 0 Vers 
ee ne hen in] > 09- 


4 
2 


Mittelebene der dicken Platte ist, während in [1] w eine mittlere Durchbiegung der dünnen 
Platte bedeutet. Natürlich bedingt die Berücksichtigung der Plattendicke Änderungen in den 


_ von z abhängigen Koeffizientenfunktionen. Da jedoch diese Verteilungsfunktionen unabhängig 
von den Funktionen in x und y sind, kann man die in [1] und nachfolgenden Arbeiten erzielten 
- Lösungen für die dünne Platte auch für die dicke Platte nutzbar machen. 


- Über die Dehnungsgruppe werden in [1] keine Aussagen gemacht. Die Differentialglei- 
chung (28) zeigt, daß im Fall Ap = 0, d.h. bei relativ gleichmäßiger Belastung, die Dehnungs- 
größen lediglich über die Querkontraktion etwas von der Belastung p spüren. 

Eine Näherungstheorie, die im wesentlichen zu denselben Resultaten wie [1] führt, wird 
in [2] gegeben. Die Funktionen S, und S, aus den Rotationsanteilen im Ansatz für die Ver- 
schiebungen werden in [2] über den „zur Mittelfläche antimetrischen Spannungszustand einer 
eben bleibenden Platte‘ eingeführt. 

In [3] wird darauf hingewiesen, daß man zu den bisherigen Lösungen des Plattenproblems 
Glieder hinzufügen kann, die den R,-Gliedern (6) entsprechen. Wie oben gezeigt wurde, liegt 
der Grund hierfür einfach darin, daß die bisherigen Ansätze unvollständig waren. 


7. Der Spezialfall Ap = — y?p 
Der im vorangehenden behandelte Fall AAp = 0 ist ein Spezialfall von A"p=(0. Da 
die Bestimmung der Koeffizientenfunktionen und Konstanten endgültig ist, lassen sich weitere 
Funktionen und Konstanten leicht berechnen. 
Es sei nun noch kurz auf Belastungen eingegangen, bei denen Ap proportional zu p ist; 


zum Beispiel Ap=—-y?p. Ein Beispiel für diese zweite, ebenfalls sehr allgemeine Klasse 
von Belastungen ist die Rechteckplatte (Kantenlängen 2a und 2) mit doppelter Cosinusbelastung 
ER ae 


Hier ist Ap = — y? p mit 
72 
Wagener D?). 


Im Falle der Proportionalität von Ap und p kommt man mit den Koeffizientenfunktionen 7, 
I, My N, allein aus. Die Differentialgleichungen für w* und o vereinfachen sich entsprechend. 
In den Verteilungsfunktionen treten jetzt neben Polynomen in z auch hyperbolische Funktionen 
auf. Es ergibt sich z. B. 


AT 3 h (2 —.u) 2 le 
mp), Hein ES] Foont Se wil + |? A a) Aw*|, 
DEE I ZE 2(2 — u) zn 
726 [Cin2yh—2yh] »6ll—a) dp) 
N EEE Ä 
* y„2.2G[©in2yh+2yhl Plz (80) 


Sinyh— (1— 2u)Cojy hl Sinyz—yz:Cojyh-Colyz 
Ei 2Gy?[Sin2yh—2yhl 
[1 — 2u) Sinyh—yhbojyh] boijyz+Yy2 Ginyh- = ee 


IT 


= 26 [Sin2yh+2yhl | 


44 Voditcka, N  .. in zusammengesetzten Stäben na.gr me. 112 Tan fwebr. 1057 
ET h 
oder Sofyh-+ De 
[% 2 A—w)Eoiyh+Y ge a De ER 
venttAgg Anis, ToGlemeyn- 2 ee 
Ib h&inyh+2(1—w)6Cofyhl&ofyz—yzCojyh- Oinyz 1). 
T 2Gy[Ein2yh—2yhl 


2 (1 —u) Sinyh-+yhbojyhl- Sinyz—yz Surnlya, 
oe 2Gy[Sin2yh+2yhl 

Die Funktionen w*, S,, 9, S, müssen dabei den Differentialgleichungen 

._YPA—WBonh TREE 32), 


Adw 36[&n2yh—2yhl 


2 » 
AS — 50 =0 ne BE er ne JEBhe 
Eat [1—2u) Sinyh—yhbojyhl (34) 
Im mare 
1, 55-0 a a (35) 


genügen. 5 nr 

Entwickelt man in (32) und (34) die hyperbolischen Funktionen nach Potenzen von y Ah, 
so folgen für die beiden ersten Glieder die Gln. (26) und (28). 

Mit Hilfe der Partikularlösungen für (32) und (34) 

> — h h&inyh 
5 A RS ee 
y-2G [Sin2yh—2yh] 
-  [d—2u) &inyh—yh6boyhl, 
0 9%-2G[&in2yh+2yh]l 
lassen sich die Ausdrücke (30), (31) vereinfachen, da dann gerade diejenigen Glieder wegfallen, 
die hier infolge von Partikularlösungen der rekursiven Gleichungen für r(z) , /(z) und m(z) auf- 
treten. 

Die dargelegte Theorie der dicken Platte umfaßt die bereits länger bekannten Spezial- 
fälle, etwa die in [4] behandelten einfachen Beispiele konstanter oder linearer Belastung, weiter 
die ebenfalls in [4] angegebenen ‚‚verallgemeinerten ebenen Spannungszustände“ oder die in 
[5] behandelten Partikularlösungen. 


Ges 
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Biegungsschwingungen in zusammengesetzten Stäben 
Von Väclav Vodicka in Pilsen 


Mathematische Untersuchungen von Biegeungsschwingungen gleichförmiger zusammengesetzter Stäbe führt 
zu recht komplizierten Rechnungen. Eine Vorstellung davon gibt der folgende Aufsatz. 


The mathematical investigation of bending oscillations of uniform composite rods leads to rather complicated 
caleulations. This is illustrated in the present paper. 


L’examen math£matique de vibrations de flexion de bätons uniformes et composes aboutit ü des caleulations 
compliquees. L’article suivant en donne une idee. 


MarTemarnyeckoe HCcJIeNoBaHne MOINEePeUHbIX KOJleÖaHHÜ COCTABHLIX 6aloR IpPUBOMUT K 


MOBOJIBHO CJIO}KHbIM PACUETaM. HekoTopoe IIpelCcTaBlIeHuHe 060 3TOM Naect Tpemnaraemarı 
CTaTbA. 
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1. Allgemeine Theorie 

1. Ein gerader gleichförmiger Stab von der Länge I =1, +1, + :-- + I,, dessen Querschnitt 
die Fläche S und das Trägheitsmoment I hat, setzt sich aus n Teilen von verschiedenen Materialen 
zusammen. Der k-te Teil hat die Länge /,, den Elastizitätsmodul E, und die Dichte o;. 

Wir betrachten den Fall, wo die Stabenden x = 0, x = ! frei unterstützt sind und weder 
Einzelkräfte noch Einzelmomente an den Berührungsstellen zwischen irgend zwei benachbarten 
Teilen vorhanden sind. Die äußere auf die Längeneinheit bezogene Kraft sei f(x, i), die Anfangs- 
lage und Anfangsgeschwindigkeit des Stabes sei durch u(x), bzw. v(x) bestimmt. 

Schreibt man allgemein y = y(z, t) für die Ausweichungen des Stabes aus der Ruhelage 
(d.h. von der a) und setzt noch 


=, u, Ve Pe MER Er (1), 


z=1 
so geht es um die Lösung des Differentialgleichungssystems 


9, 1 ®y fr ee | 
— 0: = 1 (k=eh2e, 2 
ST RE eh en ge 
mit den Anfangssbedingungen 
9y,.(&, O 
Y.(2, 0) = u,(X), en Ev) ar kkelae..n, 00H 
und mit 
9°y1(0, 
a a a EEE ce 3.2.1) 
yo = 1> 8.2.1) 
2 Sy 9% (5, Oyılsı De Oyr+ı(Sp D) 
Y,(5,D) = Yr+1(5 D, a u = du? e Br (32:2) 
Oyr(sı i) Di Oyrrı[S% D) Er+1 N a a 
0x3 er u E, en 
®yull id 
y.d= er OD 0 ee (3.2.3). 


Die f; = fx(&, D), uy(x), vy(x) in (2) und (3.1) sind natürlich die auf den Bereich 51 <r<s; 
entfallenden ‚Teile‘ der Funktionen /(z, £), u(x) und v(e). 


2. Setzt man allgemein 


Mep)=njetine dd, Fii,D)= 5 gErofet)dir aa (4), 
ö 
so folgt aus (2)—(3.2.3) das gewöhnliche a 
d2Y Fy(&, P) 
Er ey 2 Bud =, Ha<ı<yaı = Pi, p) = et 7 a ua) v7 3 il) 
2) re ae See ee (5), 
BES IE EV (ON EB ne EEE Be (6.1), 


Yı(»P) = Yerılsp P) » Yı(sp» P) = YeıD) Y% (S% P) 
2,20. Ye. D- ar De 1,2 2 nm 1) 16.2); 
El a ea ee ae ae ie (6.3). 


Striche bedeuten Ableitungen nach z. 
q Mit der Abkürzung 


je OR N, STD. (7) 


I 
ergibt sich aus (5) leicht 
Ye, p) = E ch > x + B,sh 204 C, cos Ay x + D,sin nn Zr iS pP)» 


ZIP) > a D,(&, p) [sh A, (e — &) — sin A, (@ — Sr u enasse.n) 
+, 


Daher ist = x h 
Ze» PW= Zus, P) = Zee De > Wear. an el) 


Aus den Bedingungen (6.1) — (6.3) erhalten wir nun leicht folgendes Gleichungssystem für 


die Konstanten A,, By, Gy, Dr: 
ee ee ee (10.1), 


Arzı ch sy Aryı + Bari SR 5, Ap4ı + Cr+ı 608 5, Ag+ı + Dir sin s Ap+1 ) 
z A, ch Sp Ar + B, sh 5. Ar -H C; cos 5. Ar + D; sin 5, Az + Z(5% D) 
Arsı sh 5% Arc+1 + Br+1 ch Sr Ac+ı 47 Cr+1 sin 5% Ak+1 E= Dr+i COS S; Ar+ı1 
=h, IA: sh 8,A, + By ch 5, A, — C, sin s,A, + DyC0s 5,4, + le p)) 
k 


Aryı ch 5, Antı + Ba+ı SA 5, Aryı — Cr+ı 608 Sp Ag 4 ı — Der SiN Sy Ag+ı 


hz 
= Fr m ch s,A, + By sh 5,4, — €, 608 5. A, — D; Sin 5; A, +5 : „7 6 s»P p)| ur 
k 
Ayrı sh Sk Ar+1 + Br+ı ch Sk Ak+1 + Cr+1 sin 5% Ak+1 —— Di+i COS Sy, a 
BR 
= la „Shs,Ar . B, ch 5,4, + C,sin s, A, — DyC08 5,4, + F 1 ze, p)|; 
); 'a E; 
ELBE ee =| &r "er je er 
EURE = Ve E, ( 1,2,....,n=1) 
A„chlA,„ + B„shlA, + C„cos!A,+ D,sin I), = —Z,(l,p) 
1 RSS] 


A„chlA, + B„sh1A, — C„ cos 1}, — D, sin 1), = Zul P) 


cin 
3. Unsere Aufgabe reduziert sich so zur Auflösung der Gleichungen (10.1)—(10.3), was 


wahrscheinlich am besten in der Matrizensymbolik geleistet werden kann. Wir gehen aus von 
den Beziehungen (10.2), die sich nach einer längeren, jedoch leichten Rechnung in die Form 


KM DD en. (663) 


setzen lassen. Hierin ist 


Kö ENTER EEE (11.1), 


k 
die Elemente ml, (r, s=1,2,3,4) der Matrix M,(p) berechnen sich nach der Vorschrift 


mi = (1 4 =) (ch 5.4 ch 5; Ar+1 —h, sh S. Ar sh 5, Ak+1) 


Rz 


e 
m) = ( 7) (cos s, A, ch S.Ap+ı + Äh, sin Ss, 1, sh S, Ar+ı) 
(& Bu\S® 
m» = \1 — 2x (sin 8, A, ch 5, Apr +ı — hy COS $, A, sh 8, Ax +1) 
ee: 
n®—=\1l+ & (— ch 5, A, Sh Sy Ay+ı + hy sh sy A; ch 8, Ar) 
m = (1 + sh sa, shsyA + h,chs,A, ch s, A j 
B e k kAR+1 x Ch S, A, ch 5, Ar+ı) 
(k) — |] hi 
My; — > &; (— COS 77 Ar sh Sr Ar+ı en h; sin Sr A ch Sk Ar+1) 
an hi h 
m) 1 = (— sin s, A, sh S, Apr +ı + hy COS $, A, ch 8,441) 
hi 
m) — ( — e) (ch 5, A, C08 8, %x4ı — h, sh s, A, sin 5, Ag 41) nn 
me— 1 hi 
rl (sh 5, A, C0S 8, Ax4ı — h, ch s,A, Sin S, Ag +1) 
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Un 
m ( -4+ 75 (cos 5. A, COS S, Ar+ı + h, sin s, A, sin Ss, Ar +1) 
k 


REN: 
I y— i + (sin Ss, A, COS $, Arıı — A; C0S 5, A, Sin 5, Ar+1) 
k 
E hi 
m a (ch s,A, sin Ss, Aptı + hy sh Ss, A, COS 5, Ar +1) 


(sh s, A, sin Ss „Ar+ı = h, ch s.A, C0S 5, Ar +1) 
Br hi 
m®—=\1+ Er (cos Ss, A, sin $, Ar ıı — hy, Sin s, A, COS 5, Ax+1) 
k 
s Be 
md—|1-+ ei (sin Ss, A, sin $, Apıı + I, COS S, A, COS 5, Ar+ı) 
k 
(= 1,2, MArn- 1) 
und die einspaltige Matrix P,(p) hat die Elemente: 
1 hi 
pi = 32% DS p) (1 ur =) [ch 5, A +ı $h Ay (85, —&) — hy Sh Sp An+ı ch Ay (5, — ES) 
. k 
1 


hr 
— ( _- ee) [ch 5 Ag ıı sin A, (5, — &) — hy 5h 5, Aryı 608 A, (8, — 5) de 
k 


5% 
1 hr 


%k—1 


hi 2 
— | — “) [— sh 5, Ag yı Sin A, (5, — &) + hy ch Sy Ar +ı 608 A, (85, — 91) de 
k 


&% (11.3). 
 — 5% | D,(E, p) (1 -) [cos 5; Ar zı sh A, (5, — €) — hr Sin Sy Ay rı Ch A, (5, — O)] I 
Sk—1 
—_ ( - na [cos 5; Ag +1 sin A, (8; — &) — hy Sin 5; Aprı COS A, (S, — 5) dE 
5% 3 
Di — 5 | ,(& p) (i 2 a) [sin sy Aysı shAy (4 — + hy 008 Sy Ayzı ch Ay (5 — Bl 
®k—1 
ir ( 4 ) [sin 5; Ag sin A, (5, —&) + hy COS 5, Aryı COS A, (S, — 9) de 
a ) 
Dazu ergibt sich aus (10.3) 
1 
A,„chlA,+ B„shl}, = an IEZE p) shA,(l— &) dE 
RENTE a hr m (11.4) 


l 


1 
C,„ cos 1A, + D,sin 1A, = 958 | D,(&, p) sin A, (L— 8) dE 
Br Sn—1 
und unser Problem besteht so im wesentlichen aus der Berechnung der Beiwerte A,, B;, Cy» 
D, aus (10.1), (11)—(11.4) und Einsetzung in (8). Jedoch ist auch dies noch recht umständlich. 
Zuerst ergibt sich mit der Bezeichnung 


N%» = N%» (p) = M,(p) - Mi-ı(p) Mı-s(p) * M«(p) 
Bares rk, nl) Kalt sense (12) 
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ID.| 
g&(p) = P(p), Q®(p) = 3 2-1 N%=+D Pp) +1 Pıp) (k=2%3,...,n—1) (13.1). 
x=1 


Drücken wir nun A,, B,, C„, D, mit Hilfe von (13) durch B,, D, aus und setzen dann in 
(11.4), so kommen wir zum Gleichungssystem 
B, [n®-%D ch 1A, + na sh 1A,] + Di [Ind ch 1A, + ng "Dh 14,] 
1 


aus (11), (10.1) und (11.1) leicht 


0 
2% K,,, = N%) fr 


n—2 
= ed hl, do | Dun) hr, d— He 


= | (14). 
B, [n®—1D cos 1A, + nid sin 1A,] + D, [ni cos II, + ng sin 1A,] 


1 
Nn—2 
— — ga) cos1A,— gr sinlA, + —_ | D„(E, p) sind, I— de 
In—1 ; 
Die n'y »», 4%” sind dabei die Elemente der durch (12) und (13.1) bestimmten Matrizen 
Na, 1) a 08H : 

Aus (14) sind nun die B,, D, in jedem Einzelfall zu berechnen und sodann in (13) einzu- 
setzen. Theoretisch ist damit — s. die Formeln (10.1), (14), (11.1) und (13) — die Bestimmung 
der Beiwerte A,, B,, C,, D, zu Ende geführt. Ihre Einsetzung in (8) gibt uns die Laplace-Trans- 
formierten der gesuchten Lösungsfunktionen y,(&,D); k=1, 2,...,n unseres Problems (2) — 
(3.2.3). Die Lösung selbst findet sich daraus nach der bekannten Theorie der Laplace-Transfor- 
mation. 

4. Wenn auch unser Problem durch vorige Erwägungen theoretisch als gelöst angesehen 
werden darf, so ist doch sowie die wirkliche Durchführung als auch das Ergebnis der Rechnung 
in allgemeinen hoffnungslos kompliziert. Der Kern und zugleich auch der schwierigste Teil des 
Problems liegt in der durch (12) angedeuteten Multiplikation der Matrizen M,(p). Die Auf- 
stellung irgend einer geschlossenen Formel (analog wie es der Verfasser in seinen früheren Arbei- 
ten [1], [2] getan hat) ist hier sogar in verhältnismäßig einfachen Spezialfällen h,=1, h; = e; 
(k=1,2,...,n—1) ganz unmöglich. 


2. Der Fall eines zweiteiligen Stabes r 
1. Obwohl unsere vorige "Theorie im wesentlichen nur unüberwindbare Kompliziertheit 
des in Frage stehenden Problemkreises zum Vorschein bringt, so ist sie doch keineswegs zu ver- 
werfen. Sie gibt je dem Techniker und Physiker fertige Formeln zur direkten Behandlung wich- 
tiger Einzelfälle, wie man sie in der Praxis antrifft. Als einfaches Beispiel soll hier mindestens 
der kaum in der Literatur zu findende Fall eines zweiteiligen Stabes betrachtet werden und zwar 
unter spezieller Annahme h? = e, (was übrigens nicht wesentlich ist). 
In diesem Falle haben wir mit nur einer Matrix NÜ,D(p) = M,(p) und einer Matrix Q®(p) 
— P,(p) zu tun, so daß die oben erwähnte schwierige Multiplikation fortfällt. Die Elemente 
mp; = m,, (,s =1,2,3,4) von M,;(p) vereinfachen sich jetzt:nach (11.2) in folgender Weise: 
Mg; = Mg, = My = Mg = My = Mg = My =Mg = 
My = 2 (ch s, A, ch s, A, — h, sh sı A, sh Ss; A,) 
My = 2 (sh s,A,chs Aa, —hıch s, A, sh s, A,) 
my = 2(—chs A, shs/A;, + h,shs, A, ch s, A,) 
Mg = 2 (— shs A, shs,A, + h,chs, A, chs,A,) a N“ 
My; = 2 (C0S 5, A, C0S SA, + hy Sin s; A, sin s; A,) 
My = 2 (sin Ss, A, cos Ss, A, — hy cos s; A, Sin Ss; As) 
My; = 2 (C0S S, A, sin $) Ag — hy Sin Sı A, cos S; As) 
my = 2 (sin Sı A, sin Ss Ag + I C0S Sı A, CoS Ss; A,) 
Was die zweite im Problem beteiligte Matrix P,(p) betrifft, so hängt diese nach (11.3) und 
(5) mit dem Anfangszustand und der Belastungsart des Stabes zusammen. Um etwas Bestimmtes 


vor Augen zu haben, betrachten wir den häufig vorkommenden wichtigen Fall, wo der ursprüng- 


lich ruhende Stab plötzlich an der Stelle x’ («’ > s;) etwa seines zweiten Teiles durch eine kon- 
stante Lokallast G getroffen wird. 


+09); (k=1,2,...,n—1) a (13), 


en 
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In üblicher Vorstellungsweise geht es hier um eine auf der kleinen Strecke |e — x’| < z 
gleichmäßig mit der Intensität - verteilte Belastung. Nach (5) ist dort 


G 
AI; 
2 


und überall anderswo ®,(x, p) = P;(t, p) = 0 zu setzen. 
Es ist daher nach (11.3) P,(p) = 0 und die Gleichungen (13) reduzieren sich zur einzigen 
Beziehung 


As [0 
. B, —— Ag! B 

2. C, = M,(p) | 0 er ee ER (16). 
D, | D, 


Die Matrix M,(p) hat die Elemente (15) und für B,, D, ergibt sich in unserem Spezialfalle aus (14) 
nach einfacher Rechnung 


ONE shi, (l— ®) 

hr 2E,1I 48 [sh SA ch}, (l—s)+hı ch SAL Shds (l — sı)] E; 
D G sin A, (I — x’) Ban 
1 


is EI "28 [sin s; A, cos A, (—s,) + hy cos s; A, sind, (— sı)] 

Für die Beiwerte A,, B,, C,, D, erhalten wir nach (16) und (15) leicht 
A, = (sh s, A, ch a —hıchs Ash sh) BI» Ba = (—shsı Ah shsıda+hych SsAgch3,.A) Bi: 
C, = (sin s; A, 608514, —h, c0S $,A, sin 5A) Di, D, = (sin sı 4, sin sı A, + hı cos sı4, C0S5,A,) D, 


womit die Berechnung aller Integrationskonstanten A,, B....., D, zu Ende geführt ist. 

2. Die bisherigen Ergebnisse ermöglichen uns nun geschlossene Ausdrücke für die Laplace- 
Transformierten Y,(z, p), Ys(, p) der gesuchten Lösungsfunktionen y,(%, 2), Y, (x, t) aufzuschrei- 
ben. Wir haben nämlich nach (8) 


2.) = 02r<s; Ze pel ene; 


(18) 


sh), (c— x) — sin, (@— a), U <r<I ER ak), 


G 
Zx(X, pP) = IE, 1% [ 
und die Formeln (8), (10.1), (17)—(19) ergeben sofort 
G sin A, x sin A, (l— ’) 
Y. == 1 2 
ı(@ P) 3 E, 138 Isin sı A, c08 A; (l— 51) + hy cos 5; A, sin A, (1 — sı) 
shA, 2shA,(1—) 
sh s, A, chi, (I— 5) + hu ch 51 A, sh), (L— 5ı) € 
G sin sı A, 008 A, (@ — sı) + A cos 5; A, Sin, @—Sı) d2el20), 
y 2 4 2 1 ı ı 41 2 1 Be 
„(@ P) 3 E, 138 sin s/ A, c0sA, (I — 51) + hy cos 5; A, sin A, ( — 5) lee 
sh s)A, cha, @ — Ss) + huchsi ), shi, (2 — Sı) n 
— —_ —— — sh A, (l— W(«,p)|; 
sh SA ch2, (I— 51) + hy ch 5; A, sh, (L— 51) het Wirp) 
Wap)=0, 1 <e<e, Wa,p) = sh, @e — 2) — sn, (ce —rT), U ST< I 
Die Formeln (20) geben uns die gesuchten Laplace-Transformierten, woraus nun die Lösung 


selbst zu finden ist. Natürlich ist dies in unserem Falle nicht leicht. 
3. Wir haben in gewöhnlicher Schreibweise 


ep) en 
2ri p 


Dee, 


44% I) = Te Ve RN 


wobei die län 
sind. 


gs bekannter Wege genommenen Integrale mittels des Residuenkalküls zu berechnen 


Die dazu nötigen singulären Punkte der Integranden sind zuerst der einfache Pol p = 0 
und dann die aus den Wurzeln ß, der Gleichung 


sin hy sı ß cos (I— 51) ß + A cos haßsnlt—s)B=0...-..- (22) 
sich ergebenden einfachen Pole 
De, r1aP Gl, 2,2.) ET 
Die p;, gehören dabei zu den Integranden mit hyperbolischen, die p_, mit trigonometrischen 
Funktionen. 


E: 
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Für y,(x, {) bekommen wir dann durch eine prinzipiell klare, jedoch in Einzelheiten recht 
mühsame Rechenarbeit den Ausdruck 
Ned = Tg |) EIER) + AR +2S Lt E- DS 
2 
„8 cos hy ß, 5, sin ß, (I— 51) Sin h, ß,x sin, (I— x) 
TAB tie, sin 2ß, le) —U— 5) sin 2hyß, 5, 
ß, sind dabei positive Wurzeln der Gleichung (22). 
Noch viel komplizierter ist die Rechnung zur Aufstellung der Formel für y,(, i). Als Ergeb- 
nis erhalten wir 


G ’ 
1.) = gr Rt st + + BR + 
2 
— Ks [8 +3@— Ss? + @—s) BR + @— SI} —12(l— x] + we) 
: . : ; (25). 
122 3" sin 2 h, ß, sı sin ß,(I— x) sin f,(I— x) » 
* — B» [hi sı sin 2 8, (l— 5) — (l— sı) sin 2 h, P, sı] 
ur) else, Re a | 
Im Einklang mit der elementaren Balkentheorie müssen sich natürlich die Ausdrücke (24) 


und (25) noch weiter umformen lassen. Eine hübsche Rechnung gibt wirklich folgende Endfor- 
meln für die Biegungsschwingungen unseres zweiteiligen Stabes: 


G 
nd = TE Re ) + 2M—Dad—sl 


x cos hy ß, sı Sin B, (l — sı) sin h, P,xsin ß, (I — X’) 
4 1 Ri 1 1Pv = v 
FAR ns sn29, (=) — A ans 


cos Pi ai], 0<r<s. (24); 


cos ß, a, t}; 


cos ß, q, l, <258 Ce: 


yo, I) = . 5 ja A) [IR @ lea —e—x?) +2 (k —1)s$(l—a)] 


Q sin 2 h, ß, sı sin P, (L— x) sin ß, (I — x’) 
Dr Fear "ur- Tara = 
+ BRD en sin 28, ds) —(— sin 25] 
den im Bereich 2 <x << I geltenden Ausdruck erhalten wir durch Austausch von x und x’ in 
(26.2). 
Zu Ende stellen wir noch die Bedeutung einzelner in (26.1) und (26.2) vorkommenden Hilfs- 
größen übersichtlich zusammen: a, und Ah, berechnen sich nach (2) und (10.2) in unserem Spezial- 


falle zu 
ET yE 
er ne 
d2 Va a 


s, ist die Länge des ersten Stabteiles, ! die des ganzen Stabes und f, sind die positiven Wurzeln der 
Gleichung (22). 

4. Der Übergang von (20) zu (24) und (25) scheint vielleicht dem Leser dieses Aufsatzes ein 
wenig schnell. Darum ist es nützlich hier zwei — sonst nur selten in der Literatur vorkommende — 
Ergebnisse aus der Theorie der Laplace-Transformation mitzuteilen. 


Es seien %, c, a, ß, y, ö, h von p unabhängige Werte, u = yi xp und 
H(u) =shyuchöu+hchyushöu, Hu) =shyucheu+hchyusheu, 


cos Pr q, ı, s<ı<a (%.2); 


T(u) = sinyucosöu-+hceosyusinöu, Tı(u) = sinyucoscu+ hcosyusincu) ' Gi, 
Setzen wir weiter 
sinxusinßu shaushßu Y(p) 
Y(p) = ——- — 2, = —oet 
T(u) H(u) 8) pe" 
T,(u) sin H l nie 
Z(p) >= we a ß u (4) Ss ıp u £ g(p) Re: ZPp) ept 
(u) H(u) pu° 
so ist die Summe R[/(p)] aller Residuen von f(p) gegeben durch 


RUDI Sg RR HB HI a +) 


S sin a u, sin ß u, cos y u, sin d u Mat 
328 a Re Fe ne 4 
2 op sn2dn, Senn) a 295 


v 


Pr 


I 


u 


6: 


5 
=“ 
y 
A 


Ra = ze o + +38) na + ar 


. — + AR+3ZHY+heBY HAI —RY+h)aY+hd} 
<, sinßwsin2yasinö— 0m mt 
4 v v Ur N 
4 Ze due) NSS 292). 
In (29.1) und (29.2) sind die «, positive Wurzeln der Gleichung 
sin yıweos da -Reosyusnou=ÜW..,. ....n abe (30). 


3. Schlußbemerkungen 


1. Die vorstehenden Ausführungen sind erstens gemeint als Beitrag zur Theorie «kompli- 
zierter Schwingungsprobleme und zweitens als Beispiel für Anwendung der Laplace-Transfor- 
mation zur Lösung nicht trivialer Fragen der mathematischen Physik und Technik. Trotz aller 
Kompliziertheit läßt sich jedoch die Rechnung für kleine Werte von n(n = 2,3) genau ganz 


_ allgemein vollführen. Die technische Bedeutung unserer Theorie erhellt schon aus dem in 2 


behandelten speziellen Falle, dessen Lösung der Verfasser in der Literatur nicht gefunden hat. 

2. Die physikalische Bedeutung der im Abschnitt 2 zugrunde gelegten Annahme 4} = & 
ersieht man leicht aus den Formeln (3.2.2) und (10.2): es geht um einen zweiteiligen Stab, dessen 
Materialkonstanten durch o, E, = 0, E, verknüpft sind. 

3. Die zeitunabhängigen Glieder von (26.1) und (26.2) bedeuten natürlich die statische 
Durchbiegung des Stabes durch die an der Stelle x’ angreifende Einzellast @. 

4. Selbstverständlich verwandeln sich für A, = 1 unsere allgemeinen Formeln (26.1) und 
(26.2) in das elementare Ergebnis [3]) aus der Theorie homogener Stäbe. 

5. Sehr schön gestaltet sich die Verwendung unserer Theorie zur Berechnung eines drei- 
teiligen Stabes unter speziellen Voraussetzungen y =, = 1. 
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Beiträge zum Quellsenkenverfahren für die Berechnung 
von Überschallströmungen 


Von Johann Münch in München!: ?) 

In der vorliegenden Arbeit werden die Quellsenkenverfahren eingehend untersucht und dahingehend ver- 
allgemeinert, daß sie auch auf Drehkörper mit Unstetigkeiten in den Ableitungen der Konturkurve anwend- 
bar sind. 

In the present paper the source-sink-methods are discussed in detail and generalized for the case that the 
contour of the body concerned, which is supposed to have rotational symmetry, has a non-continuous derivative. 


Dans l’article present les: „„Quellsenkenverfahren“‘ sont examines de pres et ils sont generalises de fagon 
qu’ils sont applicables aussi aux corps de rotation avec des discontinuites aux derivees de la courbe de contour. 


B upensaraemoi pa6oTe nccHenyeTcH MeToN CTOKOB U HCTOUHUKOB, KOTOPbIH Pacııpoc- 
TpaHsiercan Ha cıyuali O0TeKAHUA TEIL BPAIICHUA, KOHTYPHAA KPUBAH KOTOPBIX MOMKET HMETB 
paspbIBHyIo IIPOU3BOAHYIO. 

$ 1. Einleitung 
In der Gasdynamik verwendet man zur Berechnung von Überschallströmungen um 
spitze und schlanke Drehkörper die Quellsenkenverfahren. Das bekannteste ist das Verfahren 
von v.Kärmän, Moore und Tsien. Es wurde von Sauer auf die Umströmung von langsam 
pendelnden Drehkörpern ausgedehnt. Ein anderes Verfahren stammt von Probstein und 
Charyk, ein weiteres von Stetter. 

Wir stellen diese Verfahren unter einem einheitlichen und allgemeinen Gesichtspunkt 
dar, indem wir Potenzansätze (mit positiv-ganzzahligen Exponenten) für die Belegungs- 
1) Mathematisches Institut der Technischen Hochschule München, Arbeitsgruppe für Gasdynamik. 

2) Auszug aus der Dissertation des Verfassers, TH München 1956. (1. Ref.: Prof. Dr. R. Sauer u. Doz. 


Dr. F.L. Bauer, 2. Ref.: Prof. Dr. J. Lense.) 
4* 
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funktion, durch die sich Potential und Geschwindigkeitskomponenten quellenmäßig dar- 
stellen lassen, betrachten. Die Bestimmungsstücke der Strömung drücken sich dann durch 
Kugelfunktionen 2. Art aus. 

Die oben genannten Verfahren ergeben nur Strömungen mit „glatten“ Stromlinien, 
d.h. sie eignen sich vornehmlich für Drehkörper, deren Kontur keine Knicke aufweist. Wir 
zeigen, daß und wie man durch Hinzunahme von Potenzansätzen für die Belegungsfunktion 
mit halbzahligen Exponenten (die Kugelfunktionen spezialisieren sich dann zu vollständigen 
elliptischen Integralen) auch Unstetigkeiten in den Ableitungen der Konturkurve des 


Drehkörpers approximieren kann. 


$ 2. Potentialgleichung, Anfangs- und Randbedingung 
In diesem Paragraphen werden die Grundlagen und Voraussetzungen für unsere späteren 


Ausführungen kurz zusammengestellt. 

Wir betrachten die stationäre, achsensymmetrische Überschallströmung um einen 
schlanken, zugespitzten Drehkörper im Rahmen der linearisierten Theorie. In einem System 
von Zylinderkoordinaten x, r läßt sich die Strömung durch ein Potential ©, das der Gleichung 


(rO),— Hr D)), = EEE ee RE 


genügt, beschreiben. Dabei ist x = yM: —1=ctgf, wo M die Machsche Zahl und ß der 
Machsche Winkel der Anströmung mit der Geschwindigkeit V, ist. Es sei r=k(x) die 
Gleichung der Meridiankurve des Drehkörpers mit der x-Achse als Symmetrieachse und der 
Spitze xe=r=(0. Wegen ihrer Linearität gilt Gl. (l) sowohl für das Geschwindigkeitspotential 
Ö als auch für das Störpotential g = ® — y mit dem durch 9 = V,x gegebenen Potential 
der Grundströmung. Linearisierungsvoraussetzung ist, daß o, und 9, klein gegenüber V, 
sind. 

Als Anfangsbedingung gilt 

D’—= MHbzw. ol ar ur 0: 
Die Forderung, daß der Geschwindigkeitsvektor in die Tangente der Meridiankurve 


fällt, liefert unter Vernachlässigung der kleinen Größe 9, gegenüber V, die Randbedingung 


in der Form 
DREEV LEE) für r—ka)=0. 


$ 3. Spezielle Lösungen der Potentialgleichung 
Es werden durch einen Separationsansatz spezielle Lösungen der Potentialgleichung (1) 
gewonnen. Diese Lösungen erhält man auch aus der bekannten Integraldarstellung für das 
Potential bei einem Potenzansatz für die Belegungsfunktion. Ferner werden einige Relatio- 
nen zwischen den Bestimmungsstücken der Strömung angegeben. 
Nach Transformation auf die durch 
2 = oltgr, ar=— 
E =. x Eintr 
definierten Koordinaten o, r erhält man aus Gl. (1) für das Zusatzpotential die Gleichung 
(*Po)e = 9 Sin? T 
und hieraus folgt mittels des Separationsansatzes 
9,(0, 7) = 0’ y, (Üig r) 
die Legendresche Differentialgleichung 
Eye) dy,(d) 
ION AR v g 
mit ö&=(tgr. Wir beschränken uns dabei auf reelle ». 


Lösungen y,, für welche die zugehörigen i i Ü i 
2 „ p, die Anfangsbedingung erfüll 
die Legendreschen Funktionen 2. Art. Wir können also ER a ne 


P=9l0,T) = Ve Q,(Ciar). Tür 22 re (2) 
und ie 0 2 Raum der ungestörten Grundströmung. 
ie in der Einleitung erwähnt, lassen sich die dort 
a ‚la S genannten Quellsenkenverf. 
unter einem einheitlichen und allgemeinen Gesichtspunkt e wenn en 
bekannten Integraldarstellung a 


BET 


u nee 


| een 
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für das Potential die Belegungsfunktion F(£) auf der x-Achse in der Form 
Fo a ee Mir >00 
0 für 2=e0 
untersucht. Die zu f(£) = &* gehörigen o, sind identisch mit den (bei v» = n) durch GI. (2) 
gegebenen. 


Für die Geschwindigkeitskomponenten u=g9,undv = häl : 
Darstellungen p Pz v = g, erhält man aus Gl. (2) die 


aloe) V wor (EIN Den A ea era), 
DOSE ED Vo ot Le (Old nd oh) 
für den Raum £ >xr, in dem sie nicht verschwinden. 


Die bekannten Rekursions- und Differentiationsformeln für die Kugelfunktionen 
liefern mit » #0, —1,—2,... die Beziehungen 


G+Dd) sa = HEDolgr re pa . 222.) 
ya 
al U,+1ı = (2 N 1) 0 &ig a 0° a A (4b), 
Belyonı (29 1yonstg ev, 700, 0 2 rn (4c) 
und 
ou, Qu; 100, 
ar = yU,—-5» + == a en BT (52), 
oo, v—1 Dr 1 
rer D, z N a (5b). 


Gl. (4a) gilt auch für v» = 0, wenn formal (v 9,-1),=0 = Yo gesetzt wird. 


Aus den Darstellungen der Legendreschen Funktionen durch die hypergeometrische 
Reihe gewinnt man die Reihendarstellungen 


u Vo ye KR" F6,v5 27,09) WR VER Le 
= Aa gi En Ela ee) a 
mit X = nr k+k®=1. 
Bei bekanntem u, und v, ist 9, gegeben durch i 
= gan, +o@— Mu) Se 


$4. Die Sonderstellung des Falles » = 1 /2 
Es wird gezeigt, daß sich bei einem Potenzansatz für die Quellbelegung im Falle halb- 
zahliger Exponenten die Bestimmungsstücke der Strömung zu vollständigen elliptischen 
Integralen spezialisieren. Im Falle des Exponenten 1/2 weisen die Geschwindigkeitskompo- 
nenten auf dem stromabwärts gerichteten Machkegel einen endlichen Sprung auf. Die Ver- 
mutung dieses Sachverhaltes entstand beim Studium einer Arbeit von Heinz [1]?) und hat 


den Anlaß zu folgender Untersuchung gegeben. 
Stromaufwärts der Machlinie «-ar =0 verschwinden Potential und Geschwindigkeits- 


komponenten. Um das Verhalten dieser. Funktionen bei Annäherung von rechts an die 
Machlinie durch den Koordinatenursprung zu untersuchen, betrachten wir für} — + 0 den 
Grenzwert von g,, u, und v, auf einer Machliniiee—ar=44A=const> (0, beir=const + 0 
(zunächst für positive Indizes v). Unter Berücksichtigung von 

Be 
A+ar-+o 


o=fifiFlor, R- 


3) Die Zahlen in eckigen Klammern verweisen auf das Literaturverzeichnis am Schluß der vorliegenden 


Arbeit. 
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liefert die Untersuchung des Faktors 0”! k2” in den Reihendarstellungen (6a) und (6b) die 


Grenzwerte 
0 für re >1j2 


- p) v J \ 
nor ar lm Dr) Bl m Vz ur» u =L[2E RE 
ee 3 ar A+or+e i ar NEE r 


o&', Tür, „u 132 


Der Vergleich der Gln. (6a, b) mit Gl. (8) zeigt, daß die Geschwindigkeitskomponenten 
im Grenzwert A — + 0 gegen 0 oder oo streben, je nachdem, ob »v größer oder kleiner als 
1/2 ist; nur im Falle des Index 1/2 weisen sie auf dem stromabwärts gerichteten Machkegel 
einen endlichen Sprung auf. Die hypergeometrischen Reihen (6a, b) sind bei v = 1/2 identisch 
mit denen der vollständigen elliptischen Integrale und man erhält 

V,k aVok 
in = Kim e— —— Eh) Ne ae 
1/2 2yo ( ) 12 2Yo 1% 

Zusammenhänge zwischen Torusfunktionen, zu denen sich die Kugelfunktionen bei 
halbzahligen Indizes spezialisieren, und vollständigen elliptischen Integralen hat Airey [2] 
untersucht. Im folgenden benötigen wir nur noch den Sonderfall, daß 2» + 2 eine negative 
ganze Zahl ist, ohne daß » negativ-ganzzahlig ist. Eine bekannte Relation zwischen den 
Torusfunktionen liefert 


2 5° 
Un+32 — 5, 4 ja Dnysp > eu 
mit ganzzahligen n und man findet: 
i V,a(2n +2)! 
lim e$ se ar+4,r) = e END, 
ee ee In + DIE @& «a r)r+1Y2ar nn 
R — V,n.(2 1)-@2 
im ar tin = en Le 
a>+0 2 [(n +1)! a ryr+1Y2ar 
ö — Van Dep 
Im 02, = > 
u ve(ear+ Ar) 2 EN I KA EN 


Die Beziehungen für die zur Bestimmung des Druckkoeffizienten maßgebliche Größe u 
können unseren bisherigen Ausführungen bzw. der am Schluß dieses Paragraphen an- 
gegebenen Tabelle entnommen werden. Zur Erfüllung der Randbedingung ist für uns nur 
noch die Geschwindigkeitsquerkomponente von Interesse, deshalb wollen wir weitere Be- 
ziehungen nur für diese angeben. Zunächst findet man (bei positivem ») 


0 für » >12 

Ss in vo I ea 

2 Son form Eee ae 
6%) für »<1/2 


und mit Hilfe der Differentiationsformeln (5a, b) die Grenzwerte 


. OD, oVv 
lim Be lim Kerr n- R = 
i>--0 x a=ar-4A i>+0 or z=ar+i 0 fürn >= - 2 
lim = hei lim Oayıa, das RE x 
4>-+0 0% s=or+i X >40 or IP PSH A 9r - 
lim Be zer 3 lan | ee 3 Ver 
> +0 TC Ig=artA IXH>+6 or a=or-A 16 ry2« r 
und daraus 
lim ze] _ Yor[7a Kar) —3] 
i>+0 I z=ar-ti 16 ry2ar . (12a), 
lim dvzıs a. _ 83 Wrlak(ar)—1] V a ; 
A +0 dx | HIERHER 34 » (12h. 
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In folgender Tabelle sind Potential und Geschwindigkeitskom i inbli 
onenten im Hinblick 
auf die praktische Anwendung der Theorie kurz eingestellt, - Be 


1 


Y e”’9,= U, (Cigr) 


T 


k3 C(k) 


r-&tgr—1 


Sa) Ci) — Kih)] 


.- 
1, Br 1), Eigr 


1 1 

era rn ig) | m, 0 vo, = Duni (Üigr) 
— &jr 

En 
DR: 


E(k) 


— Gint 


k3 


1 [ER) Ch) — KM] 


1 
.2Cigr— 2 ER ge 
r.2Ctgr er Cr 


Dabei ist C(k) durch CC = (2 —R) K-—2E definiert. Mit Hilfe der Rekursions- 
formeln (4a, b, c) läßt sich die Tabelle beliebig weit fortsetzen. *) 


$ 5. Umströmung eines schlanken und spitzen Drehkörpers 


Zunächst werden die in der Einleitung erwähnten Quellsenkenverfahren kurz beschrie- 
ben. Weiterhin wird gezeigt, wie sich durch einen Potenzansatz für die Quellfunktion der 
Strömungsverlauf in der Umgebung eines beliebigen Konturpunktes dergestalt berechnen 
läßt, daß die Strömung nur von den Krümmungseigenschaften der Kontur in diesem Punkt 


abhängt. 
Zur praktischen Bestimmung der Belegungsfunktion setzen Probstein und Charyk [3] 


diese als Polynom n-ten Grades 
d=-23%%; a, = const 
»—1 


an. Zur Festlegung der Konstanten a, werden n verschiedene Punkte (X, r;) der Kontur- 
kurve r = k(x) herausgegriffen und die Randbedingung liefert mit Gl. (3b) für die a, das 


Gleichungssystem 


a, [ed (ig NM, = Ol, r 5 Er 


2 
IMs 


4) Nach Drucklegung der vorliegenden Ausführungen wurde ich auf die zwei Arbeiten ‚‚Practical caleulation 
of second-order supersonic flow past nonlifting bodies of revolution; NACA TN No. 2744“ und ‚‚First- and second- 
order theory of supersonic flow past bodies of revolution; Journ. Aeron. Sc. 18, 161—178 (1951) von M.D.van 
Dyke aufmerksam, in denen Potenzansätze für f(&) mit den Exponenten v = — 1/2, 1/2 und 3/2 zur Approxi- 
mation von Sprüngen in der Konturkurve (auf die Diskussion dieses Falles sind wir nicht eingegangen), deren 
Tangente und Krümmung verwendet werden. Die dort für Potential und Geschwindigkeitskomponenten mittels 
der vollständigen elliptischen Integrale E(t) und K(t) angegebenen Beziehungen gehen durch die Landen-Trans- 
formation t = (1 — k)/(1 + k’) bzw.k = 2t/(1 + t) in die entsprechenden bei uns über. 
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Nach v. Kärmän und Moore [4] wird der Kurvenverlauf der Belegungsfunktion durch 
einen Polygonzug approximiert. Mit der durch ’ 
| 1 für2>0 
I 0 fürz<o 
definierten Heavisideschen Sprungfunktion stellt sich die Quellbelegung analytisch in 
der Form 


N 
105) == = a) VAce &, d,; & —= const, & = 0 
v1 
dar. Zur Bestimmung der Konstanten a, werden wieder n verschiedene Punkte (z;, r,) der 


Konturkurve herausgegriffen, die Punkte £, auf der Abszisse durch Er = u, — ar, fest- 
gelegt und man erhält für die Konstanten a, das gestaffelte Gleichungssystem 


I VE TE] ran Won; = bion. 


Bild 15) 


Eine Verallgemeinerung dieses Verfahrens stammt von Stetter [5]. Der Kurvenver- 
lauf der Belegungsfunktion wird stückweise durch Parabeln k-ter Ordnung ersetzt und die 
Parabelstücke werden so gewählt, daß f(£) samt den k—1 ersten Ableitungen auch an den 
„Flickstellen‘ stetig ist. Die Konstanten b, in dem Ansatz 


k—1 n 
fd) = 2 Dee 2, (EEE, VYEe—8) 5b,c,&,=cont, 4=0 


werden nun folgendermaßen festgelegt: Die Konturkurve r = k(x) sei in einer Umgebung der 
Spitzerr = r== Oidureh (dievPotenzreihe 5 y,%,y, = const, gegeben. Es wird nun mit der 
Reihe v(«) BEN v,|z, k(x)] ein Koeffinientenresteielt so durchgeführt, daß alle b, durch 
die k—1 sten Ahleitinzen von k(x) an der Stelle x = 0 bestimmt sind. Man findet 

En e ee 


b, = ——— 2 5 .. 
1 —o2y? Yı 140°, Mol, 


Die n Konstanten c, kann man wie bei dem Verfahren nach v. Kärmän und Moore 
festlegen. n verschiedene Punkte (2, r,) der Konturkurve liefern mit &11 = 2; — «ar, wie 
dort für die c, das gestaffelte Gleichungssystem 


n DER ji k—1 
ZA En Yin) = | SA, ig]; 
v= v=1 % N, 
a 
Den vorausgehenden Untersuchungen nach handelt es sich also bei den Quellsenken- 
verfahren im Grunde genommen um eine Entwicklung des Potentials einer achsensymmetri- 


schen Überschallströmung nach Kugelfunktionen 2. Art, d.h. nach einem die Anfangsbe- 
dingung erfüllenden System von Lösungen der Grundgleichung. 


Bei dem Verfahren nach Stetter hängt der Strömungsverlauf um die Spitze nur von 
deren Krümmungseigenschaften ab. Es kann also bei analytischer Vorgabe der Kontur- 


5) In Abb. 1 soll es richtig x,, r, statt x,, r, und ar; statt «&, r, heißen. 
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Kürver— Kr) = = BER CUNST: I einer Umgebung der Spitze = r = 0 der Strö- 
mungsverlauf durch ne Ansatz /(£) = 5 b,& für die Belegungsfunktion, wobei die Kon- 


stanten b, wie oben angegeben zu ne: sind, theoretisch in dieser Umgebung hinreichend 
genau berechnet werden. Das ist quasi ein Analogon zur Approximation einer analytischen 
Funktion durch ihre Taylorreihe. 


Die Frage ergibt sich, ob eine ähnliche Entwicklung nicht auch in der Umgebung eines 
vom Nullpunkt verschiedenen Punktes (x*, r*) möglich ist und zwar dergestalt, daß sich bei 
analytischer Vorgabe der Konturkurve in dieser Umgebung der weitere Strömungsverlauf 
berechnen läßt, indem man die Belegungsfunktion in einer Reihe nach Potenzen von &-&*, 
&* = x* — ar*, ansetzt und die Konstanten dieser Reihe direkt aus den Ableitungen der 
Konturkurve an der Stelle (2*, r*) bestimmt. Potenzen für die Belegungsfunktion mit positiv- 
ganzzahligen Exponenten n sind dabei sicher nicht brauchbar, da, wie gezeigt wurde, die zu- 
gehörigen v, an der Machliniex—ar= &* verschwinden und ihre ersten Ableitungen an dieser 
Machlinie 0 oder oo werden, je nachdem ob n>1 oder n=1 ist; entsprechendes gilt für 
die höheren Ableitungen. 

Die vorangehenden Untersuchungen zeigen aber, daß sich der Strömungsverlauf von 
dem Punkt (x*, r*) ab stromabwärts (ohne Rücksicht auf die ankommende Strömung) wie 
angegeben berechnen läßt, wenn man die Belegungsfunktion in der Form 


ME Er = Da, (Erin Ye Er),  d,= const 
n—1 


ansetzt. Die Geschwindigkeitsquerkomponente schreiben wir bei dieser Quellbelegung in 
der Form 


riet 24 mın@—&*,r) für 222% 
0 aurgerr au 
Die Randbedingung liefert nun mit Gl. (11) sofort 
v(x*, 7*) = d, v112(a 7%, r*) = Vo [Kor vr - 
Ebenso folgt wegen 


dv, [x, k@)] _ K'(@) % nn k 
re a)? I, K(a)] +1 —& a [2 k(&)] 
aus der verallgemeinerten Randbedingung 
dm v[x, k(x dm +1 k(x 
| a Bar a 0 Br 


mit m = 1 sofort 
d 2 Ir 
ann], „= Volk’Oken- 


Ist die Konturkurve in einer Umgebung des Punktes («*, r*) durch die Potenzreihe 


Ka) 2; (x — x*)”, 6, = const, gegeben, dann ergibt sich die Konstante d, nach Gl. (11) 
v1 


in der Form 


„_ _2ayar 
7 & 
und die Gln. (12a, b) liefern 
1. Vor@e&—d RE Vor (a2, u a 
16 r*Y2 «a r* 


also 
16 6, — ö, (7 & 6, 3) 


ee: nr 


2 


rare 
ie weiteren Konstante „ erh: 
 bedingung. x: fr 
MorkEsiwärd er für die ee er Sinn habe 
mittels solcher Potenzreihenansätze für die Belegungs: 
b„ bzw. d, zu berechnen. Ganz abgesehen vom hohen Recl 
ob genügend hohe Ableitungen der Konturkurve überhaupt praktisch 
können und zweitens ist es fraglich, ob ein solcher Aufwand im | Rahme 
Theorie physikalisch noch sinnvoll ist. Es wird deshalb genügen, wenn a 
Tangentenrichtung nur noch die Krümmung an geeigneten Stellen der Konturkurve 
sichtigt. Im Sonderfall der Approximation an der Spitze allein (Stetter) heißt da, - 
die N stückweise aus Parabeln dritter Ordnung aufbaut. ek; 


86. Vorschlag eines neuen Quellsenkenverfahrens 


Auf Grund der Untersuchungen in $5 wird das von Stetter angegebene OS He 2. 
verfahren dahingehend verallgemeinert, daß der Strömungsverlauf nicht nur an der Spitze, 
sondern auch an den anderen (zur Bestimmung der Konstanten in der Belegungsfunktion 
herausgegriffenen) Konturpunkten durch die Krümmungseigenschaften der Konturkurve B 
in diesen Punkten festgelegt wird. 

Bei unseren Betrachtungen im vorangehenden Paragraphen haben wir zunächst v(z, r) = 
für x < x* vorausgesetzt. Das ist im allgemeinen aber nicht der Fall, d. h. die Konstanten d, 
lassen sich dann nicht mehr in der oben angegebenen einfachen Weise berechnen. Wir nehmen 
an, daß die Strömung in dem Punkt (x*, r*) bereits festgelegt ist und dort die Randbedingung 
der tangentialen Umströmung erfüllt, was durch irgendeines der drei beschriebenen Quellsenken- 
verfahren erreicht worden sein mag. Die Konstante d,, die durch die Tangentenrichtung der 
Konturkurve in (x*, r*) festgelegt wurde, wird dann null. Zur Bestimmung der Konstanten d, 
durch die Krümmung in (x*, r*) ist die Ableitung von v[x, k(x)] an dieser Stelle zu berechnen. 
Es darf dann jedoch kein Term a*v, (e — &*, r), a* = const, in v (x, r) vorkommen, da dessen 
Ableitung bei a, an die Machlinie —ar = £&*, wie früher gezeigt wurde, gegen 
oo strebt. 

Wir werden nun ein neues Verfahren zur praktischen Bestimmung der Quellfunktion 
angeben. Wie bei dem Verfahren nach Stetter soll der Strömungsverlauf um die Spitze des 
Drehkörpers nur von deren Gestalt abhängen, wobei jedoch außer der Tangentenrichtung nur 
noch die Krümmung berücksichtigt werden soll. Außerdem soll aber auch die Stromlinie längs 
der Kontur in den zur Festlegung der übrigen Konstanten herausgegriffenen Punkten (,, r;) 
gleiche Krümmung wie diese besitzen. Demgemäß setzen wir die Belegungsfunktion um die 
Spitze als Parabel dritter Ordnung an, wählen aus der Konturkurve n verschiedene Punkte 
(x, r;) in hinreichend enger, sonst aber beliebiger Aufeinanderfolge aus und bestimmen damit 
die Punkte & = 0, &41 = %, — ar, auf der Abszisse. Zwei Konstanten der genannten Parabel 
werden durch Tangentenrichtung und Krümmung an der Spitze, die dritte Konstante durch 
die Tangentenrichtung im Punkt (x,, r,) festgelegt. Nun nehmen wir zur Quellbelegung in jedem 
der Punkte Z,,, zwei weitere Potenzen mit den Exponenten 3 und 3/2 hinzu und legen den ersten 
Term jeweils durch die Tangentenrichtung in den Punkten (&;11, r;+1) und den zweiten Term 
durch die Krümmung in den Punkten (x, r,) der Konturkurve fest. Analytisch stellt sich also 
die Belegungsfunktion /(£) in der Form 


Ebern nr oelsg on 
fe) = bE+b2tca®+olE— 8) + d, (E— £8,)°2 für B<E<SE 


DE+bERH 2% e—EP4 Sa., le ae 


dar. Die Konstanten b, und b, können wir dabei als bekannt voraussetzen, sie sind identisch 
mit den bei der Beschreibung‘ des Quellsenkenverfahrens nach Stetter angegebenen. 
Setzen wir zur Abkürzung 


5 ki) 


d: am 


: 
; 
ü 
. 
E 
j 


BARS: für ml 
Zur IR”«,) Für mern 


und 


— [2 — &, k@)] 
dam 


a I (&, — 8) für m=0 
Zw tu v; (X, — &,) für m=1 
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dann lautet das gestaffelte Gleichungssystem für die Konstanten c, und d, folgendermaßen: 


2 
e b, v,(%) + C Val) = Vok'(&ı) ; 


2 
3 b, D,(&) + 1 d3(&) + dy V3j2(a 7) = Vok''(&); 
2 2 S 
P> b, D,(X) 1 2% Dz; (X =Iz =) 3 d, LEIE (X; a &,) — Vo k'(&,) 9 


Me» 


2 
b, v,(&;) 2 C,Dd5 (&, — £&,) + di va (& — &) + d; v3j2 (& 75) = Vo Kk’(&;) 


1 


x 
|) 


Blaue ae iente sta ale al we hu.te else! oje, Kohrei area.” si Te Yarile. Be es see 


2 n—1l n—1 
> b, A) an 2 C, v3; (1—&)+ 2 ur v3[: (a1 — 8) + dn-ı D38 (I) FVo Kl) 


1 v 


Dabei ist v3j2(& r,) nach Gl. (12b) gegeben durch 
Aral 
4 De 


v3j2(&r,)= 
Die v; (x, — £,) sind mittels der Relation 


A|) 


T, 


KT, 
( ] D)—_ı (X, A &,) 
T, 
zu berechnen. Nach Bestimmung der Konstanten durch Auflösung des Gleichungssystems hat 
man die Darstellung 


2 N N 
v2, r) = 2d,v,@ Tyer 2% DASS, DW.) a Yan (&— Eur) Y@—ar—5,) 


in welcher an Stelle von y nun 9, U, eingesetzt werden kann. 


$ 7. Drehkörper mit Knieken in der Konturkurve 


Die Ausführungen des $ 4 haben ergeben, daß die Geschwindigkeitskomponenten im Falle 
des Index 1/2 auf dem stromabwärts gerichteten Machkegel einen endlichen Sprung aufweisen. 
Wir verwenden diese Tatsache dazu, um Unstetigkeiten in den Ableitungen der Konturkurve zu 
approximieren. 

Hat die Konturkurve an der Stelle (x, r,) einen Knick, so kann man diesen folgender- 
maßen approximieren: Zur Bestimmung der Strömung sei irgendeines der beschriebenen Quell- 
senkenverfahren ausgewählt und die sich ergebende Darstellung der Geschwindigkeitsquer- 


Bild 2 


komponente sei (bei noch nicht festgelegten Konstanten) mit d(z, r) bezeichnet. Den Potenz- 
ansatz für die Quellbelegung ergänzen wir nun durch Hinzunahme des Terms d, (£ — £,)!!? und 
aus 
v(% +06, 7) =0(, +07) + ddus (BETSR,) 
und 
tn )=- en —0,r) 
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folgt wegen D (x, + 0, r,) =D (2, —0, r,) sofort 
dvds lar, tr) =ple, Or) ve, 09). 
Zur Festlegung der Konstanten d, kann man nun 
vs, +09, r)= WE, +0) 

und damit er 
dvn rs, r)=V Ka +) —0,;r,) 
als weitere Bedingungsgleichung zu denen für die anderen Konstanten hinzunehmen. 

Eine andere Möglichkeit ist, durch 

v (x, +0,17) —v(,— 0,7) = Volk (, +0) — K (a, — 0)] 
die richtige Sprunghöhe vorzuschreiben. Die Konstante d, ist dann nach Gl. (11) unabhängig 
von den anderen Konstanten direkt gegeben durch { 
ddr (oT, 1) = Vol +) —- KO], 


= 2 Ve OR +0) 
IT & 


Dabei ist zu beachten, daß v(x, r) bei Verwendung dieser Formel an der Stelle (x,, r,) zwar 

den Sprung von der richtigen Höhe hat, daß jedoch dann 

v(+0,17)= Vok (+0) 
ist, falls 

v(% —0,r,)= Vok (x, — 0) 
ist. Es wird deshalb zweckmäßig sein, die anderen Konstanten so festzulegen (oder zu dem 
Quellbelegungsansatz einen weiteren Term hinzuzunehmen und die zugehörige Konstante so zu 
bestimmen), daß v (x, —0, r,) = V, k (x, — 0) wird. 

Auf Grund der Ausführung in den vorangehenden Paragraphen liegt nun der Gedanke 
nahe, auch noch die Krümmung der Konturkurve an der Stelle (x, + 0, r,) zu berücksichtigen. 
Kommt in v(z, r) kein Term mit v, (@ — £,r) vor, also in /(&) kein Glied mit (£ — £,), dann 
fügen wir zu dem Potenzansatz für die Belegungsfunktion noch den Term d, (£ — £,)®'2 hinzu. 
Aus 

v (X, +0, 1A) = v (2, +0, T,) == d, vij2 (& U) T,) fr d, D372 (& Ts T,) 


erhält man mit der Forderung 
v’ (+0, 1r)=V,k'@, +0) 
für die Konstante d, die Gleichung 
d, van & ru) = Vok’ +) —V (+0,17) —dı die (& Tr, r,) 


als weitere Bedingungsgleichung zu denen für die anderen Konstanten. Dabei sind vi, (ar, r,) 
und v3) (& Tr, r,) gemäß den Gln. (12a, b) gegeben durch 


i V,[7ak 0) —3 
en Ne 
F,VLaAo&T, 
3aV lkw +) —1ı/a 

4 2% 


v3/2 (& Tu T)= “3 (14bR: 

Wir wollen nun zur Veranschaulichung dieser Ausführungen die Approximation von 
Knicken in der Konturkurve in das in $6 angegebene Quellsenkenverfahren einbauen. An der 
Stelle (x,,r,) der Konturkurve befinde sich ein Knick. Wir setzen die Belegungsfunktion in 
der Form 


IE), Ib, nr 26 Ga lerne) 2 du, (EEE YE—E) + E— Er) VE YlEE,.1), 


wobei die Konstante d durch Gl. (13) gegeben ist, an und erhalten für die Konstanten c, und 
d, das Gleichungssystem 


2 u—l u—1 
SbvlE)+ 3,0; —EJ)YeÜ-W)+3 bad y —E)Ylii-W)=V, Kl): 
v1 v=2 


v=1 


ne, 


ie 1 ee 


un 
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2 ’ ei Br ” & "£ 4» 
aD a Pr a ei, 
el, 2 eu le 


2 7 t 
BL v,(2,)+ 2 €, Dz (X, 5) m & dd Dajg (&, Dar &,) == Vo k’ (& = 0) ’ 
Ei v,(&,)+ 2 c,d3 (X, — &,) are > ee v3 (u — &,) + d- viplar,) = Vok’ (u +0), 


> b,v,(&) + x Bla, E)YÜ— + 2 d,_ı Days (,; — &,)Y (e— 9) 


2 d: Da — Eur) = Vo Eko i=u+la+2, 
zb, v,(x;) + 3:4 v3 (©, — &)Y(i —v) + Zi v3,2 (X; N +1—») 
+d: ee u+lu+2,....n—1. 


$ 8. Illustration der Theorie an Hand eines 


Wir betrachten einen Drehkörper, dessen Meridiankurve mit einer Geraden beginnt. An 
diese ist eine Parabel so angesetzt, daß an der Ansatzstelle ein Knick vorhanden ist. Die Parabel 
wird weiterhin mit stetiger Tangente in eine Parallele zur Symmetrieachse fortgesetzt. Gemäß 
der Randbedingung v[x, k(x)] = V, k'(x) vergleichen wir die Ableitung k’(x) der Konturkurve 


mit der Funktion ee [%, k@)]- 
0 


Alae) 


Bild 3 


Es zeigt sich deutlich, daß die bekannten Verfahren zur praktischen Bestimmung der 
Quellfunktion durch Hinzunahme der Terme mit halbzahligen Exponenten zu /(&) wesentlich 
verbessert werden. Um den Zusammenhang zwischen v und f(£) besser zu veranschaulichen, 
ist jeweils der Verlauf der Quellfunktion gestrichelt mit eingezeichnet. 

Die Knickstelle befindet sich im Punkt (x, r,). Für die auftretenden Konstanten sind 


folgende Bedingungsgleichungen vorgeschrieben: 


b: [v(z, k(x))]x=0 = 1% [k’@)]z=0 =; Vo tg 0; 
d,: BE k))] z=2, = V, K(o)k=,—0 „12,3, 
c,: [ol k@))h=, = Volk, = 0 a 
di: [v(x, Ka))lo=a,r0 = Vo K Ok-n+0; 
d,: [v/(@, K&)k=-2+0 = Vo IK O)l=a.+0 a 


1. Verfahren nach Probstein und Charyk. 


a) = ImE+aE&] YO; 


en 
re 


je u 
7 N k PAR 
& s 


BE 
u ‘@; Verfahren nach v. Kärmän und M 


a) j&) = = ns. YO+ aueh) ve 
Zvoe) 
FE) 


Bild 7 


2) IO-BEYO + E-E)ALE-ENIYE-E). 


3. Verfahren nach Stetter. kai < ; 
a) IO=bEYO +E-HYE-H), 2 


Bild 9 


b) (9 =bEY(H) + I, (E— Er +d(E x: A YE-E: 


Bild 10 


YO+BE- HE ER + REM YE—E. 


Bild 11 


$ 9. Schlußbemerkung 


. Die vorliegenden Ausführungen lassen sich ohne wesentliche Schwierigkeiten verallge- 
meinern etwa auf den von Sauer [6] behandelten Fall des langsam und mit geringer Amplitude 
(räumlich) schwingenden Drehkörpers. Es ist lediglich eine sinngemäße Übertragung unserer 
Untersuchungen auf die Verhältnisse hier erforderlich. Als Spezialfälle sind in der Arbeit von 


Sauer enthalten a) der von Heinz [7] (allerdings mit einem anderen Ziel) untersuchte Fall des 


langsam und mit geringer Amplitude in einer Ebene pendelnden Drehkörpers, b) der von Tsien 
[8] betrachtete Fall des schwach angestellten und c) der in dieser Arbeit behandelte Fall des 
axial angeströmten Drehkörpers. 
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Reflexion ebener Druckwellen an einer festen Wand 
Von C. Heinz in Lörrach!) 


Es wird ein einfaches Verfahren angegeben, um die Reflexion einer aus einem Verdichtungsstoß und einer 
nachfolgenden Verdünnungswelle bestehenden ebenen Druckwelle an einer festen Wand zu berechnen. Hierbei 
werden die von Hadamard angegebenen Lösungen der Darbouzxschen Gleichung, die bei dem hier behandelten 
Problem als Potentialgleichung in der Geschwindigkeitsebene auftritt, in iterativer Form zur gleichzeitigen 
Bestimmung des Stoßes und des Strömungsfeldes zwischen dem Stoß und der Wand benutzt. Auch die inverse 
Aufgabe, bei der der Druck an der Wand gegeben und daraus die ankommende Welle zu berechnen ist, wird be- 
handelt. Als Beispiel wird die Reflexion einiger Wellen durchgerechnet. Ein Vergleich mit dem Ergebnis der 
wesentlich umständlicheren Charakteristikentheorie ergibt keinen merklichen Unterschied. 


The reflection at a solid wall of a plane pressure wave consisting of a shock wave and a subsequent expamsion 
wave is treated: by a simple method. The two-dimensional velocity potential of the problem satisfies the Darboux 
equation, and the solutions of this equation given by Hadamard are used in reiterated form to determine both 
the shock wave and the velocity field between wave and wall. T'he inverse problem where the pressure on the wall 
is given and the incident wave has to be found is also discussed. T’he reflexion of several waves is compu- 
ted as example. A comparison with results obtained by using the far more complicated theory of characteristics 
shows that there is no noticeable difference. 

Un simple proced& est indique afin de calculer la reflexion d’une onde de pression plane consistant d’un 
choc de compression sui vie d’une onde de depression-sur um paroi solide. Dans ce cas les solutions, indi- 
quees par Hadamard, de l’equation de Darboux qui se produit dans le probleme present comme equation 
potentielle dans le plan de vitesse sont appliquees sous une forme iterative afin de determiner en meme temps 
le choc et les lignes de courant entre le choc et le paroi. De meme le probleme inverse est traite, c’est ü dire la 
pression sur le paroi est donnee et l’onde arrivante en est ü determiner. Comme exemple la reflexion de quelques 
ondes est calculee. Une comparaison avec le resultat de la theorie des characteristiques qui est beaucoup plus 
circonstanciee ne prouve pas de difference sensible. 


1) Mitarbeiter am Service t6chnique de l’Arme&e Frangaise. 


Einleitung : 
Eine in ruhende Luft nach rechts, d. h. in x-Richtung, laufende 
stoß und einer nachfolgenden Verdünnungswelle bestehende ebene Drucl 
feste Wand. Um die Vorstellung zu präzisieren, nehmen wir an, daß der ze 
an der Stelle x = 0 der Wand, wenn die Welle nicht reflektiert würde, durch 


il; ee | nip 
Paul je * Dee Er N ; | x ET 
| Dm 7 to | Ex 
gegeben ist. Hierbei bedeuten allgemein p den Druck und p =p—p, den Überdruck über Er © 


den Ruhedruck p,, der in der ruhenden Luft herrscht. Für t = 0 ist p = p,„, nämlich der Maxi- 


Stoß, und für!=t, ist f =(, 
d.h. p=p,- Die reflektierte Welle 


die Strömungsgeschwindigkeit u an 
der Wand, d.h. für c=(, ver- 
schwinden soll. Dadurch liegt der 
reflektierte Stoß zunächst für x =0, 
t=0 fest. Die nachfolgende Ver- 
dünnungswelle durchkreuzt den 
Stoß und lenkt ihn dabei ab (Bild1). 
Ferner wird der Stoß aber auch 
von der an der Wand reflektierten 
Verdünnungswelle, die, da die 
Schallgeschwindigkeit hinter dem 
Stoß größer als die Strömungs- 
geschwindigkeit ist, den Stoß ein- 
holt, beeinflußt. 

Das Strömungsfeld hinter dem 
reflektierten Stoß kann in guter 
Näherung als isentropisch ange- 
nommen werden. Allenfalls läßt 
sich die Entropieänderung iterativ 
berücksichtigen. Während nun in 
früheren Arbeiten [1], [2] die Strö- 

Eee mung nach der Charakteristiken- 

theorie berechnet wurde, wobei der 

Stoß in das Charakteristikennetz eingegabelt werden muß, werden wir die von Hadamard 

für Gase mit 2n +1 Freiheitsgraden (für n —=1 oder 2 ist n gerade die Anzahl der Atome), d.h. 

y=2n + 3/2n + 1, angegebenen Lösung der Potentialgleichung benutzen und den Stoß gleich- 

zeitig iterativ mitbestimmen. Die Hadamardschen Lösungen sind im übrigen in letzter Zeit 

einerseits auf beliebige Werte von y [3], [4], [5], andererseits auf nicht-isentropische Strömungen 
[6], [7] ausgedehnt worden. 


1. Grundgleichungen 
Die ebene nichtstationäre Strömung eines Gases wird durch die Gleichungen 


(( +p),-+ (eu, =0 (Erhaltung des Impulses), 
(u),+, =0 (Erhaltung der Masse) , AEEERLT), 
us,+Ss, =() (Erhaltung der Entropie) 


beschrieben. Hierin bedeuten o die Dichte, p den Druck, s die durch die Differenzen Gi ardei 
spezifischen Wärmen bei konstantem Druck und Volumen dividierte Entropie, u die Strömungs- 


ergibt sich aus der Bedingung, daß 


ar 
ns 


 maldruck unmittelbar hinter dm 


ee 
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geschwindigkeit und x und £ die Raumkoordinate und die Zeit. Nach Einführung der charak- 
teristischen Variablen A, u und der Schallgeschwindigkeit a: 


re R; eure ran 1 a 

ud Dep Lea,, A ER 12) 
5 % 1 : 2 RR 
m) — s Ta I Ta Pa P(s, 0) 
2, u Ir 5 Arm, ET 


(a, ist irgendeine Bezugsgeschwindigkeit, z. B. die Schallgeschwindigkeit der ungestörten Luft) 
gehen die Gln. (1.1) in 


dx 


di = 3, ; E, ds längs der Charakteristik Ct: Fre +a, 
Ina dx 
du= — '—ds A Sr) Fr, "ii =U— rer. 
N u kr ee Ge 
dx 
dse=NV ER) „ „ er 
C FT, u 


> [8], (S. 140). Besonders einfach werden die Gleichungen für die isentropische Strömung, 
ei der 
UA=O längs C*, 
du — a) ” Cr D 
s = (0 im ganzen Strömungsfeld 
ist. 
Die Gleichungen des Verdichtungsstoßes lauten, indem wir die Größen vor dem Stoß 


durch p, a, usw. und nach dem Stoß durch p, a, usw. sowie die Geschwindigkeit des Stoßes 
relativ zur Luft mit c bezeichnen [8], (S. 247) 


2 De er 

DE rn > 

zn ye y—1\ y—1 RL 

ne etz) er, rn 
2 il, Br P(o) 
ee AB 


Das obere bzw. untere Zeichen gilt dabei für einen nach links bzw. nach rechts laufenden Stoß. 


3. Die ankommende Welle 


Die ankommende Welle besteht aus einem Verdichtungsstoß und einer nachlaufenden 
Verdünnungswelle, die den Stoß einerseits einholt und dabei krümmt, andererseits aber auch 
am Stoß reflektiert wird. Insofern ist die ankommende Welle nicht einfach. Die Stoßgeschwin- 
digkeit ist eine Funktion w(A,) des charakteristischen Parameters A,, indem nämlich im Punkte‘ P 
der Stoßfront mit dem Wert » die C+ mit dem Wert A, den Stoß trifft (der Index v soll andeuten, 
daß es sich um die Strömung vor dem reflektierten Stoß handelt). Dabei ist die Entropie- 
änderung in der Strömung nach dem ankommenden Stoß vernachlässigt, was dann berechtigt 
ist, wenn & sich nicht stark ändert. Es ist nach (1.4) 


ER 1 a Ir 1 
=5 en 5 U(o) St A(o). 


ae 


Andererseits ist 
1 1 
Mm =— 5 00) + Al). 


dh 
do 
sehr klein bleibt. Wir wollen also zunächst in der ankommenden Welle u, = u, = const setzen, 
d.h. sie als einfache Welle betrachten, um so mehr, als der tatsächliche Verlauf von u, ja erst 
durch die speziellen Anfangsbedingungen der ankommenden Welle festgelegt ist. Für u, + const 
vergleiche Ziffer 7. 


Hieraus folgt 


0, und man stellt fest, daß du,/dw, solange w nicht allzu groß wird, 


Qu 
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Ist o, die Geschwindigkeit des ankommenden Stoßes beim Auftreffen auf die Wand, so 
ist nach (1.4) unmittelbar hinter dem Stoß 


pe or Ä ee U 
Dr a P(o,) B) A, A(o,) D Ay 7 (©) 9 


4 =; Uw,) + TS] Aw), M=75 Uw,) + Pe "A(00) 


Die Werte von u, a und p sind in der ankommenden Welle als Funktionen von A, durch (1.2) 
gegeben, also = 


u a Fe] 
Fi u na htm. 
i 2Y a CO PR N re (2.2) 
Pe _ Be. Bi _ (BY. Pia) (2) "gan 
Po Pı Po a, 
Für die C+ gilt: 
dx 
re 
z x t 
woraus mt&= —-, T=- 
Ayo to, 
Be ER BEN A a Wale Sa 2.3 
E = ao (7 9(A,)) ’ a, — p) A, 5) DE Naeee ( . ) 


folgt. Dabei ist i, irgendeine Bezugszeit, etwa die Zeit des Nulldurchganges des Überdruckes 
in der ankommenden Welle. g(A,) ist der Wert von 7 für & = 0, der sich einstellen würde, wenn 
die Wand nicht da wäre. Ist in der ankommenden Welle für &=0 


B=/@. 1O=1 


gegeben, so wird nach (2.2) 


2 


! 1 + 2). Po) —1 FR: 


wo) 


Hierdurch ist die Funktion g(A,) bestimmt. Der Wert von &, bestimmt sich aus p, durch 


1 2 
a 


Tr ah RT ER 
0 


Pı ist der Maximaldruck p,, in der ankommenden Welle. Durch Vorgabe von f(r) und Pm _ Pı 
: : p pP 
ist also die ankommende Welle vollständig beschrieben: Aus (2.5) berechnen wir @,, aus (2.1) ih 
(2.4) liefert die Funktion g(A,),-(2.2) die Zustandsgrößen als Funktionen von A, und (2.3) die 
zugehörigen Werte von & und r, die in der einfachen Welle natürlich durch A, (und u, = ı) 
nicht eindeutig gegeben sind. 


3. Der reflektierte Stoß 


Vor dem reflektierten Stoß gelten die Beziehungen aus Ziffer 2 für die ankommende Welle. 
Um den Stoß festzulegen, brauchen wir irgendeine Zustandsvariable nach dem Stoß oder die 
Stoßgeschwindigkeit als Funktion des charakteristischen Parameters A, der ankommenden Welle, 
Die Zustandswerte unmittelbar nach dem reflektierten Stoß bezeichnen wir durch A (A,)s Pn(A,) 
usw. Es gilt für die Stoßlinie, die durch &,(A,) und 7,(A,) beschrieben wird: IE 


dr; var TO We ne Lan RE (3.1). 


u, und a, sind als Funktionen von A, durch (2.2) gegeben. Mit Hilfe von (2.3) folgt daraus: 


dt +1 d [u 
1 1 H, 8 Y or 
ON ee a 


ga). Be 
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Die Anfangsbedingung für diese Gleichung verlangt }, =}, für = 0, wobei A, durch (2.1) 
gegeben ist. Die Lösung von (3.2) schreiben wir in der Form 


%, (ot % ;) 
ae +1 AR 
a ny=--%7 re % ee EIN 
4, eo 2a 4 


Hiernach könnten wir r,, und damit nach (2.3) &,, berechnen, wenn ® = w(A,) bekannt wäre. 
Für ,=0,&,=0 ist o(},) = w, dadurch bestimmt, daß an der Wand, d.h. fürr=0,8=0 
unmittelbar nach dem reflektierten Stoß, u = 0 sein muß. Das liefert nach (1.4): 


Wu ___ I ab 
a, a, % 4a A(o,) : 


Die Lösung dieser quadratischen Gleichung ist bei Berücksichtigung von w, > 0: 
r 290 — (y— — en 


Po)+1-o Y@W-Dai+2’ Y—-DP@)+Y +1 


Den weiteren Verlauf von w(A,) werden wir später in Ziffer 5 iterativ zusammen mit der Strö- 
mung zwischen dem reflektierten Stoß und der Wand ermitteln. Aus w(A,) folgt 


2y 2y 
nl). P ee (de + my 
Po Pı Pı Po a (00) Fo) Sa, At ; 
ia Zur 0 & Fine 
rd am 49) 2 Arm) 
a Urt +r—m f 


Ay A, do Ay 


= Po) 0, + m,  Flo)=— ZU) +5 AD), 


ee) 


m— = Po) +)» Frlo)= +7 Ua) +5 AD) 


für die Zustandswerte unmittelbar nach dem Stoß. Die letzte dieser Gleichungen werden wir 
benutzen, um aus w,(A,) die Funktion w(A,) zu berechnen. F,(w) ist in dem uns interessierenden 
Bereich fast geradlinig, so daß sehr bequem zu interpolieren ist. Vor allem sieht man aus der 
vierten Gleichung (3.5), daß sich A, bei einer Änderung von u, bzw. w nur wenig ändert; denn 
die Entwicklung von F, nach Potenzen von — 1 beginnt erst bei (o — 1)? (F, drückt die 
Entropieänderung beim Durchgang durch den Stoß aus). Kennen wir also angenähert o(A,); 
so wird A,(A,) schon eine ziemlich gute Näherung sein. 


4. Das Strömungsfeld hinter dem Stoß 


Ein eindimensionales isentropisches Strömungsfeld läßt sich durch eine Lösung der 
Potentialgleichung 


(a — 9) 9 Pu 2m fr 0 ::- meer (4.1) 
beschreiben, indem nämlich 


ala, Sb IS run tee Ban ee) 


ist. Vermittels einer Legendre-Transformation geht (4.1) über in 
A Da N a 


= a En re Rat, Koran ARTE TE 4.3), 
während sich (4.2) in 
Eee Pi a NE CR 
Ay Ag 


übersetzt [8] (S. 161). Die Gl. (4.3) hat gegenüber (4.1) den Vorzug der Linearität. Andererseits 
erfordert die Erfüllung der Randbedingungen bei (4.3) im allgemeinen besondere Überlegungen. 
5* 
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ame ist die allgemeine Lösung von (4.3) von Hadamard in ge- 


Für Gase mit y = 


2n+l1 ö ’ 
schlossener Form angegeben worden, vgl. etwa [8], (S.164). Fürn=2,y= 5 ist 
von + MW _L@ + MW 
Tarp + 
mit den beiden willkürlichen Funktionen L(A) und M(u). (4.4) wird damit zu 
| u 292) —=MW) LU'W)—M"(w) =; 
is (+ m® (+ m? SEEN 
a _ 192) + MW EM + MU), 270) + M"(w 
a Aue + m° + m)? 


Geben wir als Anfangsbedingungen vor, daß für & = 0 die Geschwindigkeit u verschwinden soll, 
während dort a eine willkürliche Funktion von 7 ist, also 
aa u 
Ee=(: A, Teen) rn 
so folgt aus (4.5) 
LA) —M'@ 
L’(A) PER M''()) ( ) 7 (A) = 


LA) +M'@) , „UM +MO) _8,, 
en +3 ee 


Von diesen gewöhnlichen Differentialgleichungen brauchen wir nur eine partikuläre Lösung 
etwa 


0, 


L"(4) +M"() — 3 


L() = MO) = Aly(A) — A L(2) 

mit 
A j 
2 > 2 ee 
LM) =— 5 | vgow)di, L,(A) = — 5 | v? g(v) dv, 

Ar 2# 
wo die im übrigen willkürliche Integrationskonstante A* etwa der Wert von/ fürr = 0:0 = g(A*) 
ist. Das Strömungsfeld wird gemäß (4.5) durch 
LA) — (u) 


E—- A — u)T= 42 g(A) — u? gu) Ida) — Il) 


De ae NE ee = 
2270) 4 u Hu) ee en 
ET a na 


beschrieben. 
5. Iterative Bestimmung von &(4,) und g(A) 


Wir kennen nun weder die Stoßlinie, d.h. »(A,) als Funktion des charkteristischen Para- 
meters A, der ankommenden Welle, noch die Verteilung von A auf der r-Achse, d.h. die Funk- 
tion 9(A). Wir kennen nur die Anfangswerte o, für r = 0 aus (3.4) und den Wert A* =}, ,, den 
wir aus der vierten Gleichung (3.5) für = w, und A, = A, entnehmen. 

Zunächst wollen wir annehmen, wir kennten &(A,) und wollen hierzu 9(A) bestimmen 
(Bild 2). Wir berechnen nach (3.3) und (2.3) die Stoßlinie S, die die ankommende Welle be- 
grenzt, also die Funktionen &, = &,(A,), 7, = r,(A,) des charakteristischen Parameters A, der 
ankommenden Welle. (3.5) liefert dann auch die Zustandsgrößen A,(A,), u,(4,) unmittelbar nach 
dem reflektierten Stoß. 

Angenommen, wir hätten eine Näherung g(®%(A) von g(A). Indem wir in (4.6) für g(A) die 
Funktion g(®(2) und für /, u die obigen Funktionen A,(A,), u,(4,) einsetzen, erhalten wir in dem 
durch g‘%(A) gegebenen Strömungsfeld eine Kurve K: = &4), y = t.(A,). Wäre K mit S 
identisch, dann wäre g(A) = g(%(A). Ist das aber nicht der Fall, so sind K und S durch gleiche 
Parameterwerte/, aufeinander bezogen. Nehmen wir an, daß die Charakteristiken, die von 
einem Punkte P, von S ausgehen, angenähert parallel zu denen im Strömungsfeld g(%(A) sind, 
die von dem entsprechenden Punkte P, von K ausgehen (auf jeden Fall sind die Werte von A 
und « unmittelbar nach K gleich den Werten, und u, unmittelbar nach S), so können wir 
eine weitere Näherung g(%(A) von g(A) erwarten, indem wir die Charakteristiken von P, die 
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wir ja kennen, parallel nach P, verschieben und für £ = 0 in diesem Strömungsfeld 7 = g"(A,) 
setzen. Sind 7; = g%(u,) bzw. 7; = g%(A,) die Werte von r, in denen die Charakteristiken 
C- bzw. C+ von P, aus die r-Achse schneiden, so 
gilt also für die so bestimmten entsprechenden 
Wertes. %: 


I EA ee Or 

B (5-1). 
y-ur E, F—%); = gm) 

k 


Wir setzen also als nächste Näherung 

Le A en IHN 
Konvergiert das Verfahren gegen eine Funktiong(A) 
und ist die Stoßlinie S richtig, dann muß im Grenz- 
wert von selbst neben =; =g(4,) auch 
u = u = Au.) gelten. 

Ist S nicht die richtige Stoßlinie, so werden 
wir also die erste Gl. (5.1) zur gleichzeitigen Kor- 
rektur von S heranziehen. (Das entspricht der 
Tatsache, daß für den Verlauf des Stoßes nicht 
die Verteilung von A auf der r-Achse im ganzen 
Bereich 0<r<1 verantwortlich ist, sondern 
nur in dem verhältnismäßig kleinen Bereich die- 
ser Achse, in dem die C- vom Stoß aus diese 
Achse schneiden.) 

Die erwähnte Korrektur kann dadurch ge- 
schehen, daß wir, nachdem wir 941))4,) aus (9.2) Bild 2 
ermittelt haben, aus 

ee (On) en. ee) 


ein neues ı,, als u, nach dem Stoß berechnen und sodann aus diesem u, nach der letzten Gl. (3.5) 
mit u, als u, ein neues » = wU(A,) bestimmen. Die Stoßlinie selber ändert sich dabei erfahrungs- 
gemäß sehr wenig, genau so wie das zu @o(®(2,) nach der vierten Gl. (3.5) berechnete An, wie 
schon im Anschluß an (3.5) betont wurde. 
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Insgesamt ergibt sich also folgendes Schema für die Iteration: 
Gegeben ist eine Näherung o%(A,), geawm). 


a) w(®(A,) liefert nach (2.3) und (3.3) die Stoßlinie S mit &,(A,), T(Ay)s Any)» Un(Av)- 
b) 5%) gibt uns nach (4.6) mit A, und 4, als A und u die Kurve K: £,(A,), Ty(A,)- 


c) Nach der zweiten Gl. (5.1) und (5.2) bestimmen wir die nächste Näherung g®(A) von g(2). 


d) Nach der ersten Gl. (5.1) und (5.3) sowie mit dem unter c) ermittelten g(’(A) bestimmen 


wir u, und hieraus nach der letzten Gl. (3.5) w®(A,). Das Iterationsschema ist in Bild 3 dar- 


gestellt. 


P-Druck an der Wand 
—o Kontrollrechnung für /m_ P 
2 77 


2 


Bild 4 


Su 


2,5 


Bild 5 


Als Ausgangsnäherung können wir o%(},) = @, 
— const nach (3.4) und g(®(}) aus einer vorge- 
gebenen angenäherten Druckverteilung für &= 0 
in der reflektierten Welle wählen, etwa 


en 
” 2,5 BD 3 Di r /@) a ET oe rn (5.4) 
7 mit demselben /(r) wie für die ankommende Welle, 


j d.h. einfach den entsprechend den Maximaldrücke 
affin verzerrten Verlauf des Überdruckes. Dabei ist der reflektierte Maximaldruck p De 
m 


Pm _ Puh) 
De DB: = P(w,) P(®,) bestimmt. Nach (2.2) ist dann ähnlich wie (2.4), mit u =/—u=0 


ey 


Y— ap = 


} En 
ER (Oo)! Bar %) 2 a 


IN 


P) 


Br warn, 
? % Aue 
 rechtfe igte 


Be u "$ lg 5 V, 
asselbe (isentropisch) nach der Charakteristikentheorie berechnet, 


" sich merkbarer Unterschied ergab. Die Stoßgeschwindigkeit der reflektierten Welle 


\ itgehend konstant, so daß sich die isentropische Rechnung hierdurch nachträglich 
. Bild 6 zeigt den Impuls der reflektierten Welle. Verstehen wir unter „Impuls“ 
2 j DIUHRN A 


beide Male, für die ankommende und für die reflektierte Welle, den Ausdruck J z | = d S 
Erg : 20) 
irn m e £ _ £ » to > x 
so zeigt sich, daß für J als Impuls der ankommenden und J als Impuls der reflektierten Welle 
sehr genau we 
E I _Pn 
J J _Pm 
gilt. 


6. Das inverse Problem 


Er wollen jetzt umgekehrt aus dem Druckverlauf der reflektierten Welle für E=0 die 
ankommende einfache Welle berechnen. Dieses Problem stellt sich bei der Bestimmung der 


 auftreffenden Welle aus dem etwa durch ein Kondensatormikrophon gemessenen Druckverlauf 


der reflektierten Welle. Per 
Es sei also p,, gegeben, wodurch die Funktion f(r) nach 


I) = 5, ee ie ir (6.1) 


bestimmt ist. Wir müssen jetzt zunächst &, und @, ausrechnen: Es ist 


P(o,) P(o,) = Er 


und vermöge der Reflexionsbedingung (3.4) 


Po) =3-2 —, VG % 2) a I 1 T. r ee all (6.2). 


1 | Ba) Y+D@y—D .Pn | 5] 
re + Pr 
BEN Y se 
Die Funktion 9(4) ist bekannt und durch j 
| e 


re 1 
DA 0) 


[ 2 Pr Te) ER E ID FUN (6.3) 


bestimmt. 

Hiermit liegen auch die Anfangswerte A, und nach (2.1) im Punkte &e=0,7=0 fest. 
Gesucht ist nun &(A,) und g(A,). Angenommen, wir hätten eine Näherung »®(A,) von @(A,) 
und eine Näherung g(®(A,) von g(A,). Mit g(%(A,) und w®(A,) können wir nach (3.3) und (2.3) 
& = &,(4,) und 7, = T,(A,) bestimmen. Diese hierdurch bestimmte Kurve S$ ist i. a. nicht die 
Stoßlinie. Berechnen wir nach (3.5) A, und u, als Funktionen von A, und tragen diese Funk- 
tionen A,(A,) und u„(A,) als A und u in (4.6) ein, so wird sich nicht wieder die Kurve S, sondern 
eine Kurve K ergeben, die mit S nicht identisch ist. Längs S sind aber jedenfalls die Stoß- 
beziehungen erfüllt, längs K dagegen die Beziehungen zwischen A, und w,, die für & = 0 die 
Randbedingung u = 0 zur Folge haben. 

Wir verschieben nun, ähnlich wie in Ziffer 5, die von einem Punkt P, von K ausgehenden 
Charakteristiken C und C* zu dem entsprechenden Punkt P, d.h. zu dem Punkt von S, der 
denselben Wert von A, und damit auch von 7}, und w,„ hat. Durch den Schnittpunkt der ver- 
schobenen C+ mit der r-Achse wird diesem Punkte P, vermittels r — g(A,) ein neuer Wert von 
7}, und, durch A, — A,(A,) ein neuer Wert von A, zugeordnet, während die Charakteristik GN 
entsprechender Weise ein neues 4, und damit naeh der letzten Gl. (3.5) ein neues © = a (A,) 
liefert. Die nächste Näherung von g(A,) finden wir, indem wir die Charakteristik GE in der an- 
kommenden Welle von P, aus konstruieren und diese mit der r-Achse schneiden (Bild 3), d.h. 


durch 


a ta, vr] 3—y 
gda)=T ne rn. = a at (6.4). 
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Daraus ergibt sich folgendes Iterationsverfahren: 
Gegeben ist eine Näherung o%(A,) und g%(A,). ef 
a) wO(A,), g’(A,) liefern nach (3.3) und (2.3) die „Stoßlinie‘“ S mit &,(A,); TA); And); 
M A ) i ” . . . 
je b) Die Funktionen A,(A,), An(A,) ergeben, als A und u in (4.6) eingesetzt, eine Kurve K mit 
&A,) TelA,). 
ec) Durch 
PA & H va = Fi 
rt = A > 
k 


und die vierte Gl. (3.5) mit w%(A,) , A, und A, als @, A, und A, ergibt sich eine neue Zuordnung 
A,(35). g®(A,) ist dann durch (6.4) mit A, als A, gegeben. 
d) Durch 


Ku) nt im). tee 


gewinnen wir zu A,, A, ein neues i,. Die letzte Gl. (3.5) gibt uns sodann o®(A,). 

Als Ausgangsnäherung können wir wieder o®%(A,) = @, = const nach (3.4) nehmen. Für 
g‘®(A,) wählen wir wieder die Verteilung von A, auf der r-Achse, die dem den Maximaldrücken 
entsprechend affin verzerrten Überdruckverlauf entspricht, also nach (2.4) 


! 


Pi PR, Es 
14), BP) —1. 


Die praktische Erprobung ergab genau so gute Konvergenzverhältnisse wie bei dem in 
Ziffer 5 beschriebenen Verfahren. 


2y 
y—1 a! 
a en 


7. Allgemeine auftreffende Welle 


Obwohl die Annahme w„, = const = u 
in der auftreffenden Welle i. a. nicht nur gut 
erfüllt ist, sondern auch, solange man den 
Ursprung, d.h. die speziellen Anfangsbedin- 
gungen, dieser Welle nicht kennt, nicht weni- 
ger plausibel als irgend eine andere Annahme 
ist, sei hier noch der Fall skizziert, daß wir für 
ein & < 0 sowohl A, als auch u, als Funktionen 
von r kennen, womit natürlich auch der Verlauf 
des Stoßes für & <<< 0 gegeben ist. Die An- 
nahme «, = const ist dann durch die Vertei- 


lung A, = g(T), #, = h(r) in einem gewissen Be- 
reich der Geraden & = £&, der nach oben und 
unten durch den ankommenden und reflek- 
tierten Stoß begrenzt ist, ersetzt und präzi- 
siert (Bild 7). Die Rechnung ist dann nicht 
wesentlich komplizierter: Die Strömung, in 
dem durch die Verteilung von A, und u, für 
& = £, bestimmten Strömungsfeld kann analog 
wie in Ziffer 4 berechnet werden, und liefert 
an der Stelle &=0%=g(r), #u=h(r) und 
damit p’/pı = f(r). Das Iterationsverfahren 
Bild 7 von Ziffer5 ist dann auch für die ankom- 
mende Welle, also für die Strömung vor dem 
Stoß anzuwenden. Da prinzipielle Schwierigkeiten hierbei nicht mehr auftreten, sei hier auf 
die detaillierte Durchführung verzichtet. 
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KLEINE MITTEILUNGEN 


Interpolations- 
polynomen. 
Üblicherweise wird die Differenz zwischen einer reell- 


wertigen Funktion f(x) der reellen Variablen x und ihrem 


zugehörigen ‚Interpolationspolynom P(x) unter Be- 
nutzung einer höheren Ableitung von f(x) abgeschätzt. 
Da diese höhere Ableitung aber oft nur schwer zu be- 
rechnen ist, soll sie im folgenden unter einer bestimm- 
ten Voraussetzung ersetzt werden durch die entspre- 
chende höhere Differenz. 

Es sei B=&%Htih = 0,0, Ion.s st el 
Weiter sei f(x) eine (r + 1)- mal stetig differenzierbare 
Funktion im Intervall [x,, &r+1] und es sei P(x) das 
Interpolationspolynom r-ten Grades (r>1) mit 
Fa) ta) h fürrs=0, 1,0. ., 7). Wir setzen 

Ri; —h und entsprechend die höheren Dif- 
ferenzen Vkf;. 

Satz 1: Wenn die (r + 1)-te Ableitung /r+D(x) im 
Intervall [x,, &+1] stets > 0 oder stets S 0 (also ohne 
Vorzeichenwechsel) ist, so gilt für, SS ar: 


If) — P(@)| 
„ (em): —%)'- 
r!. m 
Beweis: Es sei P,(x) das Polynom (r — 1)-ten Grades 
mit P,(&;) en, 2.57. Dann ist 


di .» (2 — 27) R vr+Hf Er, 


Pie) = Pia) + pa). ea Prhr 
und für 
P,(«) = P&@a)+ 7, 2)... (a) her 


ist Pij= = fhfüri=1, 2,..., r+ 1, wie sofort aus 
der Newtonschen Interpolationsformel für aufstei- 
gende bzw. absteigende Differenzen folgt. 


Es ist nun bekanntlich für z, sx = a;: 


ee) 
(X == & = Xr) ’ 
H&) — Py() 
= (x — %rH) !! a Er Dr L an) Oo) 


M<N<Lr): 


Ersichtlich haben die rechten Seiten verschiedene Vor- 
zeichen und es ist also entweder 


Pix) = (@) = Pl) 


P,(«) < f(x) < Pe) 
Also ist sicher 


oder 


Ita) — Pla)] < |Ps(a) — Pla) 

und schließlich ist . 

Pla) — Pla) = 1 ea) (am Wh Ph} 
EN 2 Bee 


Wegen seiner Bedeutung für die numerische Inte- 
gration soll diese Methode noch auf das Stirlingsche 


Interpolationspolynom dritten Grades angewendet 


werden. Es sei f(x) eine viermal stetig differenzierbare 


Funktion im Intervall [x,, x,] und es sei Q(x) das Poly- 


nom zweiten zus mit 9(x;) = ifüri = 1,2, 3. Es 
wird 

ge) = — 2) (&— 8) (025) VPfzrr 
(k=1, 2) 


Ur) + le 


gesetzt. 
Q: (“)= fh für ı = l, 2, 3, 4 


Satz2: Wenn fIV(x) im Intervall [x,, 24] stets > 0 
oder m < 0 ist, so gilt fürx, Ss <S%: 


IN 


(a — %ı) (& — 2%) ( — 2) V*fıl 


. 
1 

<—Y Sale 

= 3 W/ fa 


Beweis: PB ist 


fe) — Qıla) = ea) («— 2) (&—%,) (© 25) VE) 
und 
I) — Qya) = 2; @— 2) @- 2) (e—) ea) WVn). 


Da die rechten Seiten verschiedenes Vorzeichen haben, 
folgt daraus entweder Qı,(x) < f(x) Ss Q;(x) oder 
0,2) <a) SQ). Also ist 


+] - 


und es ist 


Qua) — Qua) = m 


Durch die Berechnung der Extrema sieht man schließ- 
lich, daß 


1 
anlle—a) (2 — X) (C — %,)| S 7773 


a —% a) — 0) 


— 2) @—%) @—%) Ph. 


ist. 
Aus Satz2 läßt sich unmittelbar eine Fehlerab- 
schätzung für die Simpson -Regel gewinnen. 


Satz 3: Wenn fIV(x) im Intervall [x,, x,] stets > 0 
oder stets < 0 ist, so gilt 


sn IM *Fa 


Staa + eh +) SP. 


Beweis: Mit den Bezeichnungen von Satz 2 ist 


%g %a 
h 
| Ber (Aade=zhtthth: 


ad Ds 


Dann ist Q,a)=fhfürı=0,1,2,3und 


EURUOE ee; 


er % 


1 


Werner Uhlm ann. 


Eine Verallgemeinerung des Hornerschen 
Schemas 


Das Horner:sche Schema wird mit Vorteil benutzt, 
um Polynome P(x) zu behandeln, die in der Form 


- 


gegeben sind. Dabei sind die Polynome P,(x) die Po- 
tenzen von x, bei einer allgemeineren Taylor-Ent- 


P(a) - 3,4 De 


wicklung die Potenzen von 2 —x,, bei einem New-. 


tonschen Interpolationspolynom noch allgemeiner die 
. ” 

Produkte II («—x,). In allen diesen Fällen gilt 
1 


= 
(gegebenenfalls mit einander gleichen x,) ein System 
von Rekursionsgleichungen 

mt 9 


Pa) =1 
Fe) (0 2,) Br 1a m 1,2. 

Das Horner-sche Schema besteht in der Durchrechnung 

der Rekursionsformeln 


Ba(‘)—An Ye) 
B&)=Ar + (E—%r+1) Br+1(d); (r=n—]1,...,0) ; 


Löst man diese nach den A, auf und führt sie in (1) 
ein, so wird wegen der Formeln (2) 


Pie) = Bl) +@—HE Bi) Pre)... (Mi 
daraus folgt x 


Pi Bien Se (5) 
und 
Bi) BR n 
= 559 Pre) Ka”, (6). 
v1 


P(£) wird also mit dem Schema (3) bestimmt, ohne 
daß die P,(£) einzeln bekannt sein müssen. Die rechte 
Seite in (6) ist wieder von der Form (1), so daß das 
Schema in bekannter Weise unmittelbar fortgesetzt 
werden kann. 

Diese Vorteile des Hornerschen Schemas können 
auch dann noch ausgenützt werden, wenn die Poly- 
nome P,(x) in (1) dreigliedrigen Rekursionsgleichungen 


N) 
Po@) = (2’) 
Pr) = (a — br) Pr-ı(®) + Cr Pr—2(R); 


Wand) 


genügen. Die bekanntesten Polynomsysteme mit einem 
solchen Bildungsgesetz sind die bezüglich einer Ver- 
teilungsfunktion orthogonalen Polynome. Diese und 
weitere Polynomsysteme der Art (2’) treten auch als 
Näherungsnenner Stieltjesscher Kettenbrüche auf. 
[1], [2]. 

Den Rekursionsformeln (3) entsprechen nun die For- 
meln 


Bn(E) — An 
Bn-ı(E) —= An-ı + (An & — bn) Bn(E) 
Br(£) —= Ar, + (ar4ı E— br-+1) Br+1(&) 
+ er+2 Br+2(&) (r=n—2,...,0) 


(9). 


Sp: em) e—a) @—aldr un 
s % = 


Die Formeln (3°) können in 
werden, daß das Schema 


va)  \, 
An f 


Ani 


Bn | | vb Bnsarkit “ 


‘+ 3.1 R Da D- 

L Ä MITTE SZ 
Spalte für Spalte ausgefüllt wird. Für manche pr 
ziellen Polynomsysteme wird dieses Schema einfache; 
z.B.sind für Hermitesche Polynome in der üblichen 
Normierung die a, einander gleich, während die br 
verschwinden. — Die Ähnlichkeit des Schemas mitdem 
von L. Collatz [3] angegebenen ist nur formal. 


W. Everling. 


Ba—ı T 


SS 


Aachen. 
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Stabilität der Gleichgewichtslage in einem 
nicht-holonomen System 


Bei Untersuchungen über die Stabilität der Gleich- 
gewichtslage eines nicht-holonomen Systems „im 
Kleinen“ linearisierte Whittaker!) von vornherein 
die nicht-holonomen Bedingungsgleichungen und be- 
hielt in ihnen nur die Glieder erster Ordnung bei. 
Danach integrierte er diese Gleichungen und setzte 
überdies, indem er die Gleichgewichtslage in den Ko- 
ordinatenursprung verlegte, die Integrationskonstante 
gleich Null. Auf diese Weise führte er das Stabilitäts- 
problem der Gleichgewichtslage eines nicht-holonomen 
Systems auf das entsprechende Problem eines holo- | 
nomen Systems zurück. ’ 

Bottema?) wies auf die Fehlerhaftigkeit der Über- 
legungen Whittakers hin. Er erklärte, daß man bei a 
richtiger Aufstellung der linearen Näherungsgleichungen 
für nicht-holonome Systeme auch die von Whittaker 
vernachlässigten Glieder zweiter Ordnung in den nicht- 
holonomen DBedingungsgleichungen berücksichtigen 
müßte. 

Er zeigte, daß die richtig aufgestellte charakteristische 
Gleichung der linearen Annäherung immer eben so viele 
verschwindende Wurzeln besitzt, wie das betrachtete 
System nicht-holonomen Bedingungen zu genügen hat.- 

Die charakteristische Gleichung, die man als Ergebnis 
der Überlegungen Whittakers erhält, hat diese ver- 
schwindenden Wurzeln nicht. Außerdem können sich 
die übrigen, von Null verschiedenen Wurzeln der rich- 
tigen charakteristischen Gleichung sogar durch das Vor- 
zeichen des Realteiles von den Wurzeln der charak- 
teristischen Gleichung Whittakers unterscheiden. 

Im Endergebnis blieb also die Frage nach der Sta- 
bilität eines nicht-holonomen Systems offen, weil bis 
jetzt die Bedingungen für die Stabilität von Systemen, 
bei denen die zur linearen Näherung gebildete charak- 
teristische Gleichung eine beliebige Zahl von ver- 
schwindenden Wurzeln hat, im allgemeinen Fall immer 
noch nicht gefunden worden sind. Die Aufgabe läßt 
sich jedoch zu Ende führen, weil man sie auf einen 
Spezialfall zurückführen kann, der vollständig von 


FR 


_ 


a 
e 


f und I. 6. Malkin untersucht wor- 


Zu ‚diesem Zweck betrachten wir klekıeı Schwin- 
gungen eines nicht-holonomen mechanischen Systems 


mit den Koordinaten x, 2... Yı> Ya» --- m in 


bezug auf die Gleichgewichtslageg =: --=m = 


=: > 


Y . 5 
Das System ist r differentiellen (integrablen oder 


nicht-integrablen) Bedingungen 


RR 


a a Yan  . 4 


rn 


n 
Ya= 2 Aajlaı.. an, Yı =» Yn) &; “Br. (1) 


= 
(Gl. or) 
unterworfen. 
Die Appelschen Gleichungen dieses Systems geben 
wir in der Form 


& — Fi, 2%,» ns In» Yı: +: Yn) \. (2) 

=1,2...n) ; 
an. 
An Stelle von y. führen wir neue Veränderliche ua 


N 
U = Ya — 2 As; ae else (3) 
= 


ein; hierbei sind die A%; die Werte der Funktionen 4a 
für verschwindende Argumente. 

Durch die neuen Veränderlichen werden die Glei- 
chungen (1) auf die Form 


u 2 Aa Aal &% | EEE (4) 
(& — % See r) 


gebracht, und die Reihenentwicklungen der rechten 
Seiten von (4) fangen mit Gliedern an, die bezüglich 
der Koordinaten und Geschwindigkeiten von mindestens 
zweiter Ordnung klein sind. 

Die Gleichungen (4) bilden zusammen mit den Glei- 
chungen (2), wenn wir in diesen die ya vermöge (3) 
durch die u« ausdrücken, ein System, für das die lineare 
Annäherung ein charakteristisches Polynom mit r ver- 
schwindenden Wurzeln besitzt. In diesem System ent- 
halten die ersten r Gleichungen rechts keine linearen 
Glieder, und deswegen kann man darauf das Theorem 


’s 


uchbesprechungen 


von Ljapunoff-Malkin anwenden, welches besagt, 
daß ein solches System stabil (jedoch nicht asymptotisch 
stabil) ist, wenn die übrigen der n Wurzeln der charak- 
teristischen Gleichung (außer den r verschwindenden 
Wurzeln) negative Realteile haben. 

Also ist die Gleichgewichtslage eines nicht-holonomen 
Systems stabil (mur nicht asymptotisch stabil), wenn alle 
Wurzeln (mit Ausnahme der oben aufgezeigten ver- 
schwindenden Wurzeln, deren Anzahl mit der Zahl der 
nicht-holonomen Bedingungen übereinstimmt) negative 
Realteile haben. 

Mit anderen Worten: die von Bottema entdeckten 
verschwindenden Wurzeln können bei Stabilitäts- 
untersuchungen weggelassen werden, jedoch garantiert 
das negative Vorzeichen der Realteile der übrigen 
Wurzeln nicht die asymptotische, sondern nur die ein- 
fache Stabilität. 

Alles Gesagte bezieht sich nicht nur auf nicht- 
holonome, sondern auch auf semi-holonome Systeme, 
d. h. auf Systeme, die zwar differentielle, aber integrable 
Bedingungen enthalten. Wenn die semi-holonomen Be- 
dingungen bei der Aufstellung der Bewegungsglei- 
chungen und ihrer Untersuchung als holonome be- 
trachtet werden können, die eine willkürliche Konstante 
als Parameter enthalten, so muß man mit den semi- 
holonomen Bedingungen bei Stabilitätsuntersuchungen 
genau so verfahren wie mit nicht-holonomen. 
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BUCHBESPRECHUNGEN 


Ing. B. v. Schlippe, Strömung von Flüssigkeiten 
mittemperaturabhängiger Zähigkeit (Kühlung 
von Ölen). (Forschungsberichte des Wirtschafts- und 
Verkehrsministeriums Nordrhein-Westfalen, Nr. 255) 
448. m. 12 Abb. Köln/Opladen 1956. Westdeutscher 
Verlag. Preis brosch. 11,70 DM. 

Läßt man in ein genügend langes geradliniges Rohr 
eine ölige Flüssigkeit mit einer vorgegebenen Geschwin- 
digkeitsverteilung etwa von parabolischer Gestalt ein- 
strömen, so kanı die Temperatur dieses Öls nicht 
unter einen bestimmten Minimalwert herabgedrückt 
werden, wie stark man auch die Rohrwand abkühlt. 
Dieser Effekt beruht darauf, daß mit sinkender Tem- 
peratur der Rohrwand die Zähigkeit des Öls in den 
wandnahen Zonen zunimmt-und dadurch die Ent- 
wicklung von Reibungswärme gesteigert wird. V£. 
untersucht für eine laminare Strömung die Geschwin- 
digkeits- und Druckverteilung über dem Rohrquer- 
schnitt, die sich im ausgeglichenen Zustand einstellen 
und berechnet die Minimaltemperatur und den Druck- 
verlust im Rohr. Zur Vereinfachung wird für die 
Zähigkeit in Abhängigkeit von der Temperatur ein 
Potenzgesetz der Form n — 1/0* zugrundegelegt und 
das entstandene System von Integralgleichungen für 
den Wärmefluß und die Öltemperatur graphisch gelöst. 
Außerdem wird das Verfahren auf die ebene Strömung 
zwischen zwei parallelen Platten übertragen. 


Halle (Saale). H. Schubert. 


W. Blaschke, Hamburg u. H.R.Müller, Graz- 
Ankara. Ebene Kinematik. (Mathematische Einzel- 
schriften, Hamburg, Bd. 5). 269 S. m. 100 Abb. Mün- 
chen 1956. Verlag R. Oldenbourg. Preis geb. 26,80 DM. 


Das gesamte Gebiet der ebenen Kinematik wird von 
den Verfassern in den vier Abschnitten: 
Zwangläufige Bewegungsvorgänge im Kleinen, 
Zwangläufige Bewegungsvorgänge im Großen, 
Flächenläufige Bewegungsvorgänge, 
Kinematische Abbildung 
behandelt, wobei für die Darstellung ausschließlich 
Vektoren benutzt werden. Für die Bewegungsvorgänge 
werden „Ableitungsgleichungen‘ gefunden, in denen 
lineare Differentialformen in t, sogenannte „Pfaffsche 
Formen“ auftreten. Diese Schreibweise zeigt sich 
besonders bei Betrachtungen mehrgliedriger Bewe- 
gungsvorgänge als zweckmäßig. Ausführliche Behand- 
lung erfahren die Konstruktion der Krümmungskreise 
von Bahnkurven, die Bresseschen Kreise, die Hüll- 
bahnen die zyklischen Kurven, die Darstellung ebener 
Bewegungsvorgänge mittels komplexer Zahlen, die 
Koppelkurve des Gelenkvierecks und die Geradfüh- 
rungen. Bemerkenswert sind weiter die Betrach- 
tungen über Bahnflächen, Polarplanimeter, Schwer- 
punkte der Punktbahnen, Hüllbahnflächen und Krüm- 
mungsschwerpunkte der Hüllbahnen. Den dem ersten 
Abschnitt zugrunde gelegten eingliedrigen Bewegungs- 
vorgängen werden im dritten durch Hinzunahme eines 


r 


Bewegungsvorgänge gegenübergestellt. Einen be- 
ee an an at die kinematische Abbil- 
dung ein, die 1911 von W. Blaschke und J. Grünwald 
eingeführt und 1923 von E. Kruppa verallgemeinert 
wurde. Es ist anzunehmen, daß die Betrachtungen 
dieses Abschnittes wertvolle Anregungen für die Be- 
handlung der räumlichen Getriebe geben werden, und 
man kann auf den angekündigten zweiten Band, der 
der räumlichen Kinematik gewidmet sein soll, sehr 
gespannt sein. Das vorliegende Buch ist für das Studium 
der ebenen Kinematik bestens geeignet. 


Dresden. W.Lichtenheldt. 


W. Schiebeler, Luftströmungen mit Drall im 
Kreisrohr hinter radialem Leitapparat. (Mit- 
teilungen aus dem Max-Planck-Institut für Strömungs- 
forschung, Nr. 12.) 101 S. mit 66 Abb. Göttingen 1955. 
Selbstverlag Max-Planck-Institut für Strömungsfor- 
schung. Preis brosch. 9,— DM. 

Die vorliegende Veröffentlichung bringt umfassende 
Versuchsergebnisse über rotierende Luftströmungen in 
einem Kreiszylinderrohr ohne Nabe (Durchmesser 
12,2 cm). Zur Erzeugung der Drallströmung dient ein 
Leitapparat, durch den die Außenluft dem Meßrohr 
zuströmt und in den unverwundene oder verwundene 
Leitschaufeln eingesetzt wurden. Ein besonderer Leit- 
apparat dient noch der Untersuchung der Vermischungs- 
vorgänge beim Zusammenfließen zweier verschiedener 
Drallströmungen. In allen Fällen handelt es sich um 
Strömungen, die neben der axialen Geschwindigkeit 
komponente auch eine solche in Umfangsrichtung 
enthalten. Bisweilen auftretende Komponenten in 
radialer Richtung dürfen als zu gering vernachlässigt 
werden. Druck-, Winkel- und Temperaturmessungen 
dienen der Ermittelung der Geschwindigkeitsvertei- 
lungen. 

Die Messungen der einfachen Drallströmungen 
zeigen durchweg in Achsennähe einen Kern mit sehr 
geringer Durchflußgeschwindigkeit und starker Er- 
niedrigung der Staupunktstemperatur die von der durch- 
fließenden Rotationsenergie abhängt, während Kern- 
durchmesser und Drallwinkel an der Rohrwand eine 
eindeutige Beziehung zu einander aufweisen. Die Ver- 
mischungsversuche lassen erkennen, daß ein Ausgleich 
der anfänglichen Unterschiede nach Durchlaufen einer 
Länge von zwei Rohrdurchmessern erfolgt ist. 

Eine ausführliche Besprechung gilt noch Vergleichen 
mit verschiedenen anderweitigen rechnerischen und 
versuchsmäßigen Behandlungen des Problems, darunter 
auch den Beziehungen zum „Wirbelrohr‘“‘. 


Weilburg. L. Schiller. 


Tables of Weber parabolic ceylinder func- 
tions. Computed by Scientific Computing Service 
Ltd. Mathematical Introduction by J.C. P. Miller. 
234.8. m. 19 Abb. London 1955. Her Majesty’s 
Stationery Office. Preis geb. 63 s. 


Der erste Abschnitt des vorliegenden Bandes enthält 
einen umfassenden Abriß der Theorie der Differential- 
gleichungen 


v—(La+a)y=0 PN (1) 


und 


Hier werden für beide Dgln. Hauptsysteme von 
Lösungen — U(a, x) und V(a, x) für (1), W(a, x) und 
W (a,—x) für (2) — angegeben. Es folgen Integral- 
darstellungen, asymptotische Entwicklungen, Hilfs- 
funktionen, Beziehungen zu den Besselschen, den 
konfluenten hypergeometrischen und anderen Funk- 
tionen sowie Formeln für die Gammafunktion. 


weiteren (Parkmeiere zweigliedrige oder flächenläufige 


Tp 
Formeln bringt — getrennt für be 
Abschnitt. e TREE 

Tabellen für die Lösungen von (2) we 
im vierten Abschnitt veröffentlicht. 
wesentlichen: 2 | i 

Tab.I: W(a,x) und W(a,—x) (mit i.a. 


Pas, 
& 


geltenden Ziffern) sowie deren reduzierte Ableitungen 3 
N N A = En Be) 
a — a Car W(a, + x) (außer beit mit einer Dezimal- 
n!daen Dr Dir WAR 


a 5a 
stelle mehr als W) . ee - 
a= —10 (1)10;2 = 0(0,1)1;r=1(1)9. 

Tab. II: In W(a, a und u (siebenstellig) mit 
modifizierten 2. und 4. Differenzen 
a=1(1)101; +2=0(0,1)L, L>22Ya. 

Tab. III: Hilfsfunktionen F(a, x), x(a,x), G(a, x) 


und y(a, x) (achtstellig) mit modifizierten 2. Differen- £- 

zen =. * 
a=—10(1)2; x = 0 (0,1)10 > ;; 

und j 


«=3(1)10; <= 20.1)70, 22 eye 


Zwischen W(a, x), W(a,—x) und diesen Hilfsfunk- 
tionen bestehen die folgenden Beziehungen: 


W(a,x) =\kFcosx, W(a, —x) = = F sin x, 


W'a,x)=—\kGcosy, W’(a, —x) = ee : 


Vr 
k=|\1+e2ra—ena, 


Dresden. A. Schubert. 


Dr. W. Gröbner (o. Prof. a. d. Univ. Innsbruck), 
Matrizenrechnung. (Mathematische Einzelschriften. 
herausgegeben von W. Blaschke, Bd.4) 249 S. m. 
9 Abb. München 1956. R. Oldenburg Verlag. Preis 
geb. 23,— DM. 

Mit seinem Buch über Matrizenrechnung legt der 
Verfasser, der sich auf seinen Lehrer Furtwängler als 
Vorbild beruft, eine überaus reizvolle und interessante 
moderne Darstellung der Theorie der endlichen Ma- 
trizen vor, in die er in begrüßenswerter Weise die 
Vektoralgebra mit einbezieht. Nach knapper Ein- 
führung der Grundbegriffe Vektor und Matrix und der 
einfachsten Tatsachen der zugehörigen algebraischen 
Kalküle wird vom Graßmannschen Begriff des äußeren 
Produktes aus in geradezu überraschender Weise eine 
vollständige Determinantentheorie aufgebaut. An 
gegebener Stelle wird mit Hilfe der zur Verfügung 
gestellten Sätze die Algebra der linearen Gleichungs- 
systeme und Transformationen einschließlich ihrer 
Anwendungen in der Geometrie (affine und projektive 
Geometrie) behandelt. In einem Paragraphen über 
rationale Funktionen einer quadratischen Matrix, der 
bis zu den Frobeniusschen Kovarianten einer qua- 
dratischen Matrix vordringt, wird, was für die Anwen- 
dungen besonders wesentlich ist, ausführlich auf die 
charakteristische Gleichung und die Eigenwerte einer 
quadratischen Matrix und alle wichtigen diesbezüg- 
lichen Sätze eingegangen. Der vorletzte Abschnitt 
bringt die Theorie der Äquivalenz und Ähnlichkeit 
einschließlich der Elementarteilertheorie und der bei 
Basistransformation des n-dimensionalen Vektorrau- 
mes auftretenden Fragen, wobei von dem schon vorher 
bereitgestellten Hilfsmittel der Übermatrizen Gebrauch 
gemacht wird. Zum Schluß werden die allgemeine und 
die Hermitesche Kongruenz und ihre Anwendung bei 
der Erörterung der Hauptachsenprobleme behandelt, 
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Die Darstellung ist sehr konzentriert gehalten und 
verlangt vom Leser nicht nur ein intensives Mitarbeiten, 


sondern an vielen Stellen auch schon eine gewisse 


'Vertrautheit mit dem behandelten Gegenstand. Das 
betrifft insbesondere auch die zahlreichen Bemerkungen, 
Ergänzungen und Aufgaben, die den einzelnen Para- 
graphen angehängt sind und den Inhalt des Buches in 
apper Form um ein Beträchtliches vermehren. In 
diesen Ergänzungen wird gelegentlich auch das Inter- 
essengebiet der numerischen Analysis kurz gestreift. 
- Soweit auf Seiten des Referenten gegen gewisse 
Einzelheiten Bedenken bestehen, hat er sich mit dem 
Autor selbst in Verbindung gesetzt. Sie können in 
keiner Weise die Tatsache beeinträchtigen, daß es sich 
bei dem vorgelegten Buch um eine sehr wertvolle 
Bereicherung unseres mathematischen Schrifttums 
handelt. ; 


Dresden. H. Heinrich. 


Proceedings ofthe Conference on High-Speed 
Aeronautics 1955. Herausgegeben von A. Ferri, 
N. J. Hoff und A.Libby. VI + 3928. m. 174 Abb. 
Sa Institute of Brooklyn, Brooklyn, New 

ork. 

Anläßlich der Feierlichkeiten zum hundertjährigen 
Bestehen des Polytechnischen Institutes in Brooklyn 
fand im Januar 1955 eine Tagung über Hoch-Geschwin- 
digkeitsfragen der Luftfahrt statt. 

Die Reihe der 16 ausführlich wiedergegebenen Vor- 
träge, die meistens zusammenfassenden Charakter 
haben, eröffnet Th. vonKärmän. Er spricht über 
gelöste und ungelöste Probleme der Hochgeschwin- 
digkeitsaerodynamik. Sein Bericht, der mit einer 
historischen Darstellung beginnt, führt schnell an ak- 
tuellen Fragen heran. Kärmän, dem die Aerodynamik 
den Impulsansatz zur Lösung der Prandtlschen: Grenz- 
schichtgleichung verdankt (1921), hat sich seit Beginn 
der dreißiger Jahre vornehmlich mit den Fragen kom- 
pressibler Strömungen beschäftigt. So kann nun der 
75jährige eine Darstellung an Hand eigener Arbeiten 
oder Berichte von Mitarbeitern geben. Das Singulari- 
tätenverfahren zur Lösung der linearisierten kompres- 
siblen Potentialgleichung, das Ähnlichkeitsgesetz der 
schallnahen Strömung, die Berechnung der aerodyna- 
mischen Aufheizung beim Flug mit hoher Machzahl 
sind wesentliche Beiträge. Das letzte. Problem be- 
herrscht man für mäßige Überschallgeschwindigkeit, 
jedoch sind für sehr große Machzahlen M = 10 bis 
25 die physikalischen Vorgänge des Wärmeübergan- 
ges noch ungeklärt. Überhaupt scheint die kommende 
Forschungsarbeit engen Kontakt zwischen vielen 
Wissensgebieten nötig zu machen, während in der 
Frühzeit noch die Isolierung reiner aerodynamischer 
Probleme gelang. 

H.L. Dryden, der durch seine Arbeiten über Grenz- 
schichten bekannt wurde und heute Direktor des 
N.A.C.A. ist, gibt einen Bericht zum Thema Grenz- 
schichtübergang von laminar nach turbulent bei Unter- 
und Überschallgeschwindigkeit. Neue Versuchsergeb- 
nisse über den Einfluß des Turbulenzgrades und der 
Rauhigkeitserhebungen auf den Umschlagpunkt sind 
zusammengestellt. Im Überschallbereich nehmen Ver- 
dichtungsstöße Einfluß auf den Umschlag. Die Ergeb- 
nisse von Stabilitätsrechnungen für Überschallge- 
schwindigkeit werden kritisch besprochen. einige expe- 
rimentelle Arbeiten können zu Rate gezogen werden. 

Spezielle Probleme behandeln die Vorträge von 
L. Crocco (Zusammenwirken von Verdichtungsstoß 
und Grenzschicht), K. Oswatitsch und F. Keune 
(Rotationskörper bi M =), R. Michel und M. 
Sirieix (Experimente über Tragflügelströmung beim 
Schalldurchgang), H. Schlichting (inkompressible 
Gitterströmung; später ausführlich wiedergegeben in 
dem V.D.I. Forschungsheft Nr. 447) und A. Buse- 
mann (Lärmausbildung bei Überschallgeschwindig- 


keit). 


prechungen 


Die aerodynamische Aufheizung bei Überschallge- 
schwindigkeit stellt auch besondere Festigkeitspro- 
bleme. Diese sind Themen der Vorträge von L. Brog- 
lio, H. Horton und N. J. Hoff. Einen Querschnitt 
über die zur Untersuchung der hohen Temperaturbe- 
anspruchung beim N.A.C.A. geleistete Arbeit bringt 
J. Duberg. Nach einer elementaren Darstellung über 
die Vorausberechnung der Gleichgewichtstemperatur, 
die zu einer bestimmten Machzahl gehört und die 
durch Strahlung gemindert werden kann, wird die 
Wärmeübergangszahl untersucht. Festigkeitseigen- 
schaften der Beplankung abhängig von der Tempera- 
tur werden dargestellt und schließlich Wärme- 
spannungen typischer Bauelemente aufgezeichnet. 
R. Smelt referiert über die Möglichkeiten, Windkanal- 
versuche im Bereich sehr hoher Machzahl vorzunehmen. 
A. Libby und M. H. Bloom berichten über die Mög- 
lichkeiten experimenteller Untersuchungen und A. 
Ferri über theoretische Arbeiten beim Brooklyn In- 
stitut. A. Kantrowitz hat einen kurzen Abriß seines 
Vortrages niedergeschrieben über physikalische und 
chemische Phänomene in Gasen die auf hohe Tempe- 
raturen des Schnellfluges mit 15facher Schallgeschwin- 
digkeit aufgeheizt werden. Dissoziation, Ionisation, 
Bildung von Stickstoffoxyd, hohe elektrische Leit- 
fähigkeit sind zu erwarten. Einige Versuchsergebnisse 
sind in Diagrammen dargestellt. 


Dresden. W. Albring. 


R. Rothe, Höhere Mathematik für Mathe- 
matiker, Physiker, Ingenieure. (Herausgegeben 
von Dr. phil., Dr.-Ing. E. h. W. Schmeidler, 
o. Prof. a. d. T. U. Berlin-Charlottenburg) Teil VII: 
Räumliche und ebene Potentialfunktionen, Konforme 
Abbildung. Integralgleichungen. Variationsrechnung. 
2188. m.43 Abb. Stuttgart 1956. B. G. Teubner. Preis 
brosch. 19,80 DM. 

Mit dem vorliegenden Teil VII schließt der Verfasser 
die Erweiterung des bekannten und beliebten Rothe- 
schen Lehrbuches, dessen Herausgabe er nach Rothes 
Tod übernommen hat, ab. Aus dem Untertitel geht 
bereits die Stoffeinteilung hervor. Der erste Teil 
(85 8.) gibt eine sehr schöne Einführung in die räumliche 
und die ebene Potentialtheorie einschließlich der kon- 
formen Abbildung und der elliptischen Integrale und 
Funktionen und ihrer Anwendungen. Das Ziel des 
zweiten Teiles (77 S.) ist eine Darstellung der wichtigsten 
Tatsachen über lineare Integralgleichungen und ihre 
Bedeutung für die Behandlung von Randwertaufgaben. 
Sie wird vorbereitet durch eine determinantenfreie 
Behandlung der linearen Gleichungssysteme und einen 
Paragraphen über Vektoren und Matrizen, der in 
knapper Form bis zu der praktischen Lösung von 
linearen Gleichungssystemen und von Eigenwertauf- 
gaben durch Iteration geht. Ein kurzer Paragraph 
über die Fouriersche Integralgleichung, die Fourier- 
Transformation bei zwei Variabeln und insbesondere die 
Laplace-Transformation schließt dieses Kapitel ab. Der 
dritte Teil (54 8.) führt an Hand einfacher Probleme in die 
Gedankengänge der Variationsrechnung ein, behandelt 
die wichtigsten Variationsprobleme der Differentialgeo- 
metrie und gibt einen kurzen Abriß der Weierstraß- 
schen und der Jacobi-Hamiltonschen Theorie. 
Jedem Kapitel ist in Wahrung der Rotheschen Tra- 
dition eine Sammlung schöner, den Anwendungen 
entnommener Aufgabe mit Angabe des Lösungsweges 
angehängt. 

Wer dieses umfassende Programm mit dem geringen 
Umfang des Buches vergleicht und berücksichtigt, daß 
es bei einer solchen lehrbuchmäßigen Darstellung un- 
erläßlich ist, den Lernenden durch die Gedankengänge 
möglichst lückenlos hindurchzuführen und ihm die 
Lektüre nicht durch Auslassung wesentlicher Teile der 
Zwischenrechnung unnötig zu erschweren, wird ermes- 
sen können, welchergroßen Geschicklichkeit und Unter- 
richtserfahrung es von seiten des Verfassers bedarf, 


trotz aller gebotenen Kürze dem Auffassungsvermögen 
‚des Lernenden anzupassen. Der Verfasser hat diese 
schwierige Aufgabe vorbildlich gelöst. Im einzelnen 
ist die Darstellung selbstverständlich aufs engste mit 
derjenigen in den „‚Integralgleichungen“ und den „Vor- 
trägen über Determinanten und Matrizen‘ des gleichen 
Verfassers verwandt. 


Dresden. H. Heinrich. 


Dr. rer. nat. Klaus Pöschl. Mathematische Me- 


‚thoden in der Hochfrequenztechnik. VII + 
331 S. m. 165 Abb. Berlin/Göttingen/Heidelberg 1956. 


Springer-Verlag. Preis geb. 36,— DM. 

Das vorliegende Buch schließt in der Art seiner Zu- 
sammenstellung eine gewisse Lücke in der deutsch- 
sprachigen Literatur. Es gibt auf etwa 200 Seiten 
einen recht ausführlichen doch knapp gefaßten Über- 
blick über die in der Hochfrequenztechnik häu- 
figer benutzten mathematischen Hilfsmittel. Es dürfte 
dem in der Praxis stehenden Ingenieur oder Physiker 
ein wertvolles Nachschlagewerk werden. Diesem Zweck 
entsprechend werden Beweise kurz und nicht mit voll- 
ster Strenge durchgeführt, oft auch ganz weggelassen. 

Die ersten Kapitel befassen sich mit Vektoranalysis, 
Matrizenkalkül, Funktionentheorie, -Fourierschen 
Reihen und Integralen und mit der Laplace-Transfor- 
mation. Eingefügte Beispiele (u. a. Vierpole, Modula- 
tion, Übertragungsfunktion) zeigen die Anwendbarkeit 
der angeführten Methoden. Einem Kapitel über die 
Grundbegriffe der Statistik mit Behandlung des Wider- 
standsrauschens und des Schroteffektes folgt die Be- 
trachtung linearer Differentialgleichungen 2. Ordnung 
und zugehöriger spezieller Funktionen. Ein Kapitel 
über genäherte Lösung von Randwertproblemen bei 
gewöhnlichen Differentialgleichungen beschließt den 
ersten Teil des Buches. 

Im zweiten Teil kommen einige spezielle Probleme 
der Hochfrequenztechnik zu Wort. Behandelt werden 
Hohlraumresonatoren, Wellenleiter, Antennenprobleme 
und Elektronenströmung. Eine Vollständigkeit wird 
dabei nicht angestrebt, da es dem Verfasser darauf 
ankommt, an einigen speziellen Problemen der mo- 
dernen Forschung die Verwendbarkeit der behandel- 
ten Methoden zu illustrieren. 

Das umfangreiche Literaturverzeichnis von über 
50 Büchern und 70 Originalarbeiten verdient besonders 
erwähnt zu werden. 

Sowohl für den Praktiker als auch für den an den 
betreffenden Problemen interessierten Studenten ist 
mit diesem Buche ein sehr zu empfehlendes Werk ent- 
standen. 


Dresden. K.-H. Bachmann. 


Johannes Picht, Grundlagen der geometrisch- 
optischen Abbildung (Hochschulbücher für Physik, 
Bd. 14). VIII + 187 S. m. 74 Abb. Berlin 1955. VEB 
Deutscher Verlag der Wissenschaften. Preis geb. 
25,30 DM. 

Das konzentriert geschriebene Pichtsche Buch er- 
schien als Band 14 der „Hochschulbücher für Physik“. 
In der Tat darf es für sich in Anspruch nehmen, als 
Hochschullehrbuch bezeichnet zu werden. Der Autor, 
langjährig in der optischen Industrie tätig gewesen und 
hinübergewechselt zur Lehrtätigkeit an die Pädago- 
gische Hochschule Potsdam, durch fortbestehende enge 
Beziehungen zur optischen Industrie mit den modern- 
sten Problemen sowohl der Wellen- als auch der geome- 
trischen Optik wohl vertraut, hat aufbauend auf dem 
Stoff seiner Vorlesungen ein Lehrbuch geschaffen, das, 
wohl knapp bemessen, doch den neuesten Stand auf 
dem Gebiete der geometrischen Optik soweit es die 
Abbildungsgesetze betrifft, vermittelt. Man findet die 
wichtigsten theoretischen Grundlagen der geome- 
trisch-optischen Abbildung angeführt, das Gebiet der 
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allem den auf dem optischen Gebi 
tätigen Wissenschaftler zu interessiere 
schließt mit einem interessant gestalteten Abse 
über die Eikonale, die von Bruns in die Theorie 
geführten Verallgemeinerungen der optischen Ha 


ton-Funktion, und ihre Beziehungen zu den Abbil- 


dungsgesetzen. Man wird das Buch gern in den Hän- 


den jener Studierenden sehen, die sich dem Studium 


der Physik mit stark optischer Spezialisierung widmen, 
denn die Darstellung des Stoffes ohne Ballast garantiert 
die Aneignung gediegenen Wissens. Auch dem Dozen- 
ten der geometrischen Optik dürfte das Pichtsche 
Buch als Leitfaden für eine entsprechende Vorlesung 
empfohlen werden»können. Lediglich der Normenbe- 
flissene hätte zu kritisieren, daß die in dem hierfür in 


Frage kommenden Normblatt festgelegten Begriffs- 


unterscheidungen zwischen Maßstabsgrößen und Ver- 
größerungen konsequent hätten beachtet werden 
müssen. 


Jena. \ P. Görlich. 


Prof. Dr. Georg Scheffers }, Wie findet und 
zeichnet man Gradnetze von Land- und 
Sternkarten? Zweite Aufl. verbessert und erweitert 
von Prof. Dr. Karl Strubecker, Karlruhe. (Ma- 
thematisch-Physikalische Bibliothek, Reihel, Band 
85/86). IV + 1148. m. 30 Abb. u. 12 Tafeln. Stutt- 
gart 1956. B. G. Teubner Verlagsgesellschaft. Preis 
kart. 4,90 DM. 

Die erste Auflage dieses schönen auch für Nicht- 
mathematiker lesbaren Bändchens wurde in dieser Zeit- 
schrift Bd. 14 (1934), S.318 von E. Trefftz bespro- 
chen. Die Neuauflage ist im wesentlichen eine mecha- 
nische Wiedergabe der Darstellung der ersten Auflage 
mit den instruktiven Kartenentwürfen verschiedener 
Art (flächentreu, winkeltreu, perspektivisch usw.). 
Der neue Herausgeber hat in einem Anhang mathema- 
tische Ergänzungen hinzugefügt. Durch diesen Zusatz 
gewinnt die Darstellung den Anschluß an die Entwick- 
lungen der Differentialgeometrie und wird damit für 
den Studenten der Mathematik eine interessante Er- 
gänzung dieses Gebietes. 


Dresden. Willers. 


K. Beyer, Die Statik im Stahlbetonbau. 
Ein Lehr- und Handbuch der Baustatik. Zweite, 
vollständig neubearbeitete Auflage. Zweiter berich- 
tigter Neudruck. XII + 804 S. m. 1372 Abbildungen 
im Text, zahlreichen Tabellen und Rechenvorschriften. 
Berlin/Göttingen/Heidelberg 1956. Springer-Verlag. 
Preis geb. 66,— DM. 

Von der Fachwelt wird allgemein begrüßt, daß nun- 
mehr das bekannte Statikwerk des ehemaligen Pro- 
fessors für Technische Mechanik und Statik der Bau- 
konstruktionen an der Technischen Hochschule Dres- 
den in einem Neudruck der zweiten Auflage mit einigen 
Berichtigungen herausgegeben wurde. Unter den 
wenigen Gesamtdarstellungen einer Statik der biege- 
steifen Stabwerke, Scheiben, Platten und Schalen im 
Sinne der linearisierten Theorie I. Ordnung zeichnet 
sich die vorliegende durch eine straffe Gliederung des 
Stoffes und gleichzeitig durch knappe Formulierung 
aus. Obwohl die zahlreichen bis zu Endergebnissen 
durchgerechneten Anwendungsbeispiele vorwiegend 
der Praxis des Stahlbetonbaues entnommen sind, er- 
scheint der Titel des Buches zu eng gefaßt, weil’die 
allgemeinen Entwicklungen auch für die Berechnung 
der Stahl- und Holzbauwerke Gültigkeit haben, wenn 


man insbesondere von der Theorie des Fachwerks ab- 
sieht. 
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Z. angew. Math. Mech. 
Bd. 37 Nr.1/2 Jan./Febr. 1957 


Nach einer Darstellung der grundlegenden Voraus- 
setzungen der Theorie, der Herleitung der Gleich- 
gewichtsansätze und energetischen Grundgleichungen 
wird das statisch bestimmte Stabwerk behandelt. Im 
Anschluß an eine allgemeine Herleitung der Zu- 
sammenhänge zwischen dem Spannungs- und Ver- 
formungszustand werden baupraktisch häufiger auf- 
tretende statisch unbestimmte Stabwerke mit Hilfe 
der Kraftgrößen bzw. auch Deformationsmethode 
untersucht. In dem letzten Abschnitt des Buches geht 
der Verfasser in knapper Form auf grundlegende 
Probleme der Flächentragwerke ein. 

In der vorliegenden Ausgabe sind allerdings die in 
letzter Zeit viel angewandten zeitsparenden iterativen 
Lösungsmethoden noch nicht enthalten. Ferner wurde 
die Kontinuumsstatik im Sinne des räumlichen Zu- 
sammenwirkens der einzelnen Tragelemente, wie sie 
heute aus wirtschaftlicher Notwendigkeit die Berech- 
nung des Tragwerkes bestimmt, nicht weiter ent- 
wickelt. 

Das Buch von Professor Beyer ist zweifellos mit zu 
den besten Standardwerken der klassischen Stab- 
'statik zu zählen und wird nach wie vor nicht nur als 
Handbuch sondern auch als wertvolles Lehrbuch seine 
Anerkennung finden. 

Es ist erfreulich, daß der Springer-Verlag den 
zweiten Neudruck nunmehr wieder in gewohnt ge- 
diegener Ausführung herausgebracht hat. 


Dresden. G. Bürgermeister. 


Rufus Oldenburger (Direktor of Research Wood- 
ward Governor Company), Frequency Response. 
XII + 3728. m. zahlreichen Abb. New York 1956. 
The Macmillan Company. Preis geb. 7,50 $. 

Die auf dem Frequency-Response Symposium im 
Dezember 1953 in New York gehaltenen Vorträge 
liegen nun, redigiert von Rufus Oldenburger und er- 
gänzt durch wesentliche Beiträge aus dem (amerika- 
nischen) In- und Ausland, in Buchform vor. 

Eine begrüßenswerte Neuerscheinung, die uns einen 
gründlichen Einblick in Entstehung, Entwicklung 
und Stand der Methode des Frequenzverhaltens gibt. 
Das Werk ist in neun Teile gegliedert. 

Der erste Teil behandelt, nach einigen historischen 
Ausführungen, sehr klar die Grundlagen und Be- 
zeichnungsweise der Methode. Im zweiten Teil wer- 
den rechnerische und experimentelle Hilfsmittel ge- 
zeigt. Der dritte Teil bringt vier Anwendungsbeispiele 
aus Elektrotechnik, Maschinen- und Flugzeugbau, 
während im vierten Teil die zu verfahrenstechnischen 
Regelungen gehörenden Beispiele zusammengefaßt 
sind. Die Berechnung des Zeitverhaltens aus dem 
Frequenzverhalten, die Zusammenhänge zwischen 
Zeit- und Frequenzverhalten, die Anwendung kom- 
plexer Frequenzen und Ortskurven werden im fünften 
Teil dargelegt. Zur Optimierungsfrage nimmt der 
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sechste Teil mit zwei Beiträgen Stellung. Den nicht- 
linearen Regelkreisen ist mit vier Artikeln der siebente 
Teil des Buches gewidmet. Teil acht bringt dem Leser 
zwei wesentliche Beiträge über die Frequenzmethoden 
bei der Berechnung von Abtastregelungen (sampling 
control). Der neunte Teil schließt mit einem Beitrag 
über statistische Methoden das nunmehr sehr gerun- 
dete Bild des Frequenzverhaltens von geregelten 
Systemen ab. 

Dem an der Regelungstechnik interessierten In- 
genieur und Physiker vermittelt das Buch einen fast 
vollständigen Überblick über den Stand der Frequenz- 
methoden zur Berechnung von geregelten Systemen 
zum Jahresende 1953. Da ihm auch noch zahlreiche 
Literaturangaben das weitere Studium erleichtern, 
wird er gern auf dieses Werk zurückgreifen. Aber auch 
dem jungen Mathematiker kann das Buch empfohlen 
werden, denn es zeigt ihm, wie weit und einflußreich 
die Mathematik in die Regelungstechnik eingedrungen 
ist, welches Arbeitsgebiet ein Mathematiker hier zu 
erwarten hat, und schließlich kann er auch leicht er- 
kennen, welch weites Forschungsfeld ihm, insbeson- 
dere auf den Gebieten der nichtlinearen und un- 
stetigen Regelungen, noch offen steht. 


Mittweida. N. Klamka. 


John Robert Stock, Die mathematischen Grund- 
lagen für die Organisation der elektrischen 
Rechenmaschine der Eidgenössischen Tech- 
nischen Hochschule. (Mitteilungen aus dem In- 
stitut für angewandte Mathematik an der Eidgenös- 
sischen Technischen Hochschule in Zürich, Nr. 6.) 
768. m. 20 Abb. Basel und Stuttgart 1956. Birkhäu- 
ser-Verlag. Preis brosch. 7,30 SFr/DM. 

Die vom Verfasser entworfenen logischen Schal- 
tungen für Rechenwerk und Leitwerk des Rechen- 
automaten ERMETH werden beschrieben. ERMETH 
ist eine im Dezimalsystem im Serienbetrieb mit einer 
Impulsfolgefrequenz von 32 kHz hauptsächlich mit 
Diodenschaltungen arbeitende Maschine. Im ein- 
zelnen werden zunächst Wortdarstellung und ver- 
wendete Befehle erläutert. Das Rechenwerk kann so- 
wohl 11-stellige Zahlen mit beweglichem Komma in 
halblogarithmischer Darstellung als auch 14-stellige 
Zahlen mit festem Komma verarbeiten. Jedes Be- 
fehlswort enthält 2 Befehle. Der Ablauf der einzelnen 
Operationen und die zugehörige Steuerung im Leit- 
werk werden genau beschrieben. Besonders instruktiv 
sind dabei die Prinzipschaltbilder des Anhanges und 
die in Tabellen gegebene Zerlegung der Operationen 
in Einzelschritte. 

Es ist sehr zu begrüßen, daß mit dieser Mitteilung 
weiteren Kreisen ein tiefergehender Einblick in die 
logische Struktur eines modernen Rechenautomaten 
ermöglicht wird. 


Dresden. K.-H. Bachmann. 
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182 Abb. Braunschweig 1956. Verlag Friedrich 
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G. Doetsch (o. Prof. an der Universität Freiburg 
i.B.), Anleitung zum praktischen Gebrauch 
der Laplace-Transformation. Mit einem Ta- 
bellenanhang von Dipl.-Math. R. Herschel. 198 S. m. 
12 Abb. München 1956. Verlag R. Oldenbourg. 
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NACHRICHTEN 


Max Schuler zum 75. Geburtstag 


Am 5.Februar dieses Jahres beging Professor 
Dr.-Ing. Max Schuler in Göttingen seinen 75. Ge- 
burtstag. Der’ im Jahre 1946 emeritierte Gelehrte ist 
der Fachwelt durch eine große Anzahl bedeutender 
Arbeiten bekannt geworden. Es sei nur an die Auf- 
klärung des Schlingerfehlers beim Ein-Kreisel-Kom- 
paß, an die Schaffung des Drei-Kreisel-Kompasses, 
die Konstruktion einer automatischen Schiffs- 
steuerung, die Entdeckung des ‚‚84-Minuten-Prin- 
zips“ und an die genaueste Pendeluhr der Welt, die 
‚„Schuler-Uhr‘‘ erinnert. 


Professor Schuler hat nach längerer, überaus er-' 


folgreicher Industrietätigkeit noch die Laufbahn 
eines Universitätslehrers eingeschlagen und neben 
der Lehrtätigkeit auch die Leitung und den Ausbau 
des Institutes für angewandte Mechanik an der Uni- 
versität Göttingen übernommen. Durch zahlreiche 
Veröffentlichungen, durch Vorträge und Diskussions- 
bemerkungen auf Tagungen hat er in erstaunlicher 
Vielseitigkeit immer wieder Anregungen gegeben und 


Ergebnisse mitgeteilt, die nicht nur von einem ziel- 
sicheren Blick für das Wesentliche eines Problems, 
sondern auch von einer beneidenswerten Einsicht in 
die praktischen Möglichkeiten in der Mechanik zeu- 
gen. Nach seiner Emeritierung hat sich der Jubilar 
der Herausgabe einer Reihe von Büchern gewidmet. 
Möge ihm weiterhin Kraft und Gesundheit vergönnt 
sein, seine noch immer unermüdliche Arbeit auf ver- 
schiedenen Teilgebieten der Mechanik fortzusetzen. 

Es sei an dieser Stelle auch auf die ausführlichere 
Würdigung der Verdienste Schulers im Band 32 
(1952), 8. 228 dieser Zeitschrift hingewiesen. 


Freiburg, Breisgau. K. Magnus. 


Freiburg i. Br.: Der ‚‚Deutsche Verein zur För- 
derung des mathematischen und naturwissenschaft- 
lichen Unterrichts“ hält seine 48. Hauptversammlung 
vom 14. bis 18. April 1957 in Freiburg i. Breisgau ab. 
Nähere Auskunft erteilt Oberstudienrat Fritz Raith f 
Freiburg i. Br., Bürgerwehrstr. 18. 
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‚ andernfalls wird das folgende Quartal noch geliefert. Veröffent- 
n der Deutsche Demokratischen Republik. Gesamtherstellung: 


„Thomas Müntzer‘ Bad Langensalza (V/12/6) (1). Printed in Germany. 


WILHELM CAUER 
Theorie der linearen Wechselstromschaltungen 


2. verbesserte Auflage 


= 


Aus dem Nachlaß ergänzt von WILHELM KLEIN und FRANZ M. PELZ 


> .1954. XXIV, 769 Seiten — 461 Abbildungen im Text — 1 Kunstdrucktafel — gr. 8° — Ganzleinen DM 48,— 


„Das Lebenswerk von Cauer geht von der Erkenntnis aus, daß Fortpflanzungsmaß und Wellenwiderstand 
eines Vierpols, als Wheatstonesche Brücke aufgebaut, unabhängig voneinander festsetzbar sind. Die 
daraus folgenden Probleme der Realisierbarkeit und der Äquivalenz von Vierpolen sowie die Annan 
der idealen mathematischen Forderungen durch realisierbare (positive) Funktionen sind neben der Dar- 


stellung des Verhaltens der zuerst von ihm und W. Brandt im Jahre 1934 angegebenen Frequenzweichen 


konstanten Betriebswiderstandes der Inhalt des Buches. Die neue Auflage unterscheidet sich von der 


ersten durch Zusätze und Erweiterungen. So finden auch wesentliche Beiträge anderer Autoren bis in die 


jüngste Zeit hinein ihre Würdigung. Die Herausgeber haben sich besondere Verdienste um die Anwendung 


der Theorie und die leichtere Lesbarkeit des Buches erworben. Es stellt wegen der Behandlung der Grund- 


frage nach allen Lösungen der Aufgabe, lineare Wechselstromschaltungen mit vorgegebenen Eigenschaften 


zu realisieren, die tiefste und vollständigste der bisher erschienenen Arbeiten dar. Der Charakter eines 


Standardwerkes ist auch der zweiten Auflage geblieben.‘ 


(Jahresbericht der Deutschen Mathematiker-Vereinigung) 
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